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Let  p be  an  open  bounded  domain  in  En  with  smooth  boundary  r. 
Consider  the  following  boundary-value  problem: 

yt(t,x)  = -A(x,3)y(t,x)  on  Q = (0,T)xfl 

y(o,x)  = yQ(x)  on  a 

y(t,x)  = g(t,x)  on  E = (0,T)xT 

where  g is  taken  in  Lp[0,T;L2( T) ] with  1 S p £ °°.  The  objectives  of 

this  study  are  threefold: 

(i)  to  construct  a fully  discrete  Galerkin  approximation  of  (1); 

(ii)  for  an  abstract  model  containing  problem  (1)  as  a special 
case,  to  derive  error  estimates  in  the  norm  of  Lp[0,T;L2(Q)]  for 
boundary  input  in  Lp[0  ,T;  L2(  r)  ] , where  1 £ p < °°  and 

(iii)  to  present  numerical  results  corroborating  the  theoretical 
rates  of  convergence. 
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The  distinctive  feature  of  this  work  resides  in  the  nature  of  the 


boundary  conditions,  which  are  assumed  to  be  in  Lp[0,T;L2(D]  for 
1 i p i ».  It  is  well  known  that  standard  methods  based  on  variational 
inequalities  require  a certain  smoothness  of  the  boundary  input  which  in 

V, 

the  case  of  Dirichlet  boundary  conditions  is  H ^(D.  It  is  therefore 
apparent  that  the  treatment  of  "rough"  Dirichlet  boundary  data  requires 
the  development  of  a new  approach. 

Such  a new  approach  is  based  on  semigroup  representation  of  the 
solution  to  problem  (1),  as  well  as  on  that  of  its  semidiscrete  and 
discrete  approximations,  and  on  the  approximation  of  singular 
integrals.  We  make  use  of  this  approach  to  derive  theoretical 
convergence  estimates  for  the  fully  discrete  approximation  yhk  to  the 
solution  y of  problem  (1).  These  estimates  are  shown  to  be 


(2) 


ly-yhkl 


£ c(h 


V; 


Lp[0,T;L2(n)] 


2 in  1 ♦ 
h 


k 1/4  in  l)|g| 


Lp[0,T;L2(f)] 


2 

for  1 £ p £ ®,  and  with  the  restriction  k £ ch  for  some  c>0.  These 
estimates  are  also  shown  to  be  optimal  in  the  case  p = 1 , <*>,  and  almost 
optimal  (modulo  In  ^ and  In  in  the  case  1 < p < °°. 

The  construction  of  the  fully  discrete  approximation  scheme  makes 
use  of  a finite  element  technique  for  space  discretization,  combined 
with  a basic  finite-difference  scheme  for  the  time  variable. 

The  resulting  algorithm  is  then  applied  to  the  canonical  case  of 
the  one-dimensional  heat  equation  and  the  results  of  a computer  program 
implementing  this  algorithm  confirm  the  theoretical  estimates  stated  in 
inequality  (2). 
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CHAPTER  I 
INTRODUCTION 


1.1  Notations  and  Basic  Definitions 


The  main  goal  of  this  work  is  the  derivation  of  theoretical 
estimates  for  the  error  associated  with  the  full  discretization  of  a 
parabolic  equation  with  nonsmooth  Dirichlet  boundary  conditions.  In  the 
canonical  case  of  the  heat  equation,  for  example,  one  would  consider  the 
discretization  of  the  following  boundary-value  problem: 


3t  (t,x) 

= Ay(t,x) 

on  Q = 

(0,T)xQ 

y(0,x)  = 

0 

on  n 

y(t,x)  = 

g(t,x) 

on  E = 

(O.T)xf, 

where  fl  is  an  open  domain  in  En  with  smooth  boundary  r,  and  where  g is 
arbitrarily  chosen  in  L^Cl). 

Before  exposing  the  problem  however,  we  need  to  introduce  the  basic 
spaces  and  operators  in  terms  of  which  our  problem  can  be  precisely 
stated  and  solved,  along  with  some  of  their  properties. 

Let  x = (x-|,...,xn)  and  dx  = dx-|...dxn.  We  denote  by  ( £2 ) the 
space  of  equivalence  classes  of  measurable  functions  u such  that 


-1- 


-2- 


(1.1)  |u|  = (J  |u(x)  |2dx)/2  < ® 

l2(«)  ft 

With  the  inner  product  associated  to  the  norm  (1.1),  i.e. 

(u,v)  = | u(x)v(x)dx, 

L2(n)  ft 

L2(P)  is  a Hilbert  space. 

Before  introducing  certain  Sobolev  spaces,  we  briefly  recall  the 
definition  of  distributions  on  ft.  We  define 

(ft)  = {<t>:  4>  is  infinitely  differentiable  on  ft  with  compact  support 
in  ft} 

Then 

T)  '(ft)  = dual  of  ^Xft) 

= space  of  distributions  on  ft. 

•9’(ft)  is  provided  with  the  weak  topology,  i.e., 

fn  -*•  f in  -9'(ft)  if  and  only  if  <f ^ , 4>>  -*■  <f,<}>>  V<t>€.Z?  (ft). 


Remark . Any  function  in  L2(ft)  (or  in  (ft)  for  that  matter) 
defines  a distribution  on  ft,  also  denoted  by  f. 


The  derivative  - — of  a distribution  is  defined  by 

OX  . 

1 


<ib*> 

i 


_ <f  ii_> 


V4>i)e  (ft) 
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and  is  also  a distribution.  Then  Daf  is  defined  bv  iteration,  where 


Daf  = , a = (a,,..., a ) a multi-index  and 

a,  a in 


|a|  = o1  + a2  + ...+  an. 


We  are  now  in  a position  to  define  some  Sobolev  spaces.  First,  let 
m > 1 , m an  integer.  Then  Hm(fl),  the  Sobolev  space  on  Q of  order  m is 
defined  by 

Hm(fl)  = {u:  DaueL2(fl),  Va,|a|  £ m} 

with  norm 

|u|m  * ( l |d“u|2  )V2. 

H “»  |a|Sm  L2(n) 

Hm(n)  is  a Hilbert  space  with  the  inner  product  associated  with  the  norm 
above,  i.e., 


(u,v) 


l 

I a I Sm 


(Dau,D“v) 


l2(q) 


The  spaces  Hm(T)  are  defined  in  a similar  way. 

The  space  is  defined  as  the  closure  of  2)  (ft)  in  the  norm 

of  In  general  We  will  clarify  this  point  later 

on  in  this  work. 

We  also  define  the  spaces  Lp[0,T;H]  for  1 < p < °°,  where  H is  a 
Hilbert  space,  as  the  space  of  (equivalence  classes  of)  functions  f, 
strongly  measurable  on  [0,T],  with  range  in  X and  such  that 
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T 

(J|f(t)lP  dt)1//p  < * for  1 £ p < ® 

0 

and 

ess  sup  | f ( t ) | < ® for  p = ®. 

0 < t < T H 

We  note  that  L2[0,T;H]  is  a Hilbert  space  with  the  associated  inner 
product 

T 

(f ,g)L2[0,T;H]  = J^f^t)»S(t))Hdt 


df  dmf 

If  f»  all  belong  to  Ln[0,T;H],  then  f belongs  by  definition 

dt  ..  m v 

dt 

to  Hp[0,T;H],  with  norm 


m 

if  i m • i 141 

Hp[0,T;H]  j =0  dtJ  Lp[0,T;H] 


Finally,  we  define  the  spaces  Hp,s(Q),  where  Q 
are  nonnegative  integers,  as 


(0,T)xQ,  r and  s 


Hp’S(Q) 


L CO,T;Hr(Q)]  n H®[0,T;L  (Q)] , 

r r ^ 


with  norm 


|up  = |u|2  + |u|2 

Hp ,S(Q)  Lp[0,T;Hr(fl)]  Hp[0,T;L2(Q)] . 


The  spaces  Hp’s(Z),  where  I = (0,T)xT  are  defined  in  the  same  way, 
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Essential  to  the  forthcoming  study  is  the  notion  of  intermediate  or 
interpolation  space,  of  which  we  recall  the  definition  [13,  Vol . I]. 

Let  X,Y  be  two  Hilbert  spaces  with 

(1.2)  XcY,  X dense  in  Y, 
with  continuous  injection,  i.e., 

(1 .3)  |x|Y  < c |x|x  Vx  e X. 

Let  (*,*)x  and  ( ■ , • )y  be  the  scalar  products  in  X and  Y,  respectively. 
Let  D(S)  be  the  set  of  u's  such  that  the  antilinear  form 

V + (u,v)  , V £ X 

is  continuous  in  the  topology  induced  by  Y,  i.e., 

| vn  | Y + 0 =>  (u,vn)x  + 0 Vu  £ D(S). 


Then 

( 1 .^ ) (u,v)j,  = (Su , V)y 


defines  S as  an  unbounded  operator  in  Y with  domain  D(S). 
We  note  that 


(i)  D(S)  is  dense  in  Y. 
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* 

In  fact,  since  X is  dense  in  Y,  then  so  is  Y'  in  X'  by  duality,  where 
the  duality  is  considered  with  respect  to  the  inner  product  in  X. 
Therefore,  identifying  X and  its  own  dual  X',  we  have 

Y'  c X c y, 

each  space  being  dense  in  the  next  larger  space.  Furthermore,  the  space 
D(S)  can  be  identified  by  its  definition  with  Y',  and  thus  our  claim  is 
established . 

(ii)  S is  selfadjoint. 

(iii)  S is  strictly  positive,  and  in  fact, 

(Sv , v)y  = |v|^  > c|v|y. 

// 

0 

The  powers  S of  S,  0e  E can  then  be  defined  by  using  the  spectral 
decomposition  of  selfadjoint  operators.  In  particular,  we  define  the 
operator  T:  Y + Y by 

(1.5)  T = S l72  . 

T is  positive  and  selfadjoint  in  Y,  with  domain  X.  This  can  be  readily 
seen,  as  we  deduce  from  (1.4)  and  (1.5),  that 
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(u,v)x  » (Tu,Tv)Y  Vu, v t x. 

We  remark  that  since  S depends  on  the  choice  of  the  scalar  products 

on  X'  and  Y,  T is  dependent  on  those  too,  and  therefore  is  not 

intrinsically  linked  to  X and  Y.  We  are  now  prepared  to  define  the 

intermediate  spaces  [X,Y]  . 

0 

Definition . Given  two  Hilbert  spaces  X and  Y satisfying 
assumptions  (1.2)  and  (1.3),  and  T defined  by  (T.5),  we  set 

[X,Y]  = D(T1-6)  for  0 < 9 < 1 

with  the  graph  norm,  i.e., 

c|»|*  * |t’-6u|^)V2. 

// 

It  can  be  shown  that  [X,Y]_  is  intrinsically  linked  to  X and  Y,  i.e., 

0 

that  it  is  not  dependent  on  the  particular  scalar  products  chosen  on  X 
and  Y. 

Remarks . 

(i)  from  the  properties  of  the  spectral  decomposition  we  have  that 
X is  dense  in  [X,Y]„. 

0 f 

(ii)  [X,Y]0  = X and  [X , YD-|  = Y. 


// 
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Using  the  method  just  described,  we  can  define  the  intermediate 

Hp ,S(Q)  for  r-  s 2 0. 

[L  (ft) , H1 (ft)]  for  0 i s < 1 , 

^ s 

and  similarly  for  other  nonnegative  values  of  s. 

Then 


oyauco  n \ w / , 

For  example, 


def 

HS(ft)  ^ 


def 

Hp’S(Q)  = Lp[0,T;HP(ft)]  Hp[0  , T ; L2 (ft ) ] for  r,s  > 0. 


Sobolev  spaces  with  negative  indices  are  defined  in  the  following 

way : 


H S(Q)  = (H®(ft))'  = dual  of  H®(ft),  s £ 0. 

If  we  set 


Hn.’o  = = closure  of  X)  (Q)  in  Hr,S(Q), 

p , u , o p 

then  we  can  define 


H'r’"S(Q)  = (H*/*  (Q))’, 

P P;o,o 


r,s  £ 0. 
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For  a bounded  operator  A:  X -*■  Y,  the  operator  norm  will  be  denoted 


or  if  there  is  no  ambiguity,  simply  |a|. 

The  notation  will  be  used  for  the  inner  product  in 

L2 (ft ) or  L2[0,T;L2(fi)] , the  context  resolving  any  ambiguity,  and 
similarly  <•,*>  will  be  used  for  the  inner  product  in 
L2(D  or  L2[0,T;L2(r)]. 

£ 

We  shall  sometimes  use  the  notation  Dx  in  place  of  7^— , and 

J j 

also  l.i!  in  place  of  | . j or  |.|  , whenever  the  context 

s HS(n)  HS(D 

resolves  any  ambiguity. 


Let  tl  be  a bounded  domain  in  Bn  with  smooth  boundary  r.  Let 
T > 0 be  fixed,  and  consider  the  following  parabolic  problem  with 
nonhomogeneous  Dirichlet  boundary  conditions: 


| A | = sup 

X>Y  |x|  S 1 


1.2  Orientation  and  Formulation  of  the  Problem 


(t,x)  = A(x,3)y(t,x) 


on  Q = (0,T)xQ 


(1.6)  \ y(0,x)  = 0 


on  fl 


y(t,x)  = g ( t , x ) 


on  X = (0,T)xf 


where  g £L2(E)  is  given  and 


n 


n 


A(x,3)y  = I D (a  (x)D  y)  + I b (x)D  y + c(x)y 

1J  Xj  J-1  1 i 


x. 

1 


(1  .7) 
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is  a uniformly  strongly  elliptic,  second  order  differential  operator, 
i .e . , 

n n 

y a.  .£.£.£  Y y £7  for  some  Y > 0 
ij  = 1 1J  1 J i = 1 1 

with  ajj  , bi  , c e C°°(fl) . 

We  point  out  that  the  matrix  of  coefficients  (a_)  is  not  assumed  to  be 
symmetric,  which  is  significant  in  regard  to  some  of  the  literature. 

(See  section  1.5  on  the  literature). 

This  problem  is  a special  case  of  a general  abstract  model  which  is 
the  framework  of  our  study.  In  this  general  context,  our  aim  is  to 

(i)  construct  a fully  discrete  Euler-Galerkin  approximation  yhk  of 
the  solution  y to  the  problem  (1.6)  (h  and  k are  the  parameters  of 
discretization  in  space  and  time  respectively). 

(ii)  derive  optimal  error  estimates  (with  respect  to  the  regularity 
of  the  solution  y)  for  y - yhk  in  Lp[0 ,T; L2(Q ) ] norms  for  problems  with 
"rough"  Lp[0 , T; L,2 ( D 3 boundary  data. 

(iii)  present  numerical  results  from  a computer  program  which 
implements  the  aforementioned  approximation  scheme  in  the  canonical 
case  A(x,3)  = A (Laplacian).  These  results  will  be  seen  to  confirm  our 


theoretical  estimates. 
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The  first  step  in  the  construction  of  this  scheme  is  to  select  the 
basis  functions  of  a finite-dimensional  space  Vh,  which  is  a subspace  of 
the  apace  of  solutions  to  (1.6).  One  can  then  formulate  a continuous- 
time Galerkin  approximation  of  problem  (1.6),  the  solution  of  which  lies 
in  Vh  for  each  t. 

The  next  step  is  time  discretization  for  which  we  will  use  a basic 
finite  difference  scheme.  This  scheme  produces  a sequence  of 
approximating  functions  in  Vh  which  approximates  the  solution  at 
discrete  time  values.  More  precisely,  one  obtains  a sequence  of  linear 
systems  of  algebraic  equations,  the  solutions  of  which  are  the  desired 
approximating  functions.  The  task  of  solving  these  linear  systems  is 
conveniently  handled  by  computers,  which  can  produce  quick  and  accurate 
results . 

At  this  point  we  emphasize  that  the  low  regularity  of  the  boundary 
data,  L2(I)  or  more  generally  Lp[0,T;L2(r)]  as  opposed  to  H /2  (£) 
boundary  data  as  in  certain  known  treatments  in  the  literature  [ ■U , 5 ] (or 
spaces  with  higher  regularity),  is  the  main  source  of  difficulty.  In 
fact,  standard  approximation  methods  [2,  4,  5,  8,  21]  are  based  on  the 
projection  of  the  variational  formulation  of  the  solution  onto  finite- 
dimensional subspaces.  As  we  shall  see  later,  the  variational 
formulation  of  the  problem  (1.6)  requires  a higher  degree  of  smoothness 
of  the  boundary  data  than  is  available.  Thus  an  approach  different  from 
the  variational  one  has  to  be  developed. 
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The  techniques  that  we  propose  to  deal  with  LpLO.TjL^C  r)]  boundary 
data  are  based  on  semigroup  representation  of  the  solution  and  its 
approximations.  Their  presentation  is  relegated  to  the  next  chapter. 

We  note  here  that  this  approach  is  applicable  to  a variety  of 
situations,  for  example  Neumann  boundary  conditions  with  boundary  data 
in  the  space  of  distributions  (HS(r))',  s Z 0 , or  also  the  case  of  a 
differential  operator  of  order  2m,  m Z 1 , m an  integer.  Also  of 
significance  is  the  fact  that  no  restrictions  are  placed  on  the 
nonself adjointness  of  the  operator  A(x,3).  The  only  assumption  made  is 
that  A(x,3)  with  the  proper  domain  be  the  generator  of  an  analytic 
semigroup . 

The  motivation  for  using  the  semigroup  approach  as  opposed  to 
standard  variational  arguments  comes  precisely,  as  mentioned  before, 
from  the  low  regularity  of  the  boundary  input  g.  In  fact,  we  shall 
explain  below  why  the  standard  variational  approach  is  not  adequate  to 
treat  problems  with  "nonsmooth"  boundary  data. 

A standard  way  of  treating  the  nonhomogeneous  Dirichlet  boundary 
conditions  is  to  subtract  the  effect  of  the  boundary  data  from  the 
variable  and  to  consider  the  corresponding  nonhomogeneous  equation  with 
homogeneous  boundary  conditions.  To  this  end,  we  let  g be  an  extension 
of  g into  the  interior  of  fi  and  we  define  the  new  variable  y = y-g. 

Then 


(1.6')  y(0,x)  = - g(0,x) 


yt  = A(x,3)y  + A (x , 3 )g  - gfc 


on  Q 


on  a 


on  E 
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is  the  new  version  of  problem  (1  .6)  with  homogeneous  Dirichlet  boundary 
conditions . 

Let  f = A(x,9)g  - g . In  weak  form,  problem  (1.6')  becomes 

!(y  ,v)  = - a(y,v)  + (f,v)  , . for  y,ve  h'(Q) 

r L2(n)  l2(jj)  0 

y(0,x)  = - g(0,x)  on  ft 

where  a(u,v)  = - (A(x,9)u,v)  is  the  bilinear  form  associated  to 

L~(Q) 

A (x , 9 ) . 

If  y(t),  v(t)  in  (1.8)  are  restricted  to  lie  in  some  finite- 
dimensional subspace  of  H^(ft),  then  (1.8)  defines  a Galerkin 
approximation  of  the  solution  y.  The  stability  of  the  resulting 
approximation  is  proved  by  using  the  following  standard  variational 
argument  based  on  the  coercivity  of  the  bilinear  form  a(u,v).  Let 
v = y . Then 

(yt.y)  = - a(y,y)  + (f,y) 

L2(n)  l2(q) 

Since  (y^.y)  = ^ ^ ^ |y  1 2 ^ » by  integrating  in  t we  have 

T T T 

^2/  (|y|2)dt  = J - a(y(t) ,y(t))dt  + J (f (t),y(t))dt 


and  thus 
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T T 

V2  C |(y  (T ) | 2 -|y(0)|"  ) + Ja  (y  (t ) ,y  ( t ) )dt  = J ( f (t)  ,y(t)  )dt. 

l2(o)  l2(q)  0 0 


By  H1 -coerci vity  of  the  bilinear  form  a(u,v),  i.e.: 


2 2 1 

a(u,u)  £ c | u | - a|uj  Vu  e H (a)  for  some  c,a>0 

H (0)  l2(q) 


and  by  duality  (generalized  Schwartz  inequality),  we  have 


(1  .9) 


V2  |y (T ) j2  + c J|y(t)f2  dt  - a j |y(t)  J2 

L2(n)  0 H (a)  0 L2(Q) 


dt 


J|f(t)|2_1  |y(t)|21  dt  + V2  Ii(0) I 


0 h (a)  h (n) 


L2(n) 


From  the  identity 


1 2 2 

2 ab  < - a + E b , 
e 


e > 0 


we  obtain 


T T 

i2  _ ^ „2 


V2  |y (T)  I + cj|y(t)  |2  dt<cj|y(t)|2  dt  + c-1  J|f(t)|2  t dt 
T n u ""  0 1 ""  e 0 h (a) 


L2(a)  0 h ' (a)  0 l2(b) 


J|y(t)|21  dt  1-  V2  |5(o) | 


0 h (a) 


L2(a) 


+ c*e 
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Let  0 < e < 1 . Then 


T 


(1  .10) 


^2  |y ( T ) 1 2 + c(1-e)  J I y ( t ) | 2 dt 

l2(o)  o h (n ) 


T 


o h cn> 


dt  + V2  |g(0) j2 


In  order  to  conclude  that  the  left-hand  side  of  (1  .10)  is  well-defined, 
we  need  to  assume  that  f e L2[0,T;H_1 (n ) ] . As  f = A(x,3)g  - g , it 


theorem  [13,  Vol.  I,  p.  39]  enables  us  to  conclude  that  g(t)  must  be 


of  the  bilinear  form  a(u,v)  is  inadequate  for  treating  boundary  input  in 
^p[0,T;L2( T)] , anc*  therefore  a new  approach  is  necessary.  As  we  have 
indicated  before,  we  will  successfully  use  the  semigroup  approach 
combined  with  duality  theory  to  treat  this  type  of  nonsmooth  boundary 
conditions . 

Our  techniques  are  based  on  the  approximation  of  a semigroup 
formula  and  on  results  from  the  theory  of  singular  integrals  [20] 
applied  to  the  following  abstract  initial-value  problem  on  a Hilbert 
space  H with  input  from  a Hilbert  space  U: 


follows  that  we  need  to  assume  that  ge  L^tO , T; H (Q) ] . A basic  trace 


taken  at  least  in  H1/2(p). 


It  is  now  clear  that  this  standard  approach  based  on  H 1 -coercivity 
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r u = - Au  + Bf  on  D(A* ) ' (t) 

(1.11)  ] 

L u(0 ) = U 

0 

where  f e Lp[0,T;U],  B:  U D(A*)'  is  an  operator  such  that  D(B)  is 
possibly  empty  but  A- 'Be  £ (U-»-  H),  and  -A  is  the  infinitesimal  generator 
of  a strongly  continuous  analytic  semigroup  S(t).  The  solution  u 
--understood  as  a weak  solution--tnen  satisfies  the  formal  semigroup 
formula 

t 

(1.12)  u(t)  = S(t)uQ  + j S(t-i)Bf (-t)dt. 

0 


We  point  out  that  the  kernel  S(t-t)B  inside  the  integral  in  (1.12)  is 
singular  and  behaves  "comparably"  to  the  kernel  of  the  Hilbert 
transform.  In  fact, 


|S(t— t)B|  = |AS(t-x)A  ^1  < c | AS (t~T ) | 

U+H  U-»H  H+H 

since  A_1Be^  (U-»-  H).  From  the  analyticity  of  S(t),  it  follows  that 


| AS ( t-x ) | 

H+H 


< 


c 

t~T 


Therefore  the  integral  in  (1.12)  is  to  be  considered  as  a singular 
integral  defined  from  U into  H.  In  Section  2.2  it  will  be  shown  that 


(t)The  solution  of  (1.11)  should  be  understood  as  a weak  solution,  i.e., 
in  the  sense  of  distributions: 

(ut,<f>)H  = -(Au,4>)h  + (Bf ,<|))H  V<f>e  D(A* ) 
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(1  .12)  defines  an  element  u(t)  e h for  almost  every  tc  tO,T],  i.e.,  a 
strong  solution. 

In  Lasiecka  [11],  the  continuous-time  approximation  uh  is  achieved 
by  using  Galerkin  approximations  of  tne  above  singular  integral.  Our 
goal  in  this  work  is  to  extend  the  techniques  of  Lasiecka  [11,  12]  to 
the  discrete-time  approximation  uk  and  to  the  fully  discrete 
approximation  u^k  (with  discrete  semigroups).  By  doing  so  we  will  be 
able  to  prove  stability  and  the  optimal  rates  of  convergence  of  the 
resulting  fully  discrete  scheme  (see  Theorems  5.1,  5.2). 

This  general  abstract  approach  will  then  be  applied  to  the 
following  canonical  parabolic  boundary-value  proolem 

fy  = A(x,3)y  + f on  Q 

y(0,x)  = yQ(x)  on  n 

y I r = g on  I 

where  ge  Lp[0tT;L2(D],  yoeH~1/2(fl),  f e h"3/2,'3/4(Q)  , and  A(x,3)  is  a 
uniformly  strongly  elliptic,  second  order  differential  operator  given  by 
(1.7).  The  application  of  abstract  convergence  results  will  then 
produce  corresponding  convergence  results  (with  the  rates  of 
convergence)  for  the  parabolic  problem  (1.13)  (see  Theorems  3.3,  4.4, 
5.3,  5.4). 


1.3  Regularity  of  the  Solution 


In  order  to  appreciate  results  on  the  rates  of  convergence  of 
numerical  approximations,  it  is  necessary  to  understand  well  the 
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behavior  of  the  original  continuous  solution  and  in  particular  the 
optimal  regularity  of  the  solution.  It  is  the  purpose  of  this  section 
to  present  some  of  the  regularity  results  which  are  pertinent  to  our 
further  analysis. 

We  shall  start  by  recalling  the  definition  of  a strong  solution  to 
0 .13). 

Def inition.  A function  y e L^CQ)  is  said  to  be  a strong  solution  of 
(1.13)  if  there  is  a sequence  {y^}  in  i?(Q)  such  that 

|y  - y|  +o 

l2(q) 

|y  - s|  -o 

l2(D 

|D  y - A(x,3)y  - f|  - 0 

t n n l2(q) 

|y  (0)  - y | +0  as  n -*■  ®. 

n O . / \ 

L (n) 

// 

The  following  basic  result  on  well-posedness  and  regularity  for 
(1.13)  is  recalled  from  Lions-Magenes  [13,  Vol.  II,  p.  78]. 

Theorem  1.1  [13].  For  any  given 

- 1/0  -3/2  -3/4 

geL^E),  yQ  e H ^(fl),  f eH  ’ (Q),  there  is  a unique  strong 

solution  y to  problem  (1.13).  Moreover, 
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|y| 


H V2  ’ (Q) 


* c|g| 


l2(i) 


+ Cly0l  - Vo 

H 2 


(n) 


c|f|  . 


3/2, -3/4 


(Q) 


// 


In  the  sequel  we  shall  need  a regularity  result  of  a similar  nature 
to  Theorem  1.1,  but  with  the  more  general  data  g£  Lp[0,T;L2( r)]  for 

1 £ p £ 

Below  we  provide  the  statement  of  our  regularity  result  in  the 
context  of  problem  (1.6). 


Theorem  1.2.  Given  any  g e Lp[0,T;L2(  T)  ] for  1 < p < =>,  there  is  a 
unique  strong  solution  y to  problem  (1.6).  Moreover, 


v,  * c |g  I 

4 (Q)  L [0,T;L2(D] 


f or  1 < p < °° 


and 


lyl  V2-2e,  Vj,-E  * Clgl 


(Q) 


Lp[0,T;L2(f)] 


for  p = 1 , “,  and  e > 0. 


// 


The  assertions  of  Theorem  1.2  follow  immediately  from  Theorem  2.2  proved 
in  Chapter  II. 


1.4  Main  Results  for  the  Parabolic  Problem  (1.6) 


In  this  section,  we  shall  provide  for  orientation  purposes  only 
qualitative  statements  of  the  main  results,  while  the  precise 
formulation  with  all  the  involved  assumptions  is  relegated  to  chapter  V. 
Also,  for  the  sake  of  concretness,  we  shall  state  the  main  results  in 
the  context  of  the  parabolic  problem  (1.6). 


-20- 


The  main  results  in  this  work  are  of  two  kinds:  those  relating  to 

the  discrete-time  approximation  and  those  relating  to  the  fully  discrete 
approximation . 

Before  stating  our  results,  we  shall  need  to  introduce  some 
preliminary  material. 

First,  we  define  the  abstract  operator  A by 
D(A)  = {u:  A(x,3)u € 1^(8) , u | r ^ 0} 

and 

Au  = -A(x,3)u  for  ueD(A), 

where  A(x,3)  is  a uniformly  strongly  elliptic,  second-order  differential 
operator  given  by  (1.7).  It  is  well-known  [6]  that  -A  is  the 

infinitesimal  generator  of  a strongly  continuous  analytic  semigroup 

S(t). 

We  also  define  the  so-called  "Dirichlet"  map  from  the  boundary  r 
into  the  interior  of  a by 

Dg  = v if  and  only  if  A(x,3)  v = 0 

and  v | r = g . 


It  is  well-known  [14]  that 

(0)) 


D e £ (HS(f)  - HS+  1/2 


for  all  real  s. 
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The  time-discrete  approximation  yk  is  based  on  a backward  Euler 
difference  scheme  and  in  the  case  of  problem  (1.6)  has  the  formula  (as 
we  shall  prove  in  Chapter  IV) 


(1.14) 


Ji  - 

yk(nk)  = k l AE^+  JD(Pkg)(jk)  n = 1 N 

J-1 


t yk(t)  " yk(nk) 


for  (n- 1 ) k < t S nk 


where  Ek  = 1 R(l,  -A),  Pk  is  an  orthogonal  projection  on  Lp[0,T;L2(fl) ] 

or  Lp[0 ,T; L2 ( T) ] (we  shall  define  Pk  precisely  in  Chapter  IV),  and  g is 

the  boundary  input.  Here  and  thereafter  R(A,-A)  stands  for  the 

resolvent  of  A.  We  point  out  that  AEk  is  well  defined  since 

AR(A,-A)  e Jl  (L^Cft)  + L2(q))  (see  section  4.3  for  a more  precise  and  more 

general  statement,  as  well  as  the  necessary  justification). 

The  fully  discrete  approximation  yhk  is  based  on  the  same  backward 

Euler  difference  scheme,  with  the  difference  that  the  operator 
11  11 

k R^k’  cflan8e<3  to  E)lk  = — R(-,  -A^),  wliere  h parameter  of 

space  discretization,  and  where  A^  is,  for  example,  a Galerkin 
approximation  of  the  operator  A.  Ah  is  defined  on  a family  of  finite- 
dimensional subspaces  V^,  0 < h < 1,  with  suitably  chosen  approximation 
properties  among  which  in  particular  the  property  that 

# Jt 

D A v^tL^d’)  Vvh  £ Vh  (we  shall  see  later  that  this  is  equivalent  to  the 
property  that  ^ vh  6 *"2^  VvheVh’  whicfl  is  a natural  assumption  for 
splines).  The  projection  operator  Ph:  L2(n)  -*■  Vh  is  also  needed,  where 
Ph  is  orthogonal  on  L2(fl). 
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The  approximation  yhk  then  has  the  formula 


(1.15) 


yhk(nk)  = k l E[Jk1'jPh(AD(Pkg)(jk)),  n-1,...,N 

j-1 


yhk(t)  = yhk(nk) 


for  (n-1 )k  < t £ nk . 


Since  D(AD)  = 0 (the  domain  of  A only  contains  functions  which  vanish  on 
the  boundary),  the  term  Ph(AD(Pkg) ( jk ) ) requires  some  explanation.  In 
fact,  since  Ph  is  an  orthogonal  projection  on  we  have 

(Ph(AD(Pkg)(t)),vh)  = (AD(P  g)(t),v  ) VvheVh* 

l2(q)  l2(q) 

Then 

(1.16)  (AD(Pkg)(t),vh)  = <(P . g)(t),D*A*v.  > VvheVh. 

l2(q)  k h l2(D 

By  virtue  of  the  assumption  D*A*vh e L2( r) , the  right-hand  side  of  (1.16) 
is  well-defined.  This  shows  that  Ph(AD(Pkg) , (t) ) is  also  well- 
defined  . 


Remark.  We  shall  see  in  Chapter  V that  problem  (1.15)  is  in  fact  a 
sequence  of  linear  algebraic  system  of  equations  which  can  be  solved 
numerically  by  using  standard  subroutines. 

// 

We  are  now  in  a position  to  state  our  main  results,  and  we  start 
with  those  directly  related  to  the  discrete-time  approximation  scheme. 
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Lemma  1.3.  (See  Lemma  M.1  in  Chapter  IV  for  a more  general  and  precise 
statement ) . 

(a)  |(E^-S(nk))x|  < - |x|  n > 1 

l2(p)  n l2(q) 

(b)  |A(Ej-S(nk))x|  < |xj  n > 1 

l (n)  nK  l (n) 

// 

In  Huang-Thomee  [9],  this  result  is  proved  in  the  special  case  when 
a(u,v)  is  a V-elliptic  bilinear  symmetric  form,  where  a(u,v)  is  the 
bilinear  form  associated  to  the  operator  A,  and  defined  on  a Hilbert 
space  H with  V a subspace  of  H.  Using  techniques  different  than  those 
of  Huang-Thomee  [9],  we  snail  prove  Lemma  H.1  in  a more  general  setting 
without  assuming  that  A is  self adjoint.  Lemma  M.1  is  a key  result  in 
the  proof  of  Theorem  iJ.2. 

Theorem  1.4.  (See  Theorem  ^.2  for  a precise  statement).  Let  y be  the 
solution  of  problem  (1  .6)  and  let  yk  be  the  solution  of  its  discrete- 
time approximation  as  given  in  ( 1 .14).  Then 

|y  - ykl  ^ ck 1/4  “e  |g| 

Lp[0,T;L2(P)]  Lp[0,T;L2(r)] 

for  e > 0 given,  and  1 < p < °°. 

// 

One  can  note  that  for  1 < p < there  is  a loss  of  e in  the  rate 
of  convergence  for  |y  - yk|  as  compared  with  the  optimal  regularity  of 
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the  solution  y (see  section  1.3).  We  shall  see  that  this  loss  is 
partially  recovered  in  the  fully  discrete  estimate  (see  Theorem  5.3). 

In  fact,  for  the  fully  discrete  approximation,  we  have  the 
following  main  result: 


Theorem  1.5.  (See  Theorem  5.3  for  the  precise  statement).  Let  y be  the 
solution  of  problem  (1  .6)  and  let  yhk  be  the  solution  of  its  fully 
discrete  approximation  as  given  in  (1.15).  Then 

|y-y  | < c(k  /4lrJ-  + h /2ln^-)|g|  , 1<p<«. 

Lp[0,T;L2(fl)]  K n L [0,T;L2(D] 

2 

The  estimates  are  subject  to  the  restriction  k ^ ch  for  some  c > 0. 

// 

This  restriction  can  be  removed  by  imposing  stronger  approximation 
properties  on  the  subspaces  V^.  The  removal  of  this  restriction  also 
carries  the  price  of  a slight  loss  in  the  rates  of  convergence  for  k and 
h.  In  fact,  we  have  the  following  theorem. 


Theorem  1.6.  (See  Theorem  5.4  for  the  precise  statement).  As  in 

Theorem  5.1,  let  y be  the  solution  to  problem  (1.6)  and  let  yhK  be  its 

fully  discrete  approximation  as  given  in  (1.15).  With  the  additional 

2 

assumption  that  Vh  c_  H (fl),  we  have 


ly-yhkl 


c(k 


Vi 


4-e, 


Lp[0,T;L2(a)] 


hV2"2e)|g| 


Lp[0,T;L2(r)] 


1<p<= 


-25- 


without  any  restriction  on  k and  h. 


// 


Remark . The  spaces  Vh  nave  the  property  of  best  approximation, 
which  means  roughly  speaking  that  the  order  of  approximation  of  a given 
element  v by  elements  of  Vh  is  equal  to  the  number  of  derivatives  of 
v.  The  regularity  of  the  solution  y is  ^ C(Q)  for  p = 1,  ® 

V?  Vn 

and  Hp  ^ H (Q)  for  1 < p < In  view  of  the  best  approximation 
property  of  Vh,  the  result  of  Theorem  5.3  is  optimal  for  p = 1 , and 
almost  optimal  (modulo  In  1 and  In  for  1 < p < Similarly,  the 
result  of  Theorem  5.-4  is  almost  optimal  (modulo  e)  for  1 < p < «>  and 
essentially  optimal  for  p = 1 , 

// 


1 .5  Survey  of  the  Literature 


In  this  section  we  shall  collect  the  results  of  the  available 
literature  which  are  relevant  to  our  work  and  we  shall  compare  our 
results  to  those  found  there. 

In  our  work  we  consider  discrete  Euler-Galerkin  approximations  of 
linear  parabolic  boundary-value  problems  with  nonhomogeneous  and 
nonsmooth  Dirichlet  boundary  conditions.  There  is  a wide  range  of 
problems,  such  as  elliptic  and  parabolic  problems,  time-dependent 
problems,  linear  and  non-linear  problems,  problems  with  Dirichlet, 
Neumann  or  mixed  boundary  conditions,  continuous-time,  discrete-time  and 
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fully  discrete  approximation  schemes,  taken  in  several  combinations,  to 
which  certain  common  techniques  are  applicable. 

Two  basic  approaches  applicable  in  the  situations  described  above 

are 

(a)  energy  methods  using  variational  arguments  based  on  H1 -coercivity 
of  the  bilinear  form  associated  with  the  operator  a(x,3) 

(b)  a semigroup  approach  based  on  the  analyticity  of  the  semigroup 
generated  by  the  operator  A(x,3)  acting  on  some  space  with 
suitable  boundary  conditions. 

The  energy  methods  have  a wide  area  of  applicability,  including  for 
example  elliptic  and  parabolic  nonlinear  problems.  By  the  nature  of 
their  arguments,  as  seen,  for  example,  in  Douglas-Dupont  [5],  however, 
they  are  suited  for  treating  only  problems  with  high  regularity  of  the 

input  data. 

The  recently  developed  semigroup  approach,  applied  to  linear 
parabolic  and  hyperbolic  problems  (one  can  also  treat  some  nonlinear 
problems  by  considering  nonlinear  terms  as  perturbations),  permits  the 
treatment  of  input  data  with  significantly  lower  regularity 
requirements.  This  is  especially  significant  in  the  area  of  control 
theory  where  control  functions  are  in  LD[0,T;L2(r)]  as  a consequence  of 
the  natural  physical  assumptions. 

Below  we  shall  collect  and  review  some  results  on  approximation 
available  in  the  literature  and  we  shall  attempt  to  relate  these  results 
to  the  present  work. 

In  Douglas-Dupont  [5],  a number  of  problems  are  treated.  The 
starting  point  is  the  study  of  a nonlinear  parabolic  equation  with 
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homogeneous  Dirichlet  boundary  conditions.  Galerkin-type  approximation 
schemes  are  then  developed  for  both  the  continuous  case  and  the  discrete 
case.  Extensions  to  more  general  second  order  parabolic  equations  are 
considered,  such  as  the  time-dependent  case,  and  nonhomogeneous  boundary 
conditions.  The  points  relevant  to  our  work  are  that  the  second-order 
differential  operator  A is  assumed  to  be  selfadjoint,  and  that  the 
methods  used  are  variational  in  nature,  based  on  H1 -coercivity  of  the 
bilinear  form  associated  to  A.  Also,  the  convergence  estimates  involve 
-norms  of  both  the  solution  and  its  approximations.  The  solution  is 
assumed  to  have  all  of  its  derivatives,  with  respect  to  the  space 
variables  ana  the  time  variable,  in  L?( [O.Tjxfl) . This  implies  for 
example  that  the  boundary  data  must  be  at  least  in  H 1/2  (r)  for  the  case 
of  Dirichlet  boundary  conditions.  As  we  have  pointed  out  earlier,  these 
methods  are  not  applicable  to  the  case  of  boundary  data  in 
Lp[0,T;L2(T)]. 

In  the  paper  by  Bramble-Schatz-Thomee-Wahlbin  [4],  semidiscrete 
Galerkin-type  approximations  for  parabolic  equations  with  homogeneous 
Dirichlet  boundary  conditions  are  considered.  As  in  Douglas-Dupont , the 
operator  A is  assumed  to  be  selfadjoint.  L^-estimates  for  the 
semidiscrete  approximation,  including  estimates  for  L2~initial  data,  are 
derived  by  using  both  energy  inequalities  and  spectral  representation  of 
the  solution. 

As  a digression,  we  point  out  that  no  convergence  results  are  known 
as  yet  for  the  approximation  of  an  elliptic  equation  with  L^-boundary 


data. 
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The  extension  of  the  results  and  techniques  of  this  paper  [4]  to 
nonhomogeneous  Dirichlet  boundary  conditions  would  again  require  that 
the  boundary  data  be  in  H ^ (T) . 

The  paper  by  Huang-Thomee  [9]  studies  the  Dackward  Euler  method 
applied  to  a general  homogeneous  parabolic  initial-value  problem  in  weak 
form. 

Their  main  result  states  that 

|lln  - u(nk)  | < ^ |uq|  for  n = 1,...,N 

where  k is  the  parameter  of  discretization,  N,  uQ  are  given,  and  u(nk) 
is  the  actual  value  of  the  solution  at  time  t = nk  with  Un  its 
approximate  value.  As  we  shall  see,  the  statement  of  our  preliminary 

Lemma  4.1  is  essentially  the  same,  except  for  the  assumptions  made.  The 
authors  assume  the  bilinear  form  to  be  time-dependent  and  V-elliptic;  in 
our  paper,  although  time-dependent  problems  are  not  considered,  our 
assumption  that  -A  generates  an  analytic  semigroup  generalizes  the 
ellipticity  condition.  As  a result,  the  techniques  employed  for  the 
proof  of  our  Lemma  4.1,  based  on  an  integral  representation  of  the 
semigroup  involving  the  resolvent,  are  very  different  from  those  used  in 
Huang-Thomee  [9]. 

The  dissertation  by  Winther  [21]  treats  the  same  type  of  problem  as 
we  do,  i.e.,  the  fully  discrete  approximation  of  a parabolic  equation 
with  nonhomogeneous  boundary  data.  The  one  major  difference  is  that  he 
considers  Neumann  boundary  conditions,  as  opposed  to  Dirichlet  boundary 
conditions  in  our  case.  The  significance  of  this  difference  can  be  seen 
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by  realizing  that  a standard  variational  approach  requires  the  boundary 
term  to  be  in  H ^2  (r)  for  the  Neumann  problem  versus  H ^ ( r)  for  the 
Dirichlet  problem.  Thus  L2  Neumann  boundary  conditions  (which  are  the 
ones  treated  in  Winther)  are  tractable  with  variational  arguments. 
Winther  also  makes  use  of  spectral  representation  of  the  solution,  which 
requires  the  operator  A to  be  self adjoint. 

Now  we  turn  to  results  and  works  in  which  parabolic  problems  are 
treated  with  a semigroup  approach. 

In  Lasiecka  [10],  a general  and  unified  treatment  of  the  regularity 
of  solutions  of  parabolic  problems  for  a large  variety  of  nonhomogeneous 
boundary  conditions  is  provided.  The  treatment  is  abstract  in  nature 
and  the  only  assumption  made  on  the  operator  A is  that  -A  be  the 
generator  of  an  analytic  semigroup.  The  ensuing  theory  which  is 
developed  forms  the  basis  for  a number  of  results  in  the  area  of 
approximation  of  singular  integrals. 

In  fact,  semidiscrete  Galerkin  approximations  for  parabolic 
problems  with  homogeneous  boundary  conditions  have  been  considered  in 
Lasiecka  [11].  Although  this  particular  paper  is  not  directly  relevant 
to  our  work  (as  it  treats  only  homogeneous  boundary  conditions),  it 
establishes  several  results  which,  like  the  uniform  analyticity  of  the 
approximating  semigroup,  are  crucially  used  to  study  nonhomogeneous 
boundary-value  problems. 

As  a matter  of  fact,  in  Lasiecka  [12],  semidiscrete  Galerkin 
approximations  of  parabolic  problems  with  nonhomogeneous  boundary  data 
are  studied.  The  distinctive  feature  of  this  paper  is  that 
Lp[0,T;L2(r)]  (resp.  Lp[0,T;H_1 (r)])  boundary  conditions  are  treated  in 
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the  Dirichlet  case  (resp.  Neumann  case).  Using  duality  and 
interpolation  theory,  the  optimal  convergence  estimates  are  derived  for 
semidiscrete  Galerkin  approximations. 

These  results  [l^j  were  recently  used  by  Mackenroth  [15]  to  derive 
the  convergence  of  a numerical  algorithm  for  the  approximation  of 
abstract  Riccati  equations  arising  in  the  area  of  boundary  control 
problems.  These  results  [12]  will  also  be  used  in  our  work  to  prove  the 
stability  of  a fully  discrete  numerical  algorithm  which  approximates  the 
solution  of  problem  (1.6). 

To  summarize,  we  have  examined  two  different  approaches  for 
approximating  parabolic  boundary-value  problems.  First,  the  standard 
variational  approach  based  on  -coercivity  of  the  associated  bilinear 
form  is  characterized  by  high  rates  of  convergence,  but  on  the  other 

hand  requires  a high  degree  of  smoothness  of  the  boundary  data,  and  is 
limited  to  differential  operators  satisfying  some  ellipticity  condition. 

The  other  semigroup  approach  is  recent  and  has  been  developed  to 
treat  problems  where  the  boundary  data  is  "rough"  or  of  low 
regularity.  Also,  the  arguments  are  more  general  in  nature  from  the 
point  of  view  of  the  differential  operator  since  one  only  assumes  that 
-A  be  the  generator  of  an  analytic  semigroup. 

Our  work  focused  on  "nonsmooth"  boundary  data  is  a continuation  and 
extension  of  the  latter  semigroup  approach  [10].  The  main  results 
contained  in  Theorems  5.1,  5.2,  5.3,  5. A on  rates  of  convergence  of  the 
approximation  to  the  "nonsmooth"  abstract  equation  (1.11)  as  well  as  the 
"nonsmooth"  parabolic  equation  (1  .6)  are  new  to  the  author's  knowledge. 
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To  obtain  these  results  we  need  to  develop  new  techniques  based  on  the 
time  and  space  approximation  of  singular  integrals. 

1.6  Outline  of  the  Work  ' 

In  the  next  chapter  (chapter  II)  we  first  represent  the  solution  of 
our  abstract  problem  (1.11)  by  a semigroup  formula. 

We  then  recover  the  parabolic  problem  (1.6)  from  the  abstract 
approach.  At  the  end  of  chapter  II  regularity  results  for  the  solution 
to  (1.6)  are  proved  (Theorem  2.3). 

In  chapter  III,  we  first  provide  a general  treatment  for  the 
"space"  discretization  of  the  abstract  problem  (1.11).  We  then  provide 
basic  estimates  for  the  error  arising  from  the  space  approximation  of 
the  solution  to  (1.11)  (see  Theorem  3.1).  Later  we  specialize  these 
results  to  the  canonical  parabolic  problem  (1.6).  For  the  (space) 
discrete  approximation  uh(t)  defined  by  (3.25)  we  derive  the  basic  error 
estimates  in  Theorem  3.3. 

Chapter  IV  deals  with  discretization  in  time.  There  again,  we 
first  formulate  the  problem  in  an  abstract  setting.  Theorem  4.2  states 
an  abstract  convergence  result  and  later  in  section  4.5  we  specialize 
this  result  to  the  parabolic  case.  The  errors  for  the  time-discrete 
approximation  in  the  parabolic  case  are  estimated  in  Theorem  4.4. 

In  chapter  V we  treat  the  full  time  and  space  discretization. 

Error  estimates  for  the  full  discretization  of  the  abstract  problem 
(1.11)  are  given  in  Theorems  5.1  and  5.2.  By  specializing  these 
abstract  results  to  the  parabolic  case  we  then  derive  the  main  results 


-32- 


on  rates  of  convergence  of  the  fully  discrete  approximation  to  the 
parabolic  problem  (1.6),  stated  in  Theorems  5.3  and  5.4. 

Some  of  the  more  lengthy  and  tedious  proofs  of  results  stated  in 
this  work  are  relegated  to  Chapter  VI.  Thus  chapter  VI  consists  of  the 
proofs  of  Theorems  3.1,  4.2,  5.1,  5.2  as  well  as  Lemma  4.1. 

Finally,  in  chapter  VII  we  present  numerical  results  from  a 
computer  program  which  implements  the  full  discretization  scheme 
described  in  section  5.2  in  the  case  of  the  canonical  case  of  the  heat 
equation. 


CHAPTER  II 

SEMIGROUP  REPRESENTATION  FOR  PARABOLIC 
BOUNDARY-VALUE  PROBLEMS 


2.1  Opening  Remarks 


As  we  have  seen  earlier,  the  treatment  of  boundary  data  in 
Lp[0 ,T; L2 ( T) J is  not  covered  by  standard  methods  and  requires  a new 
approach.  Such  a new  approach  is  based  on  semigroup  representation  of 
the  solution,  and  on  the  theory  of  singular  integrals.  This  approach 
has  been  used  in  the  context  of  abstract  parabolic  problems  [3,  10]  and 
in  the  context  of  semidiscrete  Galerkin  approximations  [11,  12].  Our 
treatment  of  fully  discrete  Galerkin  approximations  uses  some  of  the  key 
ideas  presented  in  Lasiecka  [11,  12]. 

The  basic  stragegy  is  to  reformulate  problem  (1.6)  in  the  context 
of  the  abstract  model  (1  .11).  The  solution  of  this  abstract  model  has  a 
natural  semigroup  formulation  (see  (1.12)).  As  we  shall  see  later,  the 
continuous-time  Galerkin  approximation  of  this  abstract  problem  also  has 
a natural  semigroup  formulation,  the  properties  of  which  closely 
parallel  those  of  the  original  formulation. 

It  is  thus  crucial  to  understand  well  the  properties  of  the 
original  problem  in  order  to  be  able  to  derive  the  desired  estimates  for 
the  continuous-time  approximation.  A case  in  point  is  that  of  the 
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uniform  analyticity  of  the  approximating  semigroup,  which  reflects  the 
analyticity  of  the  original  semigroup. 

This  result,  seemingly  expected  and  straightforward , is  on  the 
contrary  highly  nontrivial  and  has  been  proved  in  papers  by  Huang  and 
Thomee  [8],  Luskin  and  Rannacher  [14],  and  Sammon  [19]  for  the 
selfadjoint  case  and  rather  recently  in  Lasiecka  [11]  for  the  general 
case. 

A similar  statement  can  be  made  in  the  case  of  the  discrete-time 
approximation.  In  fact,  the  operator  e£  (see  section  (1.5)) 
approximates  the  continuous  semigroup  S(t)  for  t = nk,  and  when 
considered  as  a discrete  semigroup  displays  properties  analogous  to  the 
analyticity  of  S(t).  These  properties  are  essential  to  the  proof  of 
Lemma  4.1  (which  states  the  approximation  properties  of  e£)  which  in 
turn  is  central  in  the  derivation  of  the  estimates  for  the  discrete-time 
approximation. 

For  simplicity  and  convenience  but  without  loss  of  generality  we 
will  concentrate  on  problem  (1.13)  with  f = 0 and  uq  = 0;  i.e.,  we  shall 
study  the  problem 


ut  = A(x,3)u 


on  Q 


(2.1) 


u(0,x)  = 0 


on  n 


on  I 


n 


n 


where  A(x,3)u  = 


l Dx  (an(x)D  u)  + l b (x )D 
i . j = 1 i 1J  Xj  i-1  1 x 


u + c(x)u 


= div  (M(x)Vu)  + b(x) *Vu  + c(x)u 
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is  a uniformly  strongly  elliptic,  second  order  differential  operator, 


Also,  as  we  have  pointed  out  in  section  1 .2,  we  do  not  assume  the  matrix 
M(x)  to  be  symmetric. 


In  this  section  we  introduce  an  abstract  framework  into  which 
problem  (2.1)  can  be  incorporated.  To  this  end  we  shall  introduce 
certain  abstract  operators. 

Let  A be  a densely  defined  closed  linear  operator  on  a Hilbert 
space  H (A:  D(A)  c H +H ) . We  assume  that  -A  is  the  infinitesimal 
generator  of  a strongly  continuous  analytic  semigroup  S(t).  It  is  well- 
known  [18]  that  equivalently, 

(2.2)  p(-A)2  l = {A:  ReA  > a - b j ImA | } for  some  a,b  > 0 


i .e. , 


bi  , c c C ( P ) • 


2.2  Semigroup  Formula  for  the  Solution  of  an  Abstract 
Boundary-Value  Problem 


and 


(2.3)  |R(A,-A)  | < — — 

H+H  1+| A | 


VA  el. 


Equivalently  to  (2.3), 


(2.4)  |AnS(t ) I £ — eat 

■ .n 


H+H  t 


, t > 0,  n 2 0. 
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We  will  assume  without  loss  of  generality  that  a < 0.  Then  0e  p(-A), 
which  ensures  that  the  fractional  powers  Aa,A*a,  0 < a < 1 are  well- 
defined  [ 1 8 ] . In  fact,  if  a > 0 it  is  enough  to  translate  the  generator 
to  the  left;  i.e.,  we  can  take  A = A + II  with  X > a,  and  so 

R ( A ,-A)  = (A+A)'1  = (A+A+A)”1  = R(A+I,-A) 

which  shows  that  p(-A)  2 l with  a < 0. 

*-a 

Also,  the  notation  D(A  ) which  will  be  used  several  times  in  the 

£ _ 

sequel  is  defined  by  D(A  a)  = D ( Aa) ' , 0 < a < 1. 

We  now  let  B:  U-»D(A*)’  be  a linear  operator  from  a Hilbert  space  U 
into  D(A* ) ' , with  the  condition 

(2.5)  A_1B  e £ (u  •+  H) . 

We  consider  the  following  abstract  model; 

( u (t)  = -Au(t)  + Bf ( t ) on  D ( A* ) ' 

(2.6)  ] 

( u(0 ) = uq 

where  uQ  is  given  in  H and  feLp[0,T;U],  1 < p < ». 

Remark . Since  D(B)  considered  on  U may  be  empty  (this  will  be  the 
case  for  our  application),  problem  (2.6)  should  be  understood  in  the 
sense  of  D(A*) ’-topology , i.e.,  as 
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( (u  (t),v(t))  = -(u(t) ,A*v(t) ) + (f (t)  ,B*v(t) ) 

L n ri  U 

(2.7)  < 

t u(0)  = uQ.  Vv(t)  t D(A*),t € (0 ,T] 

// 

The  analyticity  of  the  semigroup  S(t)  implies  (as  we  shall  see)  the 
existence  of  a solution  u(t)  defined  almost  everywhere  in  the  interval 
[0,T]  for  any  given  f e Lp[0,T;Uj. 

The  input-solution  operator  T for  the  problem  (2.6)  formally 
satisfies  the  semigroup  formula 

t 

(2.8)  u ( t ) = Tf(t)  = S(t)uQ  + J S(t-i)Bf(i)dT. 

0 

Remark . The  kernel  S(t--r)B  of  the  integral  in  (2.8)  has  the 

property  that 

j S ( t—  t ) B | = |AS(t-x)A_1B|  <c|AS(t-x)| 

U+H  U-*-H  H+H 

by  (2.5). 

Then  by  analyticity  of  the  semigroup  S(t), 

|AS( t-x) I < c • -i-  , 

H>H  T 

1/(t--r)  being  the  kernel  of  the  Hilbert  transform.  Therefore  the 
integral  in  (2.8)  should  be  considered  as  a singular  integral. 

// 


-38- 


At  this  point  we  again  assume  for  simplicity,  but  without  loss  of 
generality,  that  uQ  = 0. 

We  shall  now  see  that  the  formula  (2.8)  defines  almost  everywhere 
in  t an  element  in  H for  f e Lp[0,T;U].  If  we  assume  as  ir.  (2.5)  but 
more  generally  that 

(2.9)  A ^B  e .Z?  (U-»H ) for  some  0 £ q £ 1 , 

T can  be  rewritten  as 

t 

(2.10)  Tf (t)  = Aq  j S ( t- t ) A~QBf ( t )dt . 

0 

The  following  regularity  result  for  T is  proved  in  Lasiecka  [12],  where 

the  theory  of  singular  integrals  [20]  plays  a central  role  in  the 
arguments. 

Theorem  2.1  [12].  Let  0 < q < 1 and  assume  condition  (2.9).  Then 

(a)  Tel  ( [ Lp [ 0 ,T; U]  - Lp[0 ,T; D(A*1 '*) ] ) for  1 < p < » 

(b)  Te  £ (Lp[0,T;U]  - Lp[0  ,T;  DU*1  ~q"e)  ] ) for  p = 1 , «,  e>0. 

// 

We  shall  also  need 


Theorem  2.2.  Let  0 £ q < 1 and  assume  condition  (2.9).  Then  for  T 
satisfying  formula  (2.10),  we  have 
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-*■  Hp~Q[0,T;H])  for  1 < p < ® 

- Hp“q_e[0,T;H])  for  p = 1,  »,  e>0 . 

// 


S(t-r)Bf (t)dT. 
t 

A q Hf(t)  = f(t)  + A S S(t-l)f(T)dT 
at  ^ 

where  f(t)  = A qBf(t).  By  assumption  on  B and  f,  we  know 
that  f e L [0,T;H].  Let 

t 

Tf(t)  = A j S(t— T)f(t)dT. 

0 

a 

In  the  case  1 < p < ®,  by  applying  part  (a)  of  Theorem  2.1  to  T with 
B = I and  q = 1 , we  obtain 

(2.11)  A_q  ^ ef(Lp[0,T;U]  + Lp[0,T;H]),  1 < p < 

Similarly,  in  the  case  p = 1 , by  applying  part  (b)  of  Theorem  2.1 
to  T with  B = I and  q = 1-e,  we  obtain 


(a)  T cl  (L  [0,T;U] 


(d)  T € JL  (L  [0,T;U] 


Proof.  We  have 


^(t)  = Bf ( t ) + A J 


Then 


(2.12) 


A-q-e 


£ (L  [0,T;U] 


Lp[0,T;H] ) , p = 1, 
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We  now  apply  the  intermediate  derivatives  theorem  [13,  Vol . I,  p.  15] 
whose  original  formulation  in  the  setting  of  L2  spaces  can  be  extended 
to  Lp  spaces,  to  (2.11)  and  (2.12),  which  yields  (a)  and  (b)  as 
claimed . 

// 

These  regularity  results  will  be  used  to  derive  corresponding 
regularity  results  for  our  canonical  problem  (2.1). 

2.3  Application  to  a Canonical  Example  of  a Parabolic 
Boundary-Value  Problem 


The  goal  in  this  section  is  to  derive  a semigroup  formula  for  the 
solution  of  problem  (2.1).  This  will  be  achieved  by  expressing  it  in 
the  context  of  section  2.2.  Thus  we  set  H = L2(Q),  U = L2(D.  We  also 
need  to  introduce  a few  operators. 

Let  D:  L2(T)  -*■  L2(fl)  be  the  so-called  "Dirichlet"  map,  which  is  a 
natural  extension  into  the  interior  of  fl  of  a function  defined  on  the 
boundary  T.  By  definition, 


Dg  = v if  and  only  if  A(x,3)  v = 0 


and  v I r = g. 


It  is  well-known  [16]  that 


(2.13)  D e i(HS(T)  ^ HS+ 1/2  (Q)) 


for  all  real  s. 
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We  also  define  the  operator  A:  L^(n)  L^CQ)  by 
D(A)  = If:  A (x  , 3 )f e L^fl),  f | = 0} 

and 

Af  = -A (x , 3 )f  for  f e D(A) . 

It  is  well-known  [6]  that  -A  is  the  infinitesimal  generator  of  a 

strongly  continuous  analytic  semigroup  S(t)  on  L2 ( n ) . 

We  begin  the  derivation  of  the  desired  semigroup  formula  with  the 

following  standard  argument,  in  which  we  need  to  assume  that 
2 1 

geHp’  (Z)  (later  we  shall  extend  the  formula  to  all  of  Lp[0,T;L2(r)]  by 
density ) . 

Let  y = y-Dg.  Then,  similarly  as  in  section  1.2, 

yt  = yt  + (Dg)t  and  A(x,3)y  = A(x,3)y  + A(x,3)Dg 

= A(x,3)y 

since  by  definition  A(x,3)Dg  = 0. 

Problem  (2.1)  then  becomes: 

{yt  = A (x , 3 )y  -(Dg) 

y(0,x)  = -Dg(0,x) 

L ^Ir  " 0 

Since  y satisfies  zero  Dirichlet  boundary  conditions,  problem  (2.U)  can 
be  put  into  the  form  of  an  abstract  boundary-value  problem  by  defining 
the  following  "abstract"  variables: 


on  Q 

on  fl 
on  Z. 
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y(t)  = y(t,  •)  e L2(n) , 


y(t)  - y ( t , • ) e l2(ji). 


g(t)  = g( t , - ) e L2(D  for  t > 0. 


Then  the  abstract  problem 


yt  = -Ay  - D(gfc) 


on  Q 


(2.15) 


y(0)  = -Dg (0 ) 


on  n 


is  a special  case  of  problem  (2.6)  with  B = -D  and  u0  = -Dg(0).  We  note 


note  that  the  boundary  conditions  are  explicit  in  (2.14)  whereas  they 
are  implicit  in  (2.15).  (D(A)  only  contains  functions  which  vanish  on 

the  boundary ) . 

The  variation  of  parameter  formula  gives  the  following 
representation  for  the  solution  of  (2.15): 


t 

(2.16)  y(t)  = -S(t)Dg(0)  - / S(t-T)(Dg)  (x)dx. 

0 1 


Our  objective  is  to  obtain  a corresponding  formula  for  u(t). 
Integrating  by  parts,  which  is  justified  since  gfH^d),  gives 


t=0  0 

t 


= Dg(t)  - S ( t )Dg (0 ) - J AS(t-T)Dg(T)dT. 

0 


Since  y(t)  = y(t)  + Dg(t),  we  can  see  that 
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t 

(2.17)  y (t)  = A j S(t-t)Dg( t ) d t 

0 

is  the  desired  semigroup  formula  for  the  solution  of  problem  (2.1). 
Since 

(2.18)  AaS(t)x  = S(t)Aax  for  all  real  a,  and  x£  D(Aa), 


we  can  formally  rewrite  (2.17)  by  using  (2.18),  as 
t 

(2.19)  y (t)  = j S(t-x)ADg(t)dT. 

0 


Now  we  can  see  that  y is  the  solution  of  the  abstract  problem 


(2.20) 


fyt  = -Ay  + ADg 
y(0)  = o 


on  D(A*)' 


L 


of  which  the  canonical  example  is  of  course  (2.1).  Since  Dg | p = g * 0 
in  general,  it  follows  that  D(AD)  = 0.  It  is  therefore  important  to 
note  that  both  (2.19)  and  (2.20)  should  be  understood  (as  we  have 
pointed  out  in  the  remark  of  section  2.2)  in  the  sense  of  D ( A* ) ' - 
topology.  A proper  way  of  stating  (2.20)  would  be  for  example 

C (ut(t),v)fl  = (-Au(t),v)fl  + (ADg(t)  ,v ) 


(2.21  ) 


= -(u(t),A*v)^  + <g(t)  ,D*A*v>r  for  veD(A*) 


u(0 ) = 0. 


The  semigroup  formula  (2.17)  is  defined  so  far  for  g£H^’'(Z)  but 

P 

we  shall  show  in  the  next  section  that  it  can  be  extended  by  density  to 
g c Lp[0,T;L2(r)]. 

2.M  Regularity  of  the  Solution  to  (2.1) 

Let  L:  Lp[0 ,T; L2 ( T) ] ■*  Lp[0 , T ; (P ) ] be  the  input-solution  operator 
for  problem  (2.1),  i .e . , 

t 

(2.22)  Lg(t)  = A j S(t-T)Dg(x)dT. 

0 

L is  certainly  well-defined  on  H^^d)  as  we  have  seen  in  section  2.2, 

so  that  Hp’1 (Z)  £D(L) . Since  ’ 1 ( Z ) is  dense  in  Lp[0,T;L2(r)]  it 

follows  that  L is  densely  defined. 

We  shall  now  state  and  prove  our  regularity  result  for  the  solution 
of  problem  (2.1)  given  by  the  semigroup  formula  (2.22). 

Theorem  2.3. 

(a)  L ei(Lp[0,T;L2(r)]  ->  ^ (Q))  for  1 < p < » 

(b)  L (Lp[0,T;L2(  T)]  * Hp/2  2C’  ^ £(Q))  for  p = 1 , »,  e>0. 

// 


Proof.  Knowing  that  [7] 
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(2.23)  D(Aa)  = H2a(Q)  for  0 < a < V4  , 

we  have  by  (2.1 3) 

(2.24)  De  £ (L2(  r)  -*■  D(A  ^ G))  for  every  e > 0. 

This  fact  will  be  used  repeatedly  in  the  sequel. 

We  first  prove  the  "easier"  part  (b). 

The  semigroup  formula  (2.22)  can  be  rewritten  equivalently  as 

t 

(2.25)  y (t)  = Lg ( t ) = A3/4"£  j S(t-T)A_3/4_e(AD)g(T)dT  for  e>0. 

0 

Thus  we  are  exactly  in  the  situation  described  in  (2.6),  with 

H = L2(Q),  U = L2(D,  B = AD  (formally),  and  with  q = 3/4+e  in  (2.9). 
With  these  identifications,  it  follows  from  part  (b)  of  Theorem  2.1  and 
from  (2.22)  that 

(2.26)  Le  J!  (Lp[0,T;L2(r)]  + Lp[0  ,T;Hp1/2  ~2e(fl)  ] ) , e > 0. 

The  time  regularity  of  L follows  from  part  (b)  of  Theorem  2.2  with  the 
same  identifications.  Thus  we  have 

(2.27)  Le  i (Lp[0,T;L2(D]  + "£[0  , T;  L2  (n)  ] ) , e>0. 


This  completes  the  proof  of  (b),  as 
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H v2  ~2t,  Vi|  -e 

P 


(Q) 


= Lp[0,T;H 


V2  - 2e 


(n)] n h 


V4  -e 


[0,T;L2(Q)]. 


We  now  prove  (a).  This  will  be  achieved  by  interpolating  the  two 
following  results. 


(a1  ) Le  i (Hp3/2,3/4(Z)  - H2,1(Q))  1 < p < », 

(a2)  L e i (Hp”  1/2  1/J4  (Z)  ->  Lp[0,T;L2(P)])  1 < p < ® 

where  Hpr,S(I)  = {feHp,S(I):  f(0)  = yQ|p  = 0},  i.e.,  the  elements 

of  Hp  ’°(I)  satisfy  compatibility  relations. 

We  point  out  that 


(2.28)  Hpr’S(Z)  = Hp'S(I)  for  |r|  < V2  ,|s|  < % , 1 S p < 

Thus,  in  the  space  of  solutions  Hp/2  ’ (Q),  the  compatibility  relations 

also  disappear.  This  is  due  to  the  fact  that  H 1/2  (n)  = H1/2(a);  i.e., 

V? 

the  topology  of  Hq  ^ (p)  is  too  weak  to  "recognize"  boundary  values. 

First,  we  show  (a-|).  Because  of  the  high  regularity  of  the 

boundary  input,  we  shall  use  a standard  argument. 

Let  y = y _ g,  where  g extends  g from  1 into  the  interior  of  Q,  in 

such  a way  that  the  compatibility  relations  are  preserved;  i.e., 

g(0)  = yQ  = 0 on  £2.  By  inverse  trace  theory,  we  know  that  if 
3/2  3/4  - 2 1 

g e Hp  ’ (E),  then  g £Hp’  (Q)  [13,  Vol.  II,  p.9].  Therefore,  it  is 

sufficient  to  establish  the  same  regularity  for  u. 
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We  recall  from  section  1 .2  that  y satisifes  the  boundary-value 
problem 


(2.29) 


= -Ay  + f 


on  Q 


^ y (o)  = -g(o)  on  n 

where  f(t)  = A(x,3)g(t)  - g (t).  Since  geH2,1(Q),  it  follows  that 

t p 

f e Lp[0 ,T; L2(0) ] . Then  the  solution  of  (2.29)  is  given  by  the  semigroup 
formula 


t 

y(t)  = -S(t)g(0)  + J S(t-T)f (-r)di 

0 

t 

= j S(t-T)f(-r)dT 
0 


because  of  the  compatibility  relation  i(0)  = 0.  Hence 


t 

Ay (t)  = A / S(t-T)f(x)dT. 
0 


By  application  of  Theorem  2.1  with  U = H = L2(n),  B = A and  q = i , we 
get 


(2.30)  |y| 

Lp[0 , T; D ( A ) ] 


I Ay  I 


Lp[0,T;L2(P)] 


< c|f  J 


Lp[0,T;L2(fl)] 


1 <p<« 


Using  the  identification  [7] 
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D(A)  = D ( A*  ) = n H1 (Q) 

o 

and 

M po  ^ C lx  I for  x c D(A8)  , 0 < B < 1 , 

*T8(fl)  D(A8) 

we  can  see  that 


|y|  2 s c|r|  , i < P < «. 

Lp[0,T;H  (0)]  Lp[n,T;L2(fi)] 


As  for  the  time  regularity  of  u,  we  have 


t 

y (t)  = -Ay (t)  + f(t)  = -A  j S(t-T)f(T)dt  + f(t) 

0 


and  so  by  (2.30)  we  obtain 


|u  | <c|f| 

Lp[0,T;L2(n)] 


Lp[0,T;L2(jj)] 


1 < p < 


Therefore  (a1  ) is  proved  since  HpJ(Q)  = Lp[0,T;H2(n)  ] n Hp[0,T;L2(fi)] 

We  now  need  to  prove 


(a?)  Le  i(H  1/2  ’ 1/4  (I)  -> 


L_[0,T;Lo(fi)]) 


1 < p < ®. 


Assuming  the  validity  of  (ap),  our  result  follows  directly  from  the 
application  of  a general  interpolation  theorem  [13,  Vol . I,  p.  27], 
More  precisely,  the  theorem  concludes  that 


L e 


i([H3/2’3A  (I),  H -1* 


2 


Vi* 


(T)V 


[H2,1(Q),H°'°(Q)]J 
p P o 
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for  0 < 6 < 1 . 

Since  we  are  concerned  with  ge  Lp[0 ,T ; L2 ( r ) ] = H°’°(l),  we 
take  6 = 3/^.  By  our  remark  (2.28),  we  deduce  that 

Le  I (Lp[0,T;L2(r)]  + ^ (Q) ) 

as  desired. 

As  for  the  proof  of  ( ) we  shall  use  a simple  but  useful  duality 
result.  We  state  and  prove  this  result  in  an  abstract  context  as  it 
will  be  used  repeatedly  throughout  this  work.  First  of  all,  we  define 
the  adjoint  of  an  operator. 

Definition.  Let  B:  U ■*  H be  a densely  defined  linear  operator, 
where  U,  H are  arbitrary  Hilbert  spaces.  Then  the  adjoint  B*  of  B is 
defined  by  B*v  = z,  where 

(2.31)  (Bu,v)h  = (u,z)u  for  all  u e D(B) 

and  D(B*)  is  the  set  of  all  z such  that  condition  (2.29)  is  satisfied 
for  some  z e U . 

Now  we  have 


Lemma  2.H.  Let  K be  an  operator  such  that 
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K£  £ (L2[0,T;U]  - L2[0,T;H]), 

where  U,  H are  arbitrary  Hilbert  spaces  and  T > 0 is  fixed. 
Then 


Ke  £ C Lp L 0 ,T;U]  -+  Lp[0,T;Hj)  implies  that 


K*  e £ (Lq[0,T;H]  «•  Lq[0,T;U])( 


where  — + — 

p q 


1 for  1 < p < ® and  q 


“ (resp.  1)  for  p = 1 (resp.  ®) . 

// 


Proof . The  result  is  immediate  for  1 £ p < ® by  duality  and  by 
definition  of  the  adjoint.  We  now  prove  the  case  p = ®. 


First,  by  definition  of  K*, 

t t 

(2.32)  J(Ku(t),v(t))  dt  = J(u(t),K*v(t))  dt  VueD(K),  v£D(K»). 

0 M 0 U 

Let  v £ L-|  [0 , T;H] , and  let  { } be  a sequence  in  L2[0,T;H]  such 

that  vnMin  L-|[0,T;H].  Then  K*vn  e L2[0,T;U]  for  each  n by  Lemma  2.4 

with  p = 2.  With  u e LjO.TjU],  consider 

t t 

(2.33)  J(Ku(t),v  (t))  dt  = J(u(t),K*v  (t))  dt 

0 n H o n u 


Then  by  assumption 
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(2.34)  f (Ku ( t ) , v (t))  dt  S c lu  I | v I 

i n H n1 


LJO.TjU]  L1  [O.TjH] 


S c|u|  |v| 


We  now  let  u(t)  = 


K*vn(t) 

lKV>lu 


LjO,T;U]  L1 [0 , T;H ] 


which  is  well-defined  in  L [0,T;U]  as 


K*vn(t)eU  almost  everywhere  in  t.  Hence  by  (2.33)  and  (2.34), 


(2.35)  J|K*vn(t)|  dt  S o-1  -|v  | 


L>1  [0  ,T;H] 


< c |v  | 


L1  [ 0 , T ; H ] 

as  |u(t)|  = 1 almost  everywhere  in  t. 

By  passing  to  the  limit  in  (2.35),  we  deduce  that  K*  can  be  extended  to 
a bounded  operator  from  L^O.TjH]  into  L1[0,T;U]. 

// 

Coming  back  to  the  proof  of  (a2),  it  is  therefore  sufficient  to 
show  that 

^ e-^(L[0,T;Lp(fl)]  -*•  H ^ ^ (E ) ) for  1 < p < <*>. 


We  shall  now  derive  the  formula  for  L*v(t)  in  an  abstract  context.  This 
abstract  formulation  will  be  extensively  used  in  the  subsequent  chapters 
and  its  derivation  is  straightforward. 

Let  T be  given  by  the  formula  (2.10),  i.e. 
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t 

(2.36)  Tu(t)  = j S(t-t)Bu(t)dt. 

0 


£ 

Then  by  definition  of  the  adjoint  T , we  have 


(Tu.v) 


L [0,T;H] 
2 


T T 

J(Tu(t) ,v(t) ) dt  = J (u ( t ) ,T*v(t)..dt 
0 H 0 U 


Also, 


= (u,T*v) 


L2[0,T;U]' 


T T t 

J(Tu(t),v(t)  dt  = J(J  S(t-T)Bu(T)dT,v(t))„dt 
0 0 0 H 


T t 

= J /(S(t-T)Bu(T) ,v(t) ) dtdt. 
0 0 H 


By  reversing  the  order  of  integration,  we  have 
T T T 

J(Tu(t),v(t))Hdt  = J J(u(T),B«S*(t-i)v(t))udtdT 


T T 

= J(u(t),  J B*S*(t-T)v(t)dt)..dT 
0 T U 


T 

= J(u(t)  ,T*v(t)  )..dT. 
0 U 

Hence 

T 

(2.37)  T*v(t)  - B*  / S*(T-t)v(x)dT. 

t 
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By  taking  H = U = L2(r)  and  B = AD  formally  with 

3 

q=]j-+e,  e > 0 we  deduce  the  formula  for  L*v(t),  i.e, 


T 

L*v(t)  = D*A*  j S*(T-t)v(x)dT  on  L (ft) . 

t d 


Now  we  define  the  auxiliary  operator  M*  by 


T 

M*v(t)  = j S*(T-t)v(T)dT. 
t 

Our  strategy  is  to  first  establish  a certain  regularity  for  M*  and  then 
to  quote  a standard  result  from  trace  theory  to  conclude  the  argument. 

We  recall  that  v is  assumed  to  be  in  Lp[0  ,T;  L2(ft) ] for  1 < p < <*>. 

By  applying  Theorem  2.1  with  U = H = L2(Q)  and  B = I with  q = 0,  we  have 


|M*v|  £ c|v|  , 1 < P < ®. 

Lp[0,T;D(A*)]  Lp[0 ,T; L2(n) ] 


Using  the  identification  [7] 


D(A)  = D(A*)  = H2(ft)  n (ft) 

o 

and 


Ixl  2B  < cjxj  for  xeD(A*B),  0 < B < 1 , 

fTB(ft)  D(A*6) 


we  can  see  that 


|m*v|  £ c jv  | 

L>p[0,T ;H  (ft) ] 


Lp[0,T;L2(fl)] 


1 < p < ®. 
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As  for  the  time  regularity  of  M*v, 

T 

dM#v  r 

(t)  = _v(t)  + J A*S*(T-t)v(T)dT. 


By  applying  Theorem  2.1  again  with  U = H = ( £2 ) and  B = A with  q = 1, 

we  get 


«dM*v . 

* dt  ' 


c v 


Lp[0,T;L2(fl)] 


Lp[0,T;L2(n)] 


1 < p < ®. 


Thus  we  have  shown  that 


(2.38)  M*  e 1 (Lp[0,T;L2(0)]  -►Hp’1(Q)), 


1 < p < ®. 


To  finish  the  argument,  we  claim  that 


(2.39)  D*A*v  = for  v e D(D*A* ) , 

3vA 


where 


9 v 

3v, 


n 

l 


i»j  = 1 


a.  ,(x)n. (x) 


3v 

3x 


and  n^x) 


are  the  components  of  the 


outward  unit  vector  normal  to  the  boundary  r at  x.  This  can  easily  be 
seen  by  applying  the  second  Green's  formula  to  A(x,9): 


(A(x,3)u,v) 


(u,A*(x,3)v)a 


r 


By  choosing  u = Dg  and  veD(A*)  (i.e.  v|  = 0),  we  obtain 
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-(Dg,  A* (x  , 3 )v ) = -<g  ,D*A*v>  = - <g,|^-> 

*4  1 dv.  1 

A 


for  ail  v c D (D* A*  ) , g e D(D*A* ) ' . 

The  following  theorem,  which  was  originally  formulated  in  the 
setting  of  spaces  [13,  Vol . II,  p.  9],  can  be  extended  to  Lp  spaces 
for  1 < p < oo  and  provides  us  with  the  regularity  result  we  need. 


v-  ,v. 


Theorem  [13].  The  mapping  from  Hr,S(Q)  into  H J J(E)  which  to 

p p 


. . 9^u  . 

u associates  — r is  continuous,  where 


. !j  = r-J-  Vp 

r s r 


for  r > V2  , s £ 0 . 


By  taking  r=2,s=1,j=1,we  obtain 


// 


I— I 

9vA  Hp/2  ’ ^ (E) 


- CM  2 ! 

h;’  (Q) 


and  thus  finally 


P P 


(I) 


C |M*V  I 


Hp ’ 1 (Q) 


* c|v| 


Lp[0,T;L2(£j)] 
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This  concludes  the  proof  of  (83)  • 


Remark . Theorem  1.2  is  a direct  consequence  of  Theorem  2.3  since  y 
in  Theorem  1 .2  and  L in  Theorem  2.3  are  related  by 


y(t)=Lg(t),  t>0. 


// 


CHAPTER  III 
SPACE  DISCRETIZATION 

3.1  Opening  Remarks 

In  this  chapter  we  shall  introduce  an  approximation  scheme  which, 
given  a family  of  finite-dimensional  subspaces  Vh  and  a family  of 
finite-dimensional  operators  Ah,  produces  for  any  fixed  value  of  t an 
element  uh(t)  in  Vh  which  approximates  at  time  t the  solution  u(t). 

This  approximation  scheme,  which  is  generally  called  continuous- time 
approximation  because  the  discretization  does  not  affect  the  time 
variable,  will  be  presented  in  an  abstract  context.  An  optimal  abstract 
convergence  result  will  then  be  given.  Both  the  approximation  scheme 
and  the  convergence  result  are  found  in  Lasiecka  [12]. 

The  abstract  scheme  and  its  convergence  estimates  are  then  applied 
to  the  case  of  our  canonical  example  (2.1).  a family  of  standard 
approximating  subspaces,  as  well  as  the  standard  Galerkin  approximation 
^h  of  the  operator  A,  are  shown  to  give  rise  to  an  approximation  scheme 
with  optimal  convergence  properties  by  virtue  of  the  applicability  of 
the  abstract  scheme  to  these  choices. 

We  emphasize  here  that  this  applicability,  based  on  the  uniform 
analyticity  of  the  semigroup  generated  by  -Ah,  is  a nontrivial  result 
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which  has  been  established  only  recently  in  the  selfadjoint  case  [8,  1 ^ , 
19],  and  in  the  general  case  [.11].  This  point  will  be  further 
elaborated  in  sections  3.2  and  3.^. 

3.2  Galerkin  Approximation  of  an  Abstract  Boundary-Value  Problem 

In  this  section  we  shall  provide  a general  framework  for  the 
approximation  of  abstract  boundary-value  problems. 

Let  h be  the  parameter  of  discretization,  with  0 < h < 1 and  with 
its  values  approaching  zero.  Let  V^cD(B#)  be  a family  of  finite- 
dimensional subspaces  of  H and  let  Ph  be  the  orthogonal  projection  of  H 
onto  V^.  We  shall  assume  that  the  subspaces  Vh  enjoy  the  following 
approximation  properties: 

(3.1)  |x-P  x|  < ch2ot|x|  , 0 <,  a < 1 , 

H D(Aa) 

(3.2)  |B* (I-P  )x | < ch2(1_q)|x| 

D(A* ) 

(3.3)  |B*(I-P  )x|  < ch2a|x|  ,0  < a < 1-q, 

D(A*“  Q) 

(3.^)  " c!l  ^ lvhlH  (inverse  approximation  property), 

for  some  q < q. 

Another  inverse  approximation  property  which  is  stronger  than  (3. A)  is 
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0 < a £ 1 . 


We  will  need  to  assume  (3.4')  in  the  context  of  Theorem  5.2  rather  than 
(3.4). 


As  a preview  to  section  3.4,  we  point  out  that  the  inverse  assumption 
(3.4')  is  still  reasonably  general  as  it  is  satisfied  for  example  by 
spaces  of  splines  defined  on  uniform  meshes. 

The  operator  A will  be  approximated  by  a family  of  finite- 
dimensional operators  A^:  -*•  satisfying  the  following  generic 

properties: 


where  Sh(t)  = exp  (-Aht).  Property  (3.6)  is  referred  to  as  uniform 
analyticity  of  the  semigroup  Sj^ ( t ) . 

We  also  assume  that  the  properties  (3.5)  - (3.6)  hold  for  A*. 
Equivalently  to  (3.6),  we  have  [11] 


“ 1 — 1 2 

(3.5)  | (A  -Ah  P^)x|H  - ctl  ||XIH  (convergence), 


(3.6) 


, 0 £ g £ 2 uniformly  in  h, 


(3.7)  p(-A^)  I = {A:  ReA  > a - b | ImA  j } for  some  a,  b > 0 

and 


|RCi.-Ah)P.  |„^H  S -2. 


h h “ 


< 


(3.8) 


uniformly  in  h. 
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(3.7)  and  (3.8)  can  be  compared  with  (2.2)  and  (2.3)  in  section  2.2. 

Remark . If  A is  coercive  in  the  norm  of  D(A  ^2  ) (which  is  the  case 
for  the  canonical  example  2.1),  then  the  Galerkin  approximation  A^ 
defined  by 

(Ahxh'yh)H  = (Axh’yh)H  e Vxh’yh  e Vh£  D(A  /2  ) 


satisfies  (3.5)  and  (3.6)  (and  consequently  (3.7)  and  (3.8)).  Although 
this  result  is  to  be  expected,  it  has  been  shown  only  recently  [11]  and 
its  proof  is  technical. 

// 

Let  u(t)  be  the  solution  of  the  abstract  boundary-value  problem 

(2.6).  We  shall  consider  the  following  approximating  scheme: 
find  uh(t)  Vh  such  that 


(3.9) 


, (ah(t)-  vh>H  * •<Vh(t>’vh)H  4 <e<t>,B*Vu  Vvh£vh 
|.V°>  - 0 


Problem  (3.9)  can  be  considered  as  a semidiscrete  approximation  of 
problem  (2.6).  To  see  this,  we  define  the  operator 

Th:  L2[0,T;U]  -*■  L2[0,T;Vh]  corresponding  to  the  input-solution  operator 

T for  the  problem  (2.6),  by 

t 

(3.10)  T g(t)  = u (t)  = J S (t-i)P  (Bg)(T)dT, 

0 n 


-61- 


which  is  the  (formal)  semigroup  formula  for  the  solution  of  (3.9).  Thus 
(3-10)  is  an  approximation  of  (2.10). 

Remark . The  operator  T^  is  well  defined  for  any  input  g e Lp[0,T;U] , 
1 £ p £ °>.  This  can  be  seen  by  using  the  properties  of  Pp  and  the 
subspaces  V^.  In  fact,  for  geU,  and  from  the  orthogonality  of  Ph,  we 
have 


(Ph(Bg),v)H  = (Ph(Bg),Phv)H  = (Bg,Phv)H  VveH. 

Now  let  vh e Vh c D(B*).  Then 
(Bg.vh)H  = (g.B*vh)u  VgeU,  VvheVh. 

Thus  we  have  shown  that  Ph(Bg)  is  well-defined  in  H for  geU,  which  in 
turn  shows  that  Th  is  well-defined. 

// 


3-3  Estimation  of  the  Error 


In  this  section  we  shall  provide  results  on  the  rate  of  convergence 
of  the  solution  of  (3.9)  to  the  solution  of  (2.6). 

Let  e(t)  = u(t)  - u^(t)  = (T-Tj1)g(t),  where  u(t),  u^(t)  are  the 
solutions  to  (2.6)  and  (3.9)  respectively.  As  in  the  proof  of  Theorem 
2.3  (regularity ) , we  shall  use  Lemma  2.M  (duality  result)  by 
establishing  error  estimates  for  (T*-t£).  By  virtue  of  Lemma  2.4,  the 
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same  estimates  will  then  hold  for  (T-T^).  Thus  in  order  to  establish 
the  rates  of  convergence  for  e(t)  it  is  enough  to  estimate  T*-T^. 

We  shall  now  derive  the  formula  for  the  adjoint  T*  (with  respect  to 
the  topology  of  L2[0,T;U]  - L2[0,T;H]). 

Similarly  as  for  T* , we  have: 


(Thg,v) 


L2[0,T;H] 


T 

J (Thg(t)  ,v(t) )Hdt 


Because  Ph  is  an  orthogonal  projection,  the  following  equality  holds: 


(uh>v>H  = (uh-Phv>H  VveH,uheVh. 

Hence 

T T 

J(Thg(t),v(t))Hdt  = J(Thg(t),Phv(t))Hdt 

The  semigroup  formula  for  T^  is  then  derived  as  in  section  2.3.  Thus  we 
have 

T 

(3.1D  Tjphv(t)  = B*  J S*(T-t)Phv(T)dt 

t 

# 

*h  Gan  immediately  be  extended  to  Lp[0,T;H],  1 £ p < °°  by  setting 
T*v(t)  = T^Phv(t). 

We  now  state  a preliminary  convergence  result. 
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Theorem  3 • 1 • Let  V^cD(B  ),  0 < h < 1 , be  a family  of  finite- 
dimensional subspaces  of  H satisfying  properties  (3.1)  - (3.4)  and  let 
A^,  0 < h < 1,  be  a family  of  finite-dimensional  approximations  of  A 
satisfying  (3.5),  (3.6).  Let  T*  and  T*  be  defined  by  (2.37)  and  (3.11) 
respectively.  Then,  with  q as  defined  in  (3.1)  - (3.4), 


|(T*-T*)v|  < ch2(1"q)ln  1 |v|  , P-1,-. 

L [0,T;U]  n L [ 0 , T ; H ] 

P // 


Remark.  This  result  is  derived  in  Lasiecka  [12].  For  the  sake  of 


completness  we  provide  its  proof  in  chapter  VI. 


// 

By  Theorem  3.1  and  Lemma  2.4  we  immediately  deduce  the  following 

theorem. 


Theorem  3.2.  Let  V^.A^,  0 < h < 1 , satisfy  the  same  assumptions  as  in 
Theorem  3.1.  Then 


|e| 


Lp[0,T;H] 


= |u-uhl 


Lp[0 ,T;H] 


ch2(1_q)ln  1 |g| 


Lp[0,T;U] 


P = 1 


// 


3.4  Application  to  a Canonical  Parabolic  Boundary- value  Problem 


We  recall  the  parabolic  problem  (2.1): 


on  Q 
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Yt  - A(x,3)y 

y (o ) - o 


y|r  = b 


on  n 
on  I 


where  g e Lp[0,T; ( r) ] , 1 £ p £ 00  and  A(x,3)  is  a uniformly  strongly 
elliptic,  second-order  differential  operator  as  defined  in  (1.7). 

We  shall  now  see  how  the  abstract  approximation  scheme  of  the 
previous  section  can  be  applied  to  yield  the  optimal  rate  of  convergence 
of  the  Galerkin  approximation  to  problem  (2.1). 


As  we  have  seen  in  section  2.3,  we  can  represent  (2.1)  in  integral  form, 

i .e . , 

t 

y(t)  = j AS(t-T)Dg(x)dT 
0 

with  the  operators  A,D,S(t)  as  defined  in  that  same  section  2.3. 

We  have  also  seen  in  section  2.3  that  problem  (2.1)  can  be  formulated 
within  the  context  of  problem  (2.6)  by  taking  H = 1^(8),  U = L^r),  B = 
AD  (formally)  with  q = | + e,  e > 0 and  q = |. 

We  shall  formulate  an  approximating  scheme  for  problem  (2.1)  by 
introducing  a standard  class  of  finite-dimensional  subspaces  Vh  [2]  and 
a family  of  finite-dimensional  operators  Ah  which  comply  with  the 
abstract  conditions  (3.1)  - (3.4),  (3.4')  and  (3.5),  (3.6). 

We  recall  the  definition  of  an  S^,k-system  as  given  by  Babuska  and 
Aziz  [2], 
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p k 

Def inltion.  For  0 < h < 1 , let  S^’  (8)  be  a linear  system  of 
functions  defined  on  8,  with  r > k Z 0,  satisfying  the  following 
conditions . 


(i)  S[j,k(8)  eHk(8) 

0 

(ii)  for  every  veH  (8),  l i 0 and  0 S s £ min  (S,,k),  there  exists 
vheSh,k(8)  such  that 


(3.12) 


|v-vj 


HS(8) 


chu |v | 


H4(8) 


where  y = min  (£-s,r-s)  and  the  constant  c is  independent  of  v and  h. 
The  function  is  dependent  on  v and  may  be  different  for  different 

values  of  s. 

// 

r k 

Definition.  An  S^’  -system  is  said  to  be  regular  if  the 
approximating  function  vh  can  be  chosen  independently  of  s so  that 
(3-12)  holds  for  all  s,  with  0 < s £ min(i,,k). 

// 

We  now  recall  the  Inverse  Approximation  Property.  We  say  that  the 
p k 

system  S^’  (8)  satisfies  the  inverse  approximation  property  if  there 
exists  e,  0 < e ^ k such  that  for  every  k - e £ s £ k and  every 
vh  e S[j’k(8)  we  have 


lvhl  k 

HK(8) 


0h'U"3)|vhl  s 

HS(8) 
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where  c does  not  depend  on  h or  v^. 
Then  we  have  the  following  result. 


1°  k 

Theorem  [2].  Let  S^’  (ft)  be  an  (r,k)-system  which  satisfies  the  inverse 

assumption.  Then  for  every  s-,  ,s?  with  s < s < k and  vues[,k(ft),  we 

' c 1 2 h h 

have 


(ft) 


-(s 


< ch 


2 1 


(ft) 


where  c does  not  depend  on  h or  v^  (it  depends  in  general  on  s1  and 

s2^*  // 
We  are  now  in  a position  to  define  our  family  of  finite-dimensional 
subspaces  Vh. 

2 1 

Let  Vh  be  a regular  S'  (ft)  system  satisfying  the  inverse 
assumption.  Let  Ph  be  an  orthogonal  projection  from  L^(ft)  onto  V^. 

Then 


(3.13)  |v-p  v|  < cha  S | v j , 0 < a < 2,  0 ^ s < 1,  a-s  > 0 

1 1 s ot 

We  also  require  that  the  subspaces  Vh  satisfy 
(3. 11*)  |V(Phv-v)|  < ch  |v|2, 

(3.15)  |V(Pnv-v)|  S oh“|v|(2o^c)/2,  0<a<V2. 

if^l  s oh"3/2|v|.  ¥vevh- 

A 


(3.16) 
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Remark  . It  is  well-known  that  spaces  of  splines  defined  on  uniform 
meshes  satisfy  the  properties  of  a regular  S ’ '-system  with  the  inverse 
approximation  property.  For  example,  linear  splines  defined  on  a 
bounded  domain  n with  uniform  grid  form  a regular  S^’  (ft)-system  with 
the  inverse  approximation  property.  Furthermore,  they  also  satisfy 
conditions  (3.14)  - (3.16). 

// 

As  an  approximation  of  A we  take  an  arbitrary  operator  Ah:  Vh  -*• 


such  that 

(3.17) 

l(A_1'Ah' 

ph)x|  S 

ch2 |x | , 

(3.18) 

lAhThl  £ 

cl*BvhI 

, 0 < 8 <V2 

and 

(3.19) 

l(AhVvh 

)|  S c| 

Uhll  lVhll  * 

(3.20) 

(WV 

2 

|2,  Vuh,vh£Vh. 

At  this  point  we  need  to  make  a few  remarks. 

(i)  Most  of  the  well-known  approximations  Ah  to  elliptic  problems 
comply  with  the  conditions  (3-17)  - (3*20).  For  example,  the  standard 
Galerkin  approximation  Ah,  where 

(Wvh>  ■ <AvV 


(3.21) 


Vuh’vhe  Vh. 
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satisfies  (3.17)  - (3.20).  Other  methods,  such  as  Babuska's  penalty 
method  Cl],  and  Nitsche's  method  [17],  also  satisfy  (3.17)  - (3.20),  and 
moreover  do  not  require  the  subspaces  to  satisfy  zero  boundary 
conditions  (as  the  Galerkin  method  does). 


(ii)  Using  the  identifications  [7] 

(3.22)  D(Aa)  = D(A*a)  = H^a(n),  0 < a < |, 

(3.23)  D(A)  = D(A* ) = H2  (fi ) n Hq (£J ) , 
and 

(3.214)  |x|  Sc  |x  | for  X£D(A6),  0 S B S 1, 

* D(A^) 

we  can  see  that  the  conditions  (3.13)  - (3.16)  imply  the  corresponding 
abstract  conditions  (3.1)  - (3.4). 


In  fact,  by  taking  s = 0 in  (3.13),  we  have  by  (3.214) 


| v-P  v | S ch2a|v|  < ch2a|v|  , 0 < a £ 1. 
h 2“  ' D(Aa) 


Also  with  B*  = D*A*  (and  thus  q = 3/4+e, 

av 


section  2.4  that  D*A*v 


9v, 


-,  we  have  by 


e>0),  and  recalling  from 
(3.23)  and  (3-24): 


-69- 


and 


IVCP.^v-v)  *v^  |=  |D*A*(P  v-v)  | £ ch2a|v|| 


< ch2a|v 

2a+-^+2e 


0 $ a <i  Vi)  -e , e>0 


and  finally 


(iii)  Condition  (3.V)  is  automatically  satisfied  since  Vh 
satisfies  the  inverse  assumption. 

(iv)  it  was  shown  in  Lasiecka  [11]  that  with  Ah  satisfying  (3.17)  - 

(3.20),  the  corresponding  semigroup  Sh(t)  is  uniformly  analytic,  i.e., 


which  is  precisely  the  condition  (3.6)  with  H = L2 (£2 ) . This  result 
parallels  the  well-known  result  mentioned  in  Chapter  II,  which  states 
that  if  A is  a uniformly  elliptic  operator,  then  -A  is  the  generator  of 
an  analytic  semigroup. 

(v)  Finally,  (3.5)  follows  directly  from  (3.17). 


Thus  we  have  proved  that  the  assumptions  (3.1)  - (3.6)  follow  from 
(3.13)  - (3.20). 


lAhSh(t)xhl  " lHxhl  uniformly  in  h, 

t 


// 
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From  these  remarks  we  can  see  that  the  abstract  approximation 
scheme  can  be  applied  to  our  situation.  By  taking  H = ( ^ ) » 

V = L2(D,  B = AB  (formally)  with  q = | + e,  e > 0,  and  recalling  from 
(2.39)  that 


B*v  . D*A*v  - 

3v» 

the  algorithm  takes  the  form: 
find  y^(t)  e vh  such  that 


3v, 


(3.25) 


l2(o) 


l2(u) 


A L2(D 


yh(o)  = o 


for  all  v.  e V.  . 

h h 


If  Ah  is  a standard  Galerkin  approximation  (see  3.21),  then  (3.25) 
becomes 


3vh 

<vtKVL2<flr  -a(y*(t)’V  * 4(r) 

where  a(u,v)  is  the  bilinear  form  associated  to  A. 

The  following  convergence  result  associated  with  the  algorithm 
(3.25)  follows  directly  from  Theorem  3.2. 
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Theorem  3-3.  Let  e(t)  - y(t)  - y^(t),  where  y(t)  satisfies  the 
parabolic  equation  (2.1),  and  y^(t)  is  its  continuous-time  approximation 
as  given  in  (3.25).  Then 

lel  < ch  /2ln  J |g|  , P = 1 , «. 

L [0,T;L2(n)]  n L [0,T;L  (T)] 

P // 

Remark . In  view  of  the  best  approximation  properties  of  the 
subspaces  and  of  the  maximal  regularity  of  the  solution  y(t) 

1/p -2e  1/ji-e 

(Hp  ’ (Q)  for  P = 1 , 00  and  q = 3/^  + e: ) , it  follows  that  the 

convergence  estimate  0(h1/2in^)  is  optimal. 

The  same  result  for  1 < p < ® can  be  obtained  by  interpolation. 

Such  a result  is,  however,  not  optimal  (modulo  ln^) , since  the  maximal 
regularity  of  the  solution  is  H ’ 1/j*  (Q)  for  1 < p < «.  The  optimal 
convergence  rate  0(h  ^2  ) has  been  established  and  requires  more  refined 
techniques  found  in  Lasiecka  [12]. 

// 


CHAPTER  IV 
TIME  DISCRETIZATION 


4.1  Opening  Remarks 


In  this  chapter  we  shall  introduce  an  approximation  scheme  which 
complements  the  scheme  given  in  chapter  III.  Whereas  the  previous 
scheme  approximated  the  values  of  the  solution  for  a given  value  of  t 
and  left  the  time  variable  t unchanged,  in  this  scheme  we  shall  do  the 
opposite,  namely  we  shall  approximate  the  solution  at  discrete  values  of 
the  time  variable,  and  otherwise  leave  the  range  of  values  of  the 
solution  unchanged. 

In  section  4.2  we  shall  present  the  basic  one-dimensional  backward 
Euler  difference  scheme,  which  is  used  as  a basis  for  an  abstract 
discrete-time  approximation  scheme. 

The  next  two  sections  are  devoted  to  the  development,  justification 
and  statement  of  the  convergence  properties  of  this  approximation  scheme 
in  the  context  of  an  abstract  boundary-value  problem. 

Finally,  in  the  last  section  4.5,  the  abstract  convergence  results 
are  applied  to  the  canonical  example  (2.1)  and  convergence  estimates  are 
derived  which  are  seen  to  be  essentially  optimal. 
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4.2  A One-Dimensional  Euler  Scheme 


The  time  discretization  scheme  we  shall  use  in  this  work  is  based 
on  the  backward  Euler  difference,  which  we  describe  here  in  a simplified 
context.  The  general  scheme  will  be  given  in  the  next  section  as  an 
extension  of  this  basic  scheme.  Consider  the  one-dimensional  initial- 
value  problem 


where  xej,  J an  interval  on  E,  and  a(x)  > 0 Vx  e J. 

Let  k be  the  parameter  of  discretization,  with  k > 0.  We  let 

fn  = f(nk),  and  let  yn  be  the  time-discrete  approximation  of  y(nk).  By 


3y  1 

replacing  T^-(nk)  with  - (yn~y  ) , the  scheme  takes  the  following 
form: 


find  yn,  n = 1 ,2, 


• • • 9 


where 


, n = 0, 1 


9 • • • 


yQ  given. 


By  solving  for  yn+1 , (A  .2)  becomes 


n = 0,1 


9 • • • 


If  we  set  ek  = 1/(1+ka),  then  the  recursive  formula  (4.3)  can  be  put 
into  the  explicit  form 


(4.4) 


yn 


n 

e,  y + 
k Jo 


I 


j = 1 


n = 1,2,... 


to  be  compared  with  the  formula  for  the  actual  solution  of  (4.1) 
obtained  by  the  variation  of  parameter  formula: 

t 

(4.5)  y (t ) = e'aty0  + J e~a(t_x)f ( t )dr , t > 0 

0 


For  t = nk,  we  have 

nk 

(4.6)  y(nk)  = e ankyQ  + Je  a^nk  T^f(x)dx 

n . k 

= (e  ) y + (e  ) J e f(jk+T-k)dx 

j = 1 0 

The  obvious  questions  about  stability  and  convergence  of  the  scheme 

will  be  answered  in  the  more  general  context  of  sections  4.2  and  4.3. 

One  can  nevertheless  notice  that  ek  = 1/(1 +ka)  is  a second-order 
approximation  in  k of  e"ak  for  "small"  k.  This  particular  property  is 
the  heart  of  the  proof  of  a result  which  itself  is  crucial  in 
establishing  stability  and  convergence  (see  Lemma  4.1,  and  its  proof  in 
chapter  VI ) . 
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4.3  Extension  of  the  Scheme  to  an  Abstract 
Parabolic  Boundary-Value  Problem 


From  section  2.2  we  recall  the  abstract  problem 


(4.7) 


u = Au  + Bg 


u ( 0 ) = u 


on  D(A* ) ' 


where  -A  is  the  infinitesimal  generator  of  a strongly  continuous 
analytic  semigroup  S(t),  B is  an  operator  from  a Hilbert  space  U into 
D(A*) ' with  A~qB  e I (U+H)  for  some  0 < q < 1 , and  g e Lp[0,T;U] . 

In  order  to  apply  the  scheme  described  in  the  previous  section  to 
the  problem  (4.7),  we  define  the  operator 


Pk:  Lp[0,T;H]  + Lp[0,T;H] 

by  .. 

Jk 

(4.8)  Pkv(t)  = l j v ( t )dt 

( j-1  )k 

where  T = Nk  for  some  integer  N and  j is  the  unique  integer  such  that 
(j-1  )k  < t < jk. 

The  range  of  Pk  consists  of  the  step  functions  defined  on  [0,T] 
with  (discrete)  values  in  H.  At  this  point  the  following  remarks  are  in 
order . 

(i)  Pk  e £ (Lp[0  ,T;H]  -*L  [0,T;H]),  1 £ p ^ Furthermore,  Pkv(t) 
is  defined  for  every  te[0,T]. 
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(ii)  Pk  is  an  orthogonal  projection,  whch  implies  that 
selfadjoint  operator  (with  respect  to  the  inner-product  in  L 

(iii)  Pk  can  clearly  be  defined  the  same  way  on  Lp[0,T; 
same  properties. 


We  now  define  the  discrete- time  approximation  uk  of  (4. 
solution  of  the  difference  equation 


(4.9) 


-(u  ( (n+1 )k)-u  (nk) ) = -Au  ((n+1)k)  + BP.  g( (n+1 )k) 


u,,(0)  given 


which,  after  solving  for  uk((n+1)k),  becomes 


-1 


(4.10) 


u ((n+1)k)  = (1+kA)  [u.  (nk)  + BP  g ( (n+1 ) k ) ] 


uk(0)  given. 


Following  the  same  steps  as  in  section  4.2,  we  obtain 


(4.11)  uk(nk)  = e"  uq  + ^ Ek+1_jBPkg(jk),  n = 1 ,2, 


where  Ek  = (1+kA)-1  = 1 R(i,  -A). 


We  can  define  uk  for  other  values  of  t in  the  obvious  way,  i 


Pk  is  a 
2[0,T;H]). 

U]  with  the 

// 

7 ) as  the 


uk(t)  = uk(nk),  (n-1 )k  < t < nk,  n = 1,2,... 
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Again  for  simplicity  and  without  loss  of  generality,  we  assume  that 
uQ  = 0.  We  now  make  the  following  remarks  to  justify  the  validity  of 
formula  (4.11). 

(i)  By  (i)  and  (iii)  of  the  last  set  of  remarks,  we  deduce  that 
Pkg(jk)  is  defined  for  j = 1,2,... 

(ii)  The  exponent  (n+1-j)  of  Ek  in  formula  (4.11)  is  always  at 
least  1;  moreover, 

AEk  e £ (H  •*  H)  for  k>0 

which  can  be  readily  seen  from  the  resolvent  identity 

ARU.-A)  = I - AR( A ,-A) . 


Therefore,  EkB  e £ (U  -*•  H)  since 


(EkBu,v)H  = (AEkA  ^Bu.v^  Vu  e U,  veH. 


// 


Thus  the  formula  (4.11)  makes  sense  in  the  topology  of  H. 


4.4  Preliminary  Results 

In  order  to  study  the  stability  and  convergence  of  the  abstract 
scheme,  we  shall  need  to  establish  the  approximation  properties  of  Pk. 
From  the  orthogonality  of  Pk,  we  have 
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(4.12)  | (I-P  )x  | S cjxj  , p = 1 

L [0,T;H]  L [0,T;H] 

Also,  p p 

(4.13)  | (I-P  )x|  £ ckjx|  i P = I,00. 

Lp[0,T;H]  Hp[0 , T ; H ] 


In  fact , for  p = 1 , », 


|(I-P  )x|  = ess  sup  (jx(t)-P  x(t)  |j 

LjO.TjH]  0<t<T  H 

Jk 

= max  {ess  sup  |-^  J (x( t )-x( i) )dt | } . 

j = 1 N (j-1  )k<t<jk  k (j-1  )k  ri 


By  the  Mean  Value  Theorem, 


jk 

|(I-P  )x|  £ max  {1  J k sup|*U)|  di} 

L.[0,T;H]  J-1 N (j-1  ) k 5 


< k |x  | 

LJ0,T;H] 


< k |x| 

H [0,T;H] 

CO 


For  p = 1 , we  have  similarly 


T 

|(I-P  )x|  = / jx(t)-P  x(t)  | dt 

L1 [0  ,T; H]  0 K H 


N Jk 

= I / |x(t)-P  x(t)|  dt 

j-1  (j-1 ) k k H 
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= I S J (x c t )-x( t) )dT  I dt 

j-1  (j-1  )k  (j-1)k 


N J 

I I 

J-1  (j-Dk 


jk  t 

I |J  x(  £,)d^  | dTdt 
(j-1 )k  t H 


jk 


jk 


t:  1 S I |/|x(OLde|dTdt 

j-1  (j-1  )k  (j-1 )k  t 


N 


jk  jk 


jk 


M J 


k L J J S l*U)L  d^dtdt 

j-1  (j-1  )k  (j-1 ) k (j-1  )k  H 


S k |x| 


L1 [0 ,T;H] 


* k |x| 


Hj[0,T;H] 


A well-known  interpolation  theorem  [13,  Vol . I,  p.  27]  then  enables 
us  to  conclude  that 


(4.14) 


| (I-p  )x | < cka  jx | , p-1,  0 S a S 1 

Lp[0  ,T; H]  Hp[0 , T ; H ] 


for  x e Hp[0 , T; H] . 


The  estimate  (4.14)  can  then  be  extended  by  density  to  xeHa[0,T;H]. 
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We  shall  make  strong  use  of  (4.14)  in  the  proof  of  Theorem  4.2  and 
also  in  that  of  Theorem  5.2. 

We  now  state  another  result  which  we  win  need  in  order  to 
establisn  stability  and  convergence  of  the  abstract  approximation 
scheme . 

Lemma  4.1.  Let  A:  H+H  be  an  arbitrary  closed  linear  operator  such  that 
A is  the  infinitesimal  generator  of  a strongly  continuous  analytic 
semigroup  S(t).  Without  loss  of  generality,  we  assume  that 

o(A)  £ i = {X:  4><  | arg A | ^ it  for  some  <f>>-^}\{0} 

Let  k > 0 be  fixed  and  let  R(^-,A),  where  R(A,A)  is  the  resolvent 

K K 

of  A.  Then 


(4.15) 

and 

K-S(nk))x|H  S S |X|H, 

n = 1,2,... 

(4.16) 

|A(E"-S(nk))x|H  S |x|H, 

n = 1,2,... 

The  constant  c does  not  depend  on  n or  k.  It  depends,  however, 
on  4>,  as  can  be  seen  in  the  proof.  The  proof  of  Lemma  4.1  is  given  in 
Chapter  VI. 

// 

In  order  to  state  our  theorem  on  stability  and  convergence  of  the 
discrete-time  approximation,  we  need  to  define  an  input-solution 
operator  that  we  can  compare  with  T,  i.e., 
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Tk:  Lp[0 ,T ;U]  ■>  Lp[0,T;H], 

where 

n + 1 - ' 

(^•17)  Tkg(t)  = uk(nk)  = k l E^+  JB(P  g)(jk)  for  (n-1)k<t£nk. 

j = 1 K 


Since 
operators , 


the  theorem  will  be  stated  and  proved  in  terms  of  adjoint 

£ 

we  derive  the  formula  for  Tk.  First,  we  note  that 


(Tkg,v)  _ n = (Tkg,Pkv) 

L2[0,T;H]  K K L2[0,T;H] 


by  orthogonality  of  Pk.  Then, 


(Tkg>pkv) 


L2[0,T;H] 


= J (Tkg(t) ,PRv(t) )Hdt 


nk 


= l I (Tkg(t),Pkv(t))  dt 
n=1  (n-1)k 


By  definition  of  Tk  and  Pk,  we  have 


Tkg(t)  = Tkg(nk)  and  Pkv(t)  = Pkv(nk)  for  (n-1 ) k < t £ nk. 


Therefore , 


(Tkg,pkv) 

Li 


[0 ,T;H] 


= l j (T  g(nk) ,P  v(nk) ) dt 
n-1  (n-1 )k  K H 


N 

= k l (Tkg(nk) ,Pkv(nk))H 
n=1 
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N n 


k 1 0<  l E^+1  jB(P  g)(jk),P  v(nk)) 
n=1  j = 1 k H 


N n 

r-*  r-'  n+  1 - 1 

k I k I (Ek  "B(P  g)(jk),P  v(nk)) 
n=1  j=1  k k H 


N n 


k I l ( P^g ( j k ) > B*E*n  1_JP  v(nk)) 
n=1  j = 1 K k j 


By  reversing  the  order  of  summation, 


(Tkg-pkv> 


L2[0,T;H] 


2 N N 

k l l (Pug(jk),B*E*n+1"JP.  v(nk)) 

j = 1 n=j  k k 


N N 

= k l (Pkg(jk),kB*  l E*n  JPkv(nk)) 
j=1  n=j 


g(t)dt,T  P v( jk) ) 
1 )k  K K 


= I / (g(t),T*P  v(t)Ldt 

j=1  (j-1 )k  K K u 


" (g’TkPkv) 


L2[0,T;U] 
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Therefore 


( ^ - 18 ) Tjpkv(t)  = kB*  1 E*n+1  jPkv(nk)  for  (j-1 )k  < t < jk. 

n=j 


* 

Tk  extends  to  the  whole  of  Lp[0,T;H]  in  the  obvious  manner. 
We  have  the  following  preliminary  theorem. 


Theorem  4.2.  With  T*  and  T*  as  defined  in  (2. 37)  and  (4.18) 
respectively,  we  have 


(4.19)  |(T*-T*)vj 


Lp[0 ,T; D ( A*~a) ] 


Lp[0,T;U] 


for  1 £p<°°,  0 £ a S 1-q,  e > 0,  where  A_<^B  e Jl  (U+H)  for  some  q, 

0 S q < 1 . 


The  proof  of  Theorem  4.2  will  be  given  in  chapter  VI. 

We  are  now  in  a position  to  state  our  stability  and  convergence 
theorem. 

Theorem  4.3.  Let  e(t)  = u(t)  - uk(t),  where  u and  uk  are  solutions  to 
the  abstract  problems  (2.6)  and  (4.9)  respectively,  and  given  by  the 
respective  formulas  (2.8)  and  (4.11)  with  uQ  = 0.  Then 


// 


Lp[0,T;D(Aa)] 


Lp[0,T;U] 


for  1 Sp£®,  0£a£  1-q,  e>0,  where  A e £ (U+H)  for  some 
q,  0 < q £ 1 . 

// 

Proof  of  Theorem  4.3.  We  note  first  of  all  that  u(t)  = Tg(t)  and 
uk(t)  = Tkg(t),  thus  e(t)  = (T-Tk)g(t).  Theorem  4.3  then  follows 
immediately  by  application  of  Lemma  2.4. 

// 


4.5  Application  to  a Canonical  Example 


Similarly  as  in  the  previous  chapters,  we  apply  our  abstract  result 
to  the  case  of  the  canonical  problem  (2.1),  i.e. 


!yt  = -A(x,3)y  on  Q 

y(0,x)  =0  on  fi 

y|f  = g on  E 


As  we  have  seen  earlier  in  chapter  II,  (2.1)  can  be  expressed  as  an 
abstract  boundary-value  problem  with  H = (f2 ) , U = L2(D,  B = AD 

(formally)  with  A-C,B  e i?(U->H)  and  q = 3/4  + e,  e > 0. 

Then  as  a discrete-time  approximation  of  (2.22),  i.e., 

t 

Lg(t)  = A J S(t-T)Dg(i)dT, 

0 


we  define  (see  formula  (4.17)) 
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(^•21)  Lkg(t)  = k l AEk+1  jPk(Dg)(jk),  (n-1 )k  < t < nk. 

j = 1 

By  taking  a = 0 in  Theorem  4.3,  we  have  the  following  convergence 
result  for  problem  (2.1). 


Theorem  A , A . With  L and  Lk  given  by  (2.22)  and  (4.21)  respectively,  we 
have 


JCL-L.  )g| 


£ ck 


% 


Lp[0,T;L2(a)] 


Is  I 


L [0,T;L2(D] 


1 < p 


e > 0 
// 


Remark . By  comparing  the  rate  of  convergence  in  Theorem  4.4 

1 A.  — ^ 

( 0 ( k 4 ))  with  the  time  regularity  of  the  solution  to  the  problem 

(2.1)  (which  is  Hp1/2*"e[0,T;L2(n)]  for  p - 1 , and  Hp1/4  [0,T;L2(Q)]  for 
1 < P < “),  one  can  see  that  such  a rate  of  convergence  is  essentially 
optimal  for  p = 1,  ®,  and  almost  optimal  (modulo  e)  for  1 < p < ®.  We 
shall  see  in  the  next  chapter  that  by  combining  the  two  discretization 

schemes  (continuous-time  and  discrete- time ) , the  rates  of  convergence 
obtained  for  the  fully  discrete  approximation  are  0(h1/2in-)  and 

i h 

% i 

0(k  In-)  for  1 < p < ®,  which  are  optimal  for  p = 1,  ® and  almost 

optimal  (modulo  In  ^ and  In  -)  for  1 < p < ®. 

h k K 


Finally,  the  stability  of  Lk  is  clear  by  taking  a = 1-q  =Vij-e.  Then 
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lLkgl  U-E  * lLg!  1/.  _F 

Lp[0  ,T ; D ( A /4  e)]  Lp[0,T;D(A  /4  E)] 


|(L-Lk)g|  1, 


Lp[0,T;D(A 


* c|g| 


Lp[0,T;D(A  1/4  E ) ] 


_C)] 

for  1 S p S 
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TABLE  1.  Values  of  log1  q (1  O^e^)  for  solution  Ui  in  the  norm  of  L-^  [0,T;L2  fi]. 


<—  log2  £ — > 


iog2  1 - 
1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 .58 

1 .59 

1 .59 

1 .59 

1 .59 

1 .59 

1 .59 

1 .59 

1 .59 

1 .59 

2 

1 .45 

1 .45 

1 .45 

1 .45 

1 .45 

1 .46 

1 .46 

1 .46 

1 .46 

1 .46 

3 

1 .32 

1 .32 

1 .30 

1 .30 

1 .30 

1 .30 

1 .30 

1 .30 

1 .30 

1 .30 

4 

1.18 

1 .18 

1.18 

1.15 

1.15 

1.15 

1.15 

1.15 

1.15 

1.15 

5 

1 .04 

1 .04 

1 .04 

1 .00 

1 .00 

1 .00 

1 .00 

1 .00 

1 .00 

1 .00 

6 

0.95 

0.91 

0.88 

0.87 

0.86 

0.86 

0.86 

0.86 

0.86 

0.86 

7 

0.86 

0.80 

0.75 

0.72 

0.72 

0.71 

0.71 

0.71 

0.71 

0.71 

8 

0.79 

0.70 

0.63 

0.59 

0.56 

0.56 

0.56 

0.56 

0.56 

0.56 

9 

0.73 

0.62 

0.52 

0.46 

0.43 

0.41 

0.41 

0.41 

0.41 

0.41 

10 

0.71 

0.57 

0.43 

0.34 

0.28 

0.28 

0.26 

0.26 

0.26 

0.26 
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TABLE  2.  Values  of  logi0(102ehk)  for  solution  u^  in  the  norm  of  L^O.TjL^O)]  . 

los2  i ■ 

1 

2 

3 

4 

5 

6 

7 

8 
9 


<---  Iog2  1 — > 

1234567  89  10 


1 .69 

1 .59 

1 .59 

1 .59 

1 .59 

1 .59 

1 .59 

1 .59 

1 .59 

1 .59 

1 .45 

1 .45 

1 .46 

1 .46 

1 .46 

1 .46 

1 .46 

1 .46 

1 .46 

1 .46 

1 .30 

1 .30 

1 .30 

1 .30 

1 .30 

1 .30 

1 .30 

1 .30 

1 .30 

1 .30 

1.18 

1.15 

1.15 

1.15 

1 .15 

1 .15 

1.15 

1.15 

1 .15 

1.15 

1 .04 

1 .00 

1 .00 

1 .00 

1 .00 

1 .00 

1 .00 

1 .00 

1 .00 

1 .00 

0.90 

0.87 

0.86 

0.86 

0.86 

0.86 

0.86 

0.86 

0.86 

0.86 

0.78 

0.72 

0.72 

0.71 

0.71 

0.71 

0.71 

0.71 

0.71 

0.71 

0.68 

0.59 

0.57 

0.56 

0.56 

0.56 

0.56 

0.56 

0.56 

0.56 

0.60 

0.48 

0.43 

0.41 

0.41 

0.41 

0.41 

0.41 

0.41 

0.41 

0.56 

0.38 

0.30 

0.26 

0.26 

0.26 

0.26 

0.26 

0.26 

0.26 

10 
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TABLE  3.  Values  of  log-,  q ( 1 02ehk ) for  solution  u2  in  the  norm  of  L-,  [0 , T ; L2  n]- 


<—  log2  l — > 


, 1 
1q82  H = - 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

1 .62 

2.53 

r 

1 .77 

1 .68 

1 .65 

1 .64 

1 .63 

1 .63 

r 

1 .63 

1 .63 

2 

2.08 

2.49 

1 .64 

1 .53 

1 .49 

1 .49 

1.48 

1 .48 

1.48 

1.48 

3 

2.23 

2.53 

1 .57 

1 .40 

1 .34 

1 .32 

1 .32 

1 -32 

1 .32 

1 .32 

4 

2.26 

2.56 

1 .53 

1 .28 

1 .20 

1.18 

1.18 

1.18 

1.18 

1.15 

5 

2.26 

2.56 

1 .51 

1 .20 

1 .08 

1 .04 

1 .00 

1 .00 

1 .00 

1 .00 

6 

2.26 

2.56 

1 .49 

1.15 

0.97 

0.90 

0.87 

0.86 

0.86 

0.86 

7 

2.26 

2.56 

1 .48 

1 .1 1 

0.89 

0.78 

0.73 

0.72 

0.71 

0.71 

8 

2.26 

2.56 

1 .48 

1 .08 

0.83 

0.67 

0.60 

0.57 

0.56 

0.56 

9 

2.26 

2.56 

1 .48 

1 .08 

0.78 

0.58 

0.48 

0.43 

0.41 

0.41 

10 

2.26 

2.56 

1 .46 

1 .08 

0.75 

0.52 

0.36 

0.30 

0.28 

0.26 
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TABLE  4.  Values  of  log-|  g ( 1 02ehk)  for  solution  u2  in  the  norm  of  L^CO.T^  fl]. 

<—  log 2 1 — > 


i°62  t - 
1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 .60 

2.48 

1 .71 

1 .63 

1 .61 

1 .61 

1 .61 

1 .61 

1 .61 

1 .61 

2 

2.04 

2.41 

1 .57 

1.48 

1 .46 

1 .46 

1 .46 

1 .46 

1 .46 

1 .46 

3 

2.20 

2.46 

1 .46 

1 .32 

1 .32 

1 .32 

1 .30 

1 .30 

1 .30 

1 .30 

4 

2.23 

2.48 

1 .41 

1 .20 

1.18 

1.15 

1.15 

1.15 

1.15 

1.15 

5 

2.23 

2.48 

1 .38 

1 .08 

1 .00 

1 .00 

1 .00 

1 .00 

1 .00 

1 .00 

6 

2.23 

2.48 

1 .34 

0.97 

0.87 

0.86 

0.86 

0.86 

0.86 

0.86 

7 

2.23 

2.48 

1 .34 

0.89 

0.73 

0.71 

0.71 

0.71 

0.71 

0.71 

8 

2.23 

2.48 

1 .34 

0.83 

0.60 

0.56 

0.56 

0.56 

0.56 

0.56 

9 

2.23 

2.48 

1 .32 

0.80 

0.49 

0.41 

0.41 

0.41 

0.41 

0.41 

10 

2.23 

2.48 

1 .32 

0.79 

0.43 

0.28 

0.26 

0.26 

0.26 

0.26 
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TABLE  5. 

(a)  Values  of  the  solution  at  T=  V2  calculated  for  the  boundary  input  gQ ( t ) = 1, 

g1 (t)  = 1 , 0 < t < V2  . 


X - 

i/8 

2/8 

3/8 

4/8 

5/8 

b/8 

7/8 

N = 

8 

.98121 

.96529 

.95467 

.95094 

.95467 

.96529 

.98121 

N = 

16 

.99350 

.98798 

.98430 

.98301 

.98430 

.98798 

.99350 

N = 

32 

.99582 

.99227 

.98991 

.98908 

.98991 

.99227 

.99582 

N = 

611 

.99633 

.99322 

.991 14 

.99041 

.99114 

.99322 

.99633 

N = 

128 

.99646 

.99345 

.99144 

.99074 

.99144 

.993*15 

.99646 

N = 

256 

.99649 

.99351 

.99152 

.99082 

.99152 

.99351 

.99649 

N = 

512 

.99649 

.99352 

.99153 

.99084 

.99153 

.99352 

.99649 

N = 

1024 

.99650 

.99352 

.99154 

.99084 

.99154 

.99352 

.99650 

(b)  Values  of  the  solution  at  T=  V2  calculated  for  the  boundary  input  gQ(t)=  1 , 
0 < t £ V2  and 

gl  (t)  = 0,  0 < t £ 1/8,  V24  < t S 3/8 

g1(t)-1,1/8<tS1^l  3/8  < t < V2 


x = 

1/8 

2/8 

3/8 

4/8 

5/8 

6/8 

7/8 

N = 

8 

.93124 

.86743 

.81382 

.77641 

.76303 

.78517 

.86121 

N = 

16 

.91927 

.83682 

.75243 

.66989 

.60198 

.58020 

.67304 

N = 

32 

.92285 

.83949 

.74438 

.63421 

.51243 

.40516 

.42378 

N = 

64 

.92533 

.84376 

.74900 

.63615 

.50275 

.35225 

.23835 

N = 

128 

.92598 

.84490 

.75039 

.63751 

.50369 

.3*1958 

.18274 

N = 

256 

.92614 

.84519 

.75074 

.63785 

.50398 

.34977 

.17939 

N = 

512 

.92618 

.84526 

.75083 

.63794 

.50405 

.34982 

.17940 

N = 

1024 

.92619 

.84528 

.75085 

.63796 

.50407 

.3*1983 

.17941 

CHAPTER  V 

FULL  TIME  AND  SPACE  DISCRETIZATION 


5.1  Opening  Remarks 


In  this  chapter  we  shall  introduce  a full  discretization  scheme  in 
the  context  of  an  abstract  boundary-value  problem. 

In  section  5.2,  the  procedure  of  full  discretization  will  be 
described.  The  relevance  of  the  order  in  which  the  schemes  are  applied 
will  also  be  discussed. 

In  section  5.3»  we  shall  state  our  main  convergence  result  in  the 
context  of  the  abstract  problem  (2.6). 

Finally,  in  section  5.H,  we  shall  apply  these  abstract  convergence 
estimates  to  the  canonical  parabolic  equation  (2.1)  and  we  shall  see 
that  optimal  rates  of  convergence  are  obtained. 

5.2  The  Abstract  Discretization  Scheme 


In  this  section  we  shall  see  how  the  scheme  is  constructed  as  a 
combination  of  the  two  semidiscrete  schemes,  i.e.  the  continuous-time 
and  the  discrete-time  approximation  schemes  as  described  in  chapters  III 
and  IV  respectively. 


-92- 


-93- 


A natural  way  to  formulate  the  fully  discrete  scheme  is  to  use  the 
differential  form  of  the  schemes,  more  precisely  the  starting  point  is 
the  weak  form  of  the  abstract  problem  (2.6)  as  stated  in  (2.7),  which  we 
recall  here: 

(u(t),v)H  = -(Au(t),v)H  + <g(t),B*v>u  for  all  ve  D(B* ) 

I u(0 ) = u 
v o 

One  way  to  proceed  is  to  consider  the  continuous-time 
discretization  scheme  applied  to  (2.7),  which  we  recall  here  as  stated 
in  (3.9): 

find  uh(t)  e Vh  such  that 


where  Ah:  Vh  approximates  the  operator  A. 

We  also  recall  the  discrete-time  approximation  scheme  as  stated  in 
(4.10): 


n = 0,1 


f • • • 


where  gn  = (Pkg)(nk),  and  un  is  the  approximate  of  u(nk). 
By  applying  the  later  scheme  to  (3-9),  we  obtain 
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l 

where  u^j^  approximates  uh(nk)  (and  consequently  u(nk)).  By  factoring 
uhk  + 1’>’  one  then  has 


given 


for  all  v.  e V. 

h h 


! 


(5.1 ) 


for  all  v,  e v, 
h h 


If  the  continuous-time  discretization  scheme  is  viewed  as  a 
projection  scheme  which  restricts  the  values  of  the  variable  to  a 
finite-dimensional  space  of  the  space  of  values  H,  then  the  process 
of  full  discretization  can  be  carried  out  in  reverse  order. 

In  fact,  by  applying  the  scheme  (4.10)  to  the  abstract  problem  (2.6),  we 
get 

j k^un+1-un^  = -Aun+1  + Bgn+1  in  n = °t1,... 

u given 

L 0 

where  gn  = (Pkg)(nk)  and  un  is  the  approximation  of  u(nk).  Solving  for 
un+1 » Wft  obtain 


in  D(A*) ' , n = 0,1 


9 • • • 


which  in  weak  form  becomes 
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\ ( (1+kA)un+1 ,v)H  = (un,v)H  + k <8n+1*B*v>u  for  all  veD(B«) 

/ u given 

L 0 

Then,  as  we  have  indicated  above,  the  continuous-time  discretization 
scheme  can  be  carried  out  (in  this  case  the  term  "continuous-time"  is 
very  inappropriate)  by  simply  restricting  the  values  of  un  and  v to  a 
finite-dimensional  space  Vh£D(B*)  and  by  replacing  A with  its 
approximation  A^:  Vh  -*■  Vh.  We  then  obtain  exactly  the  scheme  (5.1). 


This  scheme  is  carried  out  in  practice  by  on  the  one  hand  selecting  a 
certain  number  M of  basis  functions  for  the  space  V^.  This  number 
determines  the  value  of  the  parameter  of  discretization  h (usually  -). 
A popular  choice  is  to  select  basis  functions  with  local  supports  in 
such  a way  that  the  inner-product  matrices  (the  so-called  mass  and 
stiffness  matrices)  have  a band  structure  or  at  least  a sparse 
structure.  These  functions  form  what  is  called  a finite  element 
basis. 

On  the  other  hand,  one  has  a finite-difference  scheme  such  as  the 
weak  formulation  of  (4.10),  i.e., 


(0+kA)uk((n+1)k),v)H  = (uk(nk),v)H  + k <Pkg(  (n+1  )k)  ,B*v>u 


Vv  e D ( B* ) , 


uk(0)  given. 


As  we  have  just  seen,  the  order  in  which  the  semidiscrete  schemes  are 
applied  is  immaterial  and  actually  has  little  meaning.  It  is  more 
appropriate  in  this  situation  to  view  the  fully  discrete  scheme  as  a 
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juxtaposition  or  simultaneous  application  rather  than  a successive 
application  of  the  two  semidiscrete  schemes.  After  assembling  the  dif- 
ferent matrices  and  vectors,  one  is  faced  with  a sequence  of  linear 
systems  of  algebraic  equations  to  be  solved,  which  in  theory  is  an 
elementary  task  (but  not  always  in  practice!).  The  finite  element 
method  is  well-known  for  its  advantages  over  standard  finite  difference 
algorithms  in  most  cases,  such  as  increased  accuracy  and  flexibility  of 
application.  We  shall  use  the  finite  element  method  in  chapter  VII, 
where  it  will  be  carried  out  in  the  canonical  case  of  the  heat  equation. 

Remark . In  the  case  of  a Galerkin  approximation,  the  single  step 
of  restricting  un  and  v to  the  subspace  Vh  is  sufficient  since  by 
definition  of  the  Galerkin  approximation  A^, 

// 

(Ahuh.vh)H  = (Auh*vh^H  Vuh»vheVh- 

5.3  Estimation  of  the  Error 


Whereas  the  formulation  and  numerical  implementation  of  the  fully 
discrete  approximation  are  most  naturally  carried  out  by  using  the 
^ ffsrenti al  forms  of  the  problems,  it  will  be  to  our  advantage,  as  we 
shall  see,  to  estimate  the  error  of  discretization  by  examining  and 
comparing  the  semigroup  formulas  corresponding  to  the  different 
discretization  schemes. 

Our  first  goal  in  this  section  is  to  derive  a (discrete)  semigroup 
formula  for  the  fully  discrete  approximation  to  the  abstract  problem 
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(2.6).  As  in  section  (5.2),  such  a formula  can  be  derived  in  two  ways, 
which,  as  we  have  pointed  out  in  the  previous  section,  are  merely  two 
ways  of  describing  the  same  idea. 

As  for  the  first  way,  we  recall  from  section  3.2  that  the  semigroup 
formula  for  the  continuous-time  approximation  of  problem  (2.6),  given  by 
(3.10),  is 

t 

Uh(t)  = J"  ShU-T)Ph(Bg(-r))dt. 

The  discrete-time  algorithm  applied  to  the  solution  of  problem  (2.6)  can 
be  also  applied  to  uh(t)  without  modification,  and  thus  we  obtain  a 
formula  similar  to  (4.11),  i.e.: 

(5.2)  uhk(nk)  - EjkPhuo  ♦ kj  C,'JphB<Pk8><Jk)’  "-'.2.- 

where  Ehk  = 1 r(1,  -a^)  and  Ah  is  an  approximation  of  A.  As  in  (4.11), 
uhi<(t)  is  defined  in  the  obvious  way  for  the  other  values  of  t. 

The  second  way  to  derive  a formula  for  the  fully  discrete 
approximation  uhk  is  to  start  with  the  formula  for  the  discrete-time 
approximation  uk  of  the  solution  to  problem  (2.6),  i.e. 

uk(nk)  = Ek  uq  + k l E^+1  ^B(Pkg) ( jk) , n = 1,2,... 

J ^ 

The  fully  discrete  approximation  uhk  is  then  obtained  by  applying  the 
orthogonal  projection  Ph  to  the  input  B(Pkg)(jk)  and  to  the  initial  data 
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uQ  and  by  replacing  A with  its  approximation  Ap.  It  can  be  readily  seen 
that  this  procedure  yields  precisely  (5.2). 

We  are  concerned  with  the  estimation  of  u-upk  in  the  norms 
Lp[0,T;H]  for  input  g e Lp[0,T;'J] , p = 1 , To  this  effect  we  define 
the  operator  Thk:  Lp[0,T;U]  -►  Lp[0,T;H]  by 

n i • 

(5.3)  Thkg(t)  = uhk(t)  = k I C "JPhB(Pkg)(jk),  n = 1 ,2 

j = 1 

where  we  assume  without  loss  of  generality  that  uQ  = 0. 

It  is  clear  that  Thk  is  well-defined  in  the  topology  of  Lp[0,T;H]  for  g 
in  L>p[0 ,T; U] , and  this  for  the  same  reasons  given  for  Tk. 

Since  we  shall  estimate  the  error  of  discretization  with  the 
adjoints,  we  need  to  derive  the  formula  for  Tpk.  We  have 

T 

<Thk8’V,L2[0,T;H]  ■ J'Thk8<t).'-(t))Hdt 

T 

j[(Thkg(t)  ,Phv(t)  )Hdt 

since  Tpkg(t)  eVh  for  all  t > 0.  From  this  point  on,  the  derivation  of 

a _ % v 

Thk  follows  exactly  that  of  Tk.  Thus  the  formula  for  Tpk  is 
(5-4)  Thkphv(Jk)  - kB*  n^Ehk+1'Jpk((phv)(nk))*  J - 1.2,... 
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Reroark . The  order  in  which  and  Pk  appear  in  the  formulas  (5.3) 
and  (5.4)  is  irrelevant  since  Ph  and  Pk  commute.  Therefore  we  can 
extend  T^k  to  all  of  Lp[0,T;H]  in  the  obvious  way. 

// 

We  are  now  in  a position  to  state  a basic  convergence  result. 

Theorem  5.1.  For  0 < h < 1,  let  Vh£D( B* ) be  a family  of  finite- 
dimensional subspaces  of  H satisfying  properties  (3.1)  - (3.4),  and  let 
Ah  be  a family  of  approximations  of  the  operator  A satisfying  (3.5)  and 
(3.6).  Let  T and  T^k  be  defined  by  (2.37)  and  (5.4)  respectively. 

Then 


|(T*-T*k)v| 


Lp[0,T;U] 


,.20-qL  1 

c(h  ^ In  - 
h 


-2q 


k m -)|v 


Lp[0,T;H] 


for  p = 1,  ®,  where  q £ q and  A qB e 1 (U+H)  for  some  0 < q < 1.  The 

parameters  of  discretization  h and  k are  subject  to  the  restriction 
2 

k < ch  with  c > 0 arbitrarily  chosen. 


// 

2 

The  restriction  k ^ ch  can  be  removed  by  imposing  stronger 
approximation  properties  on  the  subspaces  . In  fact,  we  have 


Theorem  5.2.  For  0 < h < 1 , let  Vh£D(A*)  be  a family  of  finite- 
dimensional subspaces  of  H satisfying  the  properties  (3.1)  - (3.3)  and 
(3.4'),  and  let  A^  be  as  in  Theorem  5.1  a family  of  approximations  of 
the  operator  A satisying  (3.5),  (3.6).  Let  T*  and  T£k  be  defined  by 
(2.37)  and  (5.4)  respectively. 
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Then 


|(T*-T*k)v|  < c(h2(1_q)  + k 1 ~ Q ) I V | 

Lp[0,T; U]  Lp[0,T;H] 


for  p = 1 , 00 , where  A qB  e ^(U-H)  for  some  q,  0 < q < 1 . There  is  no 
restriction  on  k or  h. 


// 

The  estimates  stated  in  Theorems  5.1  and  5.2  are  also  valid  by  duality 

for  T-Thk.  in  fact,  from  Theorems  5.1  and  5.2,  it  follows  immediately 
by  Lemma  2.*)  that 


<5'5>  |(T-VEI,  roT  ^ * K^lnl^gl 

Lp[°,T;H]  Lp[0,T;U] 

for  p = 1 , ®,  by  making  use  of  the  restriction  k £ ch2  and  with  the 
assumptions  of  Theorem  5.1.  Also,  we  have 

<5-6)  |(I-T  )g|  S cOi2'1^  . 

Lp[°.T;H]  Lp[0 ,T ;U] 

for  p = 1,  ®,  with  the  assumptions  of  Theorem  5.2,  and  without 
restriction  on  k or  h. 

We  now  make  the  following  remarks. 


(i)  The  estimates  stated  in  Theorem  5.1  and  5.2  can  be  extended  to 
the  case  1 < p < ® by  interpolation. 

(n)  In  Theorem  5.1,  the  subspaces  Vh  are  taken  in  D (B* ) , whereas 
in  Theorem  5.2,  they  are  required  to  lie  inside  D ( A* ) . This  latter 
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condition  is  clearly  stronger  and  is  necessary  in  order  to  prove  the 
auxiliary  result  A*A£_1Phe  £ (H+H)  (see  section  6.6)  on  which  the 
estimation  of  T*-T^k  is  based.  In  the  next  section  we  shall  see  what 
this  stronger  condition  means  in  the  context  of  the  canonical  parabolic 
problem  (2.1). 


The  abstract  scheme  (5.1)  and  its  corresponding  solution  (5.3),  as 
well  as  the  results  (5.5)  and  (5.6),  can  be  specialized  directly  to  the 


B = AD  (formally)  with  A ^(U+H)  for  q = | + e,  e>0,  and  q = | . 

Also  the  abstract  scheme  (5.1),  specialized  to  this  particular  context, 
takes  the  form 


5.'4  Application  to  a Canonical  Example 


case  of  our  parabolic  problem  (2.1)  by  taking  H = L2(fl),  U = l2(D, 


l2(q) 


A L2(D 


n = 0, 1 


9 • • • 


(5.7) 


for  all  v,  e V, 
h n 


In  fact,  as  our  main  result  we  have  the  following  theorem 
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Theorem  5.3.  Let  Vh£D(D*A*)  be  a regular  system  satisfying  the 

inverse  assumption,  as  well  as  properties  (3.1*0  - (3.16).  Also  let  Ah 
be  an  arbitrary  approximation  of  A satisfying  (3.17)  - (3.20).  If  y(t) 
is  the  solution  of  the  parabolic  problem  (2.1)  and  yhk(t)  is  the 
solution  of  the  fully  discrete  scheme  (5.7),  then  e(t)  = y(t)  - yhk(t) 
satisfies  the  inequality 


Similarly,  we  have  corresponding  to  Theorem  5.2  the  following  result. 

Theorem  5.**.  Let  V^.A^  satisfy  the  same  assumptions  as  in  Theorem  5.3, 
with  the  difference  that  we  now  require  that  Vh£D(A*).  Also,  let  y(t) 
be  the  solution  of  the  parabolic  problem  (2.1)  and  yhk(t)  the  solution 
of  the  fully  discrete  scheme  (5.7).  Then  e(t)  = y(t)  - yhk(t)  satisfies 

the  inequality 


£ c(h  ^ 


In  - + k 
h 


Lp[0,T;L2(n)] 


Lp[0,T;L2(r)] 


2 

for  1 ^ p < ® with  the  restriction  k < ch  for  c > o arbitrarily  chosen. 


Lp[0,T;L2(fl)] 


S c(hV2-2£  ♦ k Vi*  ~E)  [g| 


Lp[0,T;L2(r)] 


for  1 < p < oo?  an(j  without  any  restriction  on  k or  h. 


// 


We  now  make  the  following  comments. 
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(i)  In  view  of  the  best  approximation  properties  of  and  the 
maximal  regularity  of  the  solution  y to  the  parabolic  problem  (2.1),  the 
estimates  0(h  ^ In  and  0(k  ^ In  are  almost  optimal  (modulo  In  ^ 
and  In  — ) for  1 < p < °°,  and  are  optimal  for  p = 1 , 

(ii)  If  we  recall  the  rates  of  convergence  obtained  for  each 
semidiscrete  approximation  scheme  (0(h1/2ln  and  0 (k  ^ 1/lJ  ~C  ^ ) , we  see 
that  in  Theorem  5.3  essentially  a power  e of  k has  been  regained  in  the 
estimates . 

The  slight  difference  in  these  rates  of  convergence  reflects  the  methods 
used  in  proving  each  estimate.  In  the  case  of  the  fully  discrete 
approximation,  use  is  made  of  the  inverse  approximation  property  (3.4) 
of  the  subspaces  to  "kill"  the  e in  0(k^  ^ e^).  This  improvement  is 

obtained  at  the  price  of  restricting  the  values  of  the  parameters  of 
discretization  h and  k by  the  condition  k £ ch2. 

As  we  have  seen,  this  restriction  can  be  removed  by  imposing  stronger 
approximation  properties  on  the  subspaces  Vh.  Unfortunately,  the 
optimal  rates  of  convergence  0(h1/2ln  ^)  and  Otk^ln  ±)  are  lost  in  the 

process  and  we  are  left  with  0(h(1/2_2e))  and  0(k(1/4_e)).  In  the  case 
p = 1,  ®,  these  latter  rates  of  convergence  can  nevertheless  be  viewed 
as  essentially  optimal  since  the  regularity  of  the  solution  of  (2.1)  is 

hv2-2E,  vE(QK 

From  the  point  of  view  of  numerical  computations,  the  restriction 
2 

k < ch  is  far  more  significant  than  the  loss  of  "e"  in  the  rates  of 
convergence.  In  the  case  of  Theorem  5.4  another  significant  restriction 
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is  Vh£D(A*).  It  is  known  that  linear  splines  defined  on  uniform  meshes 
are  in  D(D*A*)  but  not  in  D(A*4)  for  q > Therefore,  whenever  linear 
splines  are  used  Theorem  5.4  does  not  apply  and  the  restriction 
k S ch^  is  a necessity. 

As  for  the  case  1 < p < <*>,  the  rates  of  convergence  0 ( h ^2  ) and 

0(k(  4|  e))  follow  immediately  by  examination  of  the  proof  of  Theorem 

2 

5.1  anci  the  restriction  k £ ch  still  applies.  These  rates  are  almost 


optimal  (modulo  c)  in  view  of  Theorem  2.2  (for  1 < p < the  regularity 

p 


of  the  solution  is  H 1/2  ’ ^ (Q)). 


CHAPTER  VI 

PROOFS  OF  THEOREMS  3-1,  4.2,  5.1,  5.2,  AND  LEMMA  4.1 


6.1  Opening  Remarks 


In  this  chapter  we  shall  present  the  proofs  of  the  main  theorems 
which  have  been  stated  in  this  work.  Each  of  these  results  states  the 
approximation  properties  of  the  operators  corresponding  to  one  of  the 
schemes  we  have  defined.  Except  for  Theorem  3.1,  which  was  established 
in  Lasiecka  [12],  all  the  theorems  proved  in  this  chapter  are  original. 
Also,  the  proofs  contained  in  this  chapter  are  all  in  the  context  of  the 
abstract  problem  (2.6). 


6.2  Proof  of  Theorem  3.1 

We  recall  (2.37)  and  (3.11 ). 

T 

T*v(t)  = B*  J S*(i-t)v(T)dt  , 

t 

T 

v(t)  = B*  j S*(T-t)Phv(-r)dT  . 
t 

For  the  purpose  of  estimating  (T*-T£)v(t),  we  let 
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(T*  -T£)v(t)  - e*  (t)  + e*(t) 


where 


e*(t)=  B*(I-Ph)  J S*(x-t)v(x)dt 
t 


and 


e*(t)  = B*  J[PhS«(T-t)-S*(i-t)Ph]v(T)dt. 

We  first  condsider  the  case  p = ®.  For  t < T-h14 , we  have 

t+h^ 

|e*(t)|  < |B*(I-Ph)  j S*(x-t)v(x)dx| 


|B*(I-P  ) j 4 S*(x-t)v(x)dx| 
t+h4 


By  applying  (3-3)  to  the  first  term  with  a =Vn-2e,  and  (3.2)  to  the 
second  term,  we  get 

t+h4 

|e*(t) | S ch2(1_q~e)|  f A*1_eS*(x-t)v(x)dx|H 


+ Ch2(1 "q)  | / A*S*(T-t)v(T)dT| 


t+h 


Then,  by  analyticity  of  the  semigroup  S*(t), 


t+h 


|e*(t)|u  s oh2(1-<i-£>  / 


1 


W dT  M 

t (x-t)  L [0 , T ;H ] 


L [0,T;H] 
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S ch2(1'<1‘£,h1'£|y| 


. Oh21'-1!)  1 n 'I 


oh"'  ■>'  in  -j|v  | 

L C0,T;H]  h L 


< c(h2(1"q+e)+  h2(1_q)ln  l)|v| 


L [0,T;H] 


Similarly,  for  t 2 T-h  , we  have  by  application  of  (3.3) 
with  a = Vi|  -2c, 

T 

le*(t) lu  = |B*(I-Ph)  J S«(T-t)v(T)di|u 


< ch2(1  q_e)|  / A*1  ES*(T-t)v(T)dT|H 


S oh2(,-<l-£)  ; — d-r  |v  I 

t (T-t)  L [0,T;H] 


< Ch2(1"q"e)(T-t)£|vl 


LCO.TjH] 


S Ch2(1"q+e)|v| 


Lm[0,T;H] 


Thus  we  have  shown  that 

le*  I S o(h2(1-<ltEh  h2^  1),V| 

L^O.TjU]  LJO.TjH] 

We  now  consider  the  term  e*(t).  By  (3.4),  we  have 
T 

N*(t)|  = |b*  /(PhS*(T-t)-S*(T-t)Ph)v(T)dT|u 
- T 

s Ch_2q  / |PhS«(T-t)-S*(T-t)P  | dt |v | 

1 H+H  L^O.TjH] 


[0  ,T; H ] 

00 
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In  order  to  estimate  the  integrand,  we  shall  need  to  use  a result 
derived  in  Lasiecka  [11]. 


Proposition  6.1  [11].  With  A, A*  satisfying  (3.5),  (3.6),  we  have 


|Sh(t)Ph-PhS(t)|„,„  S c 


21 
i 

£ ’ 


0 < £ < 1 . 


// 


For  t < T-h  , we  have 


t+h 


le2(t)lu  S ch"2<1  I |PhS*(T-t)-S*(,-t)Ph|HiHd,|v| 


* ch'2q  / _2|PhS*(T-t)-S«(,-t)Ph|H<H<J,|v| 


t+h 


L [0,T;H] 


L_[0 , T; H ] 


By  applying  Proposition  6.1  with  £ = 1-e  to  the  first  integral  and  £ 
to  the  second  integral,  we  obtain 

2 


= 1 


t+h 


2-2e 


le2(t)lu  s oh'2q  [ / / 2T^<1t  ]|v| 

t (x-t)  t+h  L^LO.TsH] 

<c(h2(1~q~e)h2e  + h2(1"q)in(Ilt))|v| 


h L [0,T;H] 


< c(h2(1"q)+  h2(1  q)ln  1) U v | 


L [0,T;H] 


2 

For  t Z T-h  , we  also  apply  Proposition  6.1  with  £ = 1-e  to  get 


2-2  e 


|e.(t)|u  S ch"2q  / -5-^  dr  |v  | 

t ( T-t ) L [0,T;H] 


< ch2(1  q_e)(T-t)£|vj 


L [0,T;H] 


* ch2(1'q) 


M 


L [0,T;H] 


Thus  we  have  shown  that 


|e*|  S c(h2(,-c').  h2(,-q)ln  JL)|»| 

t-„[0,T;U]  LJO.T 

and  this  completes  the  proof  of  Theorem  3.1  for  p 
As  for  the  case  p = 1 , 

T T T 

/ le*(t) Ldt  = / |B*(I-P  ) / S*(T-t)v(T)dx|  dt 
0 0 n t U 

T T 

* J /|B*(I-Ph)S*(T-t)v(T)|udTdt 

0 t 

By  reversing  the  order  of  integration, 

T T t 

/ |e*(t)  | dt  < J / |B*(l-p  )s*(i-t)v(T) LdtdT 
0 u 00  u 
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S j / |B*(I-P  )S*(T-t)v(T)|MdtdT 
0 0 u 


+ Jh  ) |b*(i-p  )s*(T-t)v(x) LdtdT 

h 0 
T T 


+ Jj,  J jjB*(I-Ph)S*(T-t)v(T)|  dtdx 
h T-h 

We  note  that  in  order  to  estimate  the  first  integral  we  shall  use  the 


/ le^(t) lydt  £ ch2(1  q_e)  J /|A*1'eS*(T-t)v(T)|HdtdT 


fact  that  0 < t ^ h^ . 


Also  the  bound  x-h  in  the  second  and  third 


integrals  is  well-defined  since  h^  < t i T.  By  applying  (3.3) 

with  a = % - 2e  to  the  first  and  third  integrals,  and  (3.2)  to  the  second 

integral,  we  have 


T 


h 0 


T T 


/l,  Jh4|A*1_eS*(T-t)v(T)|HdtdT 


From  the  analyticity  of  S»(t),  it  then  follows  that 
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h4  T 


1 1 e* ( t ) Ijjdt  < ch2(1'q'e)/  (J 


dt 


0 0 (T-t) 


— ) |’(T)  |HdT 


* 0h2l'~q>J„  <J  7^)|'(T)|HdT 

h 0 T L M 


T T 


* 0h2l1"q"E)  l"'1  h-  — 5Tw)lv<t)lH^ 


Hence 


T hH 

j|e*(t)|  dt  < ch2(1"q"e)  J iefv(T)|Hdx 
0 u o H 

_ T 

* ch2('‘q)J11  ln(\)|v(T)jHdt 


’ll  lT 

h h 


* Oh2<1-‘l-C)  /4  h"E|v(t)|HdT 


< c(h2(,-^E)*  h^'-o’ln  1)|»| 


L,  [0 ,T;H] 


We  now  estimate  |e$l  • By  (3-4), 


L^O.TjU] 


T T 


/ Ie| ( t)  | Udt  = JlW(PhS»(T-t)-SS(T-t)Ph)v(T)dTjudt 


T T 


Ch  2q  / /|(PhS*(T-t)-S*(T-t)P.  )v(T)|udTdt 
0 t n n *H 


By  reversing  the  order  of  integration, 
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T X 


/ Ie* ( t)  lydt  S ch_2q  J /|PhS*(T-t)-S»(T-t)Ph)v(T)|HdtdT 


0 0 


l-l2  T 


* Ch_2q  J J|(PhS*(T-t)-S*(T-t)P  )v(T)|„dtdT 
0 0 H 


_ T x-h 

+ Ch_2q  j2  I |(PhS*(T-t)-S*(T-t)Ph)v(T)|HdtdT 


_ T T 

+ Ch_2q  l2  I 2|(PhS*(T-t)-S*(T-t)Ph)v(T)LdtdT 
h x-h^ 


By  applying  Proposition  6.1  with  l = 1-e  to  the  first  and  third 
integrals,  and  l = 1 to  the  second  integral,  we  obtain 


h2  t 


2-2e 


J|e*(t)|  <Jt  S oh'2q  / J -h.  dt)|v<T)|Hdx 

00  (T-t)1  e H 


T x-h 


* oh'2q  J2</  dt)|v(T)|  at 

h'  0 T H 


T T 


2-2e 


* ch"2q  /2  (J  2 — dt)|v(T)|  dl 

h z-n  (T-t)'  e H 

h2  t 

< ch2(1_q'e)  f 


/ TE|v(T)|HdT  + ch2(1_q)|  ln(^)|v(x)|  dx 

0 h h M 


+ ch2(1-q_e)  j h2e|v(x)|  dx 
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£ c(h2(1  q)+  h2(1"q)ln  1) |v| 

L1 [0 , T; H ] 

% £ £ 

Since  (T*  - T^)v(t)  = e -j  ( t ) + e2(t),  this  concludes  the  proof  of  Theorem 

3.1  • 

6.3  Proof  of  Lemma  4.1 


We  shall  make  use  of  an  integral  representation  for  T(t)  and  eH 

K 

(see  Pazy  [18,  pp.  32-331). 

Let  r = r^U  r2,  where 

r1  = {re1(J>,  0 < r < »)  and  r?  = {re_i<t),  0 < r < »}  , with  <j>  as 
stated  in  Lemma  4.1.  Then 

(6.3.1)  T(t)  = J-  |etAR(A,A)dA 

r 

and 

(6.3.2)  Ej  . ((I-kA)-')n  - j-  / 1_  R( X , A )dX 

r (i-kA)n 

In  order  to  establish  (i)  we  shall  need  to  prove  the  following  technical 
lemma. 


Lemma  6.2.  Let  z = x+iy  be  given  such  that  |z|  =1,  x < 0.  Then 


/ _ \ i niz  1 , z 2 i 

le  - cxnT  e for  0 < t < -r—r,  n > 1 


nix 
2 2 


(1-tz) 


ITT’ 


(b) 


h-xz|n 


- S e 


n 

2 


for  t -j — r,  n 2 1 

7*1 
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The  constant  cx  depends  on  x but  not  on  n or  t. 


// 


Proof  of  Lemma  6.2.  To  prove  (a)  we  first  need  to  prove  the 
following  inequality 


(6.3.3) 


TZ 
2 , 


1 -TZ 


for  0 £ t ^ 


w 


TZ  TX 


Since  1 1 — tz I = /i+t2-2tx  and  |e2  I = e2 


= e , it  is  equivalent  to  show 


(6.3.4)  eTX(1+T2  - 2tx)  > 1 


for  0 S t = 


It  is  clear  that  (6.3.4)  is  valid  for  t = 0.  It  is  therefore  sufficient 

to  show  that  the  derivative  in  t of  eTX( 1 +t2-2tx)  is  positive  on  the 
interval  [0,  1/|x|].  We  have 


^ (eTX(1+T2-2Tx) ) = eTX(x(1+T2-2Tx)  + 2t-2x) 


It  is  enough  to  study  the  sign  of 

(6.3.5)  x(1+t2-2tx)  + 2t-2x  = xt2  + 2(1-x2)t  - x. 


As  a quadratic  polynomial  in  t,  the  roots  of  (6.5)  are 
-(1-x2)  + /(1-x2)2+x2 


x 


At  this  point  we  make  the  following  observations: 


(i)  Since  - 1 < x < 0,  the  root 


-(1-x2)+/(1-x2)2+x2 


1-x 


[-1  ♦ 


1+- 


/ . 2.2' 
(1-x  ) 


is  clearly  negative. 

(ii)  the  coefficient  of  t is  negative,  which  makes  the  graph  (a 
parabola)  concave  downward. 

(iii)  the  other  root  satisfies 


-(1-x2)- 


77.  272  2 

/(1-x  ) +x 


- i [i-x2  ♦ A 


2,2  2, 
(1-x  ) +x  ] 


* ]4| (i_x2  + m ^ 


> 


T^T 


from  which  (6.3.3)  follows. 


The  next  step  in  proving  (a)  is  to  show 


// 


(6.3.6) 


xz 


1 -TZ  1 


C X 
X 


for  0 < x £ 


1 

w 


First  of  all,  we  note  that 
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Since  | 1 - tz | 


-2tx  £ 1 , 


it  is  enough  to  show 


(6.3.7) 


|eTZ(1-Tz)-1 


- V 


for  0 ^ t i 


w 


Using  a power  series  expansion  for  e 


•tz 


we  have 


eTZ(1 -tz) 


tz  TZ 

e - Tze 


1 


= 1 +TZ 


2 2 
T Z 

2! 


3 3 

T Z 

+ + 

3! 


n n 
t z 

n! 


2 2 

- 1 - TZ  - TZ(TZ)  - TZ  T'  - -- 


TZ 


n-1  n-1 


T z 
(n-1  )! 


- I 

k=2 


k k 
t z 


1 

( k- 1 ) I 


) 


v k k. 

= L t z (- 

k=2 


k-1 

k! 


In  modulus, 

TZ 


e ( 1 -tz)-1 


* l T 
k=2 


k k-1 
k! 


< 


I 


oo 


l 

k=1 


k 

T 

k! 


£ T(eT-1). 


We  have 

eT-1  ^ c t 
x 


for  0 < 


T 


< 
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where  of  course  cx  depends  on  x. 

Thus  (6.3.6)  follows. 

We  are  now  in  a position  to  conclude  the  proof  of  (a), 
n i 1 , we  have 


For 


mz  1 1 , xz  1 


n-1 


= e 


(1-TZ)1 


J (n-1-j )tz 1 

1_TZ"j=0  (1-tz)J 


I 


* |eTZ-T^l“n— 3~Tlle(n'1“J)TS| 

1 TZ  j =0  (1-TZ)J 


Since 


( 1 - TZ ) J | 1 — TZ 


• and  by  applying  (6.3.3)  and  (6.3.6), 


ntt  1 


s c ,2  V |e 


/.  sn  1 x . _ 
(1-TZ)  j =0 


n-1 


jxx 


2 ‘V  2 (n-l-j)TX 
Set  I e e J 

j=0 


n-1  "IH 

C „ 2 V (n-1  )TX  2 

< c t I e e 

j-0 


For  0 ^ j ^ n-1  , and  since  x < 0,  we  have 


j tx  _ n^J 

a 2 2 
e < e 


-tx 
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So 


, jtx  n-1 . 

2n;'  (n-i  )tx  - ~r<  2n;'  (n'1'— ),x 

c x ) e e S c x ) e 

X j =0  X j =0 


n-1 

, 2 2 
Sc  t n 8 
x 


-TX 


mx  _ TX 

2 


. 2 2 
S c x n e e 
x 


mx 

. 2 2 
S c nx  e 
x 


for  0 S x S 


■RT 


This  concludes  the  proof  of  (a).  As  for  (b),  we  have 

-n/2 


1 


1 


i in  2 . .n/2  ...  2 _ . -1 .n/2 

1-tz  (1+t  -2tx)  ((1+x  -2xx)e  ) 


, n = 1,2,... 


It  is  therefore  sufficient  to  show  that 


(6.3*8)  1+x  - 2xx  £ 


for  x Z 


1 

171 


We  now  make  the  following  observations 


(i)  X < |x|  < 1 < -jly 


(ii)  The  graph  of  y 


+ x -2xx  is  concave  upward  and  has  a minimum 


at  x 


x. 
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(iii)  It  follows  from  (ii)  that  1+t  -2tx  is  increasing  for  t £ x,  and 


thus  by  (i),  also  for  t Z 


T*T  • 


(iv)  (6.3.8)  is  valid  for  1 = -ry-r  since  1 +— y+2^4>e  from  which 

' ' x 

(b)  follows. 


// 


We  now  return  to  the  proof  of  Lemma  ^4 . 1 . 

From  (6.3.1)  and  (6.3.2),  we  have 

T(nk)  - e?  . ' f(e"k»  - — ! ) R ( A , A ) d A 

k 2,1  r ( 1 — k A ) n 

where,  as  we  recall,  r = r U T2  with 

f1  = (re1**1,  0 £ r < «>}  and  f2  = {re  i<*>,  0 S r < «} 


and  with  <p  as  given  in  the  statement  of  Lemma  4.1.  Thus  if  we  let  z 
id) 

= e = x+iy  and  X = rz,  then 


T(nk)  - e" 


CO 


1 

2iri 


r , nkrz  1 , . , 

J (e  - )R(rz,A)zdr 

0 (1-krz)n 


CO 


1 

2iri 


[,  nkrz 

J (e 

0 


— — )R(rz,A)zdr 

( 1 -krz)n 


From  the  analyticity  of  T(t),  we  can  apply  (2.3).  Hence, 
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|T(nk)-E"|  * = f |enkrz 


1 


( 1 -krz ) 


n 1 1 + rz 


dr 


+ c J |e 
0 


nkrz  1 


(1-krz)n  1 + |rz 


dr 


kJJT 

S o j [|enkrz 


1 


( 1 -krz ) 


+ e 


nkrz  1 


(1 -krz) 


1 ] -r- — dr 

- n 1 1 +r 


+ c j [ | e 
kj^| 


nkrz 


1 | | nkrz  1 i,1 

+ e - ] - — dr 

1 1 -.n1  1+r 


>n 


(1-krz)  (1-krz) 

The  next  step  will  require  the  following  inequality: 


(6.3.9) 


1 


1 


|l-Tz|n  |l-Tz|n"e  |1-TZ|' 


n-e 

1 2 1 
^ e . — 


for  t t -|Y|-,  n 2 1 , and  e>0 

which  can  be  readily  seen  since  1 1 - xz | i | tz | = t.  Therefore,  by  using 
Lemma  6.2  with  x = kr,  we  have 


1 


k x nkrx 


|T(nk)-E"|  S c / nk2r2e  2 1 dr  . o J enkrx  I Hr 


(n-e) 


+ c j e 2 
1 


i e e 

k r r 


dr 
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1 

k | x j nkrx 

(6.3.10)  < cnk2  j re  2 dr 

0 

® n-  e °° 

r nkrx  1 , 2 , -c  r 1 

+ c J e - dr  + ce  k — — dr. 

1 r 1 r +e 

*1*1  i<T^T 

If  we  set  s = -krx  = kr|x|,  then  for  the  first  integral  in  (6.3. 10),  we 
have 


k | x | nkrx 

cnk2  j re  2 dr 

0 


cnk2  j 


0 *T*T  6 


ns 

2 1 


k x 


ds 


1 ns 
2 


S cn  Jse  ds 
0 


ns  1 ns 

S =„<-  u e 2 | * f J e~  2 ds) 
0 0 


n 


n 


* 2^f(-f(e  2-  <»> 


n 


n 


£ £(-2ne  2 - Me  2 -1)) 


< £ 


, n > 1 


Similarly,  for  the  second  integral  in  (6. 3. 10),  we  have 
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00 


c 


I 


00 


nkrx  1 

e — dr  = c 
r 


ds 


00 


00 


n 


Finally,  for  the  third  integral  in  (6.3.10), 


S ce 


n-e 

2 


00 


/ 


ds 


£ ce 


n 

2 


< 


c 

n‘ 


This  concludes  part  (i)  of  Lemma  4.1. 

In  order  to  prove  part  (ii),  we  shall  consider  separately  AE^  and 


AT(nk) . 
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First  of  all,  it  is  clear  from  the  analyticity  of  T(t)  that 


(6.3.11)  (filCnk, |H^H  < 


nk' 


For  AE^ , we  shall  consider  the  cases  n = 1 and  n > 2.  For  n = 1,  we 
have 


I = (l-kA)d-kA)'1  = Ek  - kAEk, 
so  that 

kAEk  = Ek  ~ I e & (H  + H) 
and  thus  we  have  established  (ii)  for  n=1  . 


For  n > 2,  we  make  use  of  the  integral  representation  for  AEk , i.e., 


AE 


n = — — [ 
k 2iri  J 


X R ( X , A)dX 

( 1 -kX )n 


Similarly  as  in  part  (i),  we  let  z = e1^  = x+iy  and  X = rz.  Then 

CO  a, 

AEk  = dd  J”  rz  n R(rz*A)zdr  + ^ / rz  R(rz,A)idr. 

0 (1-krz)n  2irl  0 (1-krz)n 


Hence 

00 


lAEk  I s 0 I [- - — — |R(rz,A)|  + |R (rz , A)  | ]dr 

0 j 1 -krz | n 1 1 -krz | n 


Again  by  analyticity  of  the  semigroup  generated  by  A,  and  by  application 
of  the  inequality  (6.3-9),  we  deduce  that 
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1 

k I x I nkrx 


n- 1 - e 


|AE‘k‘|  < c / e 2 dr  + c J — A + £. 


dr  + c 

0 nk  n-tl  -e  (kr ) ' ^ 

S Ue  2 -1)  * «’  ^ k-1-c  j 

nk  J i+e 

1 r 


1 


dr 


dr 


k x 


2 ~r~ 

nk 


for  n k 2. 


Therefore 


|A(E"-  T(nk) ) | S |AE"|  * |AT(nk)| 


lH»H  " 1 k’H+H 
c 

nk’ 


H+H 


n > 1 , 


which  concludes  the  proof  of  Lemma  h . 1 . 


6.4  Proof  of  Theorem  4.2 

We  have  the  following  decomposition  for  (T*-T*k)v(t),  where 
(j-1  )k  < t 2 jk. 

(T*-T*)v(t)  = T*v(t)  - T*(P  v)(t) 

+ T*(Pkv)(t)  - T«(Pkv)((j-1)k) 

+ T»(PRv)( (J-1 )k)  - T*k(Pkv)(jk) 

Remark . The  apparent  inconsistency  in  the  last  line  of  the  above 
expression  ((j-1)k  vs.  jk)  is  purely  a matter  of  notation.  In  fact,  one 
can  see  from  the  definition  of  Pk  (see  (4.8))  that  T£(Pkv)(jk)  does 
depend  on  the  values  of  v over  the  interval  ((j-1)k,T]. 

// 

We  have  the  further  decomposition 
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(T*-Ti/  )v  (t ) = T*(I-Pu)v(t) 


3*t  j S*(Tt)Pv(i)dT  - j S*  (t~  ( j-1  )k  )P  v ( t ) d t ] 


(j-1  )K 


nk 


+ B*  i [ / (S*(t-(j-1  )k)  - S*((n+1-j)k))P  v(t)c!t] 


n=j  ( n— 1 )k 


nk 


+ B*  1 [ j S*( (n+1 -j ) k)P  v ( t )dx 
n=j  (n-1)k 


kE*n+1  jP,  v(nk)] 
k k 


At  this  point  we  observe  that 


nk 


nk 


j S*(  (n+1-j  )k)P  v(tMt  = S*((n+1-j)k)  J P v(t)<1t 


(n-1 )k 


(n-1  )k 


= k S*( (n+1 - j )k)P^v(nk) 


Therefore, 

(T*-T£)v(t)  = T*(I-Pk)y(t) 


" B * j S*(t-(j-1 )k)P  v(x)dT 


(j-1  )k 


- B*(S*(t-(j-l)k)-l)  j S*(i-t)P  v(i)dT 


N 


nk 


+ B*  l i (S*(i-(j-1)k)  - S*((n+1-j)k))P.  v(x)dT 
n=j  (n-1)k  k 


N 


+ B*  1 k(S*( (n+1-j  )k)  - E*n+1  J)P  v(nk) 
n=j  K k 
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Let  e*(t)  = T*(I-Pk)v(t)  = T*A*a(l-Pk)A*  °v(t) 

Since  Le  [ 0 , T ; U ] -*■  q[0,T ;H] ) for  p = 1 , ® (see  Theorem  2.2),  it 

follows  from  Lemma  2.3  that 

(6.4.1 ) T*  e J>  (Hp(1_q)[0,T;H]  - L [0,T;U])  for  p = 1 
Also,  we  recall  from  section  4.2  that 

|(I-P  )x|  < ck^jxj  for  p = 1,  ®,  0 < 6 < 1. 

L [0 ,T; H]  hJ[0,T;H] 

P F 

Again  by  Lemma  2.3  (and  by  self adjointness  of  Pk),  we  have 

(6.4.2)  | (I-P  )x | < ck3|x|  for  p = 1 , «,  0 S 6 £ 1 . 

HpP [0 ,T; H]  Lp[0  ,T;H] 

By  taking  @ = 1-q-a  in  (6.4.2),  x = A av(t),  and  by  combining  (6.4.1) 
with  (6.4.2),  it  follows  that 


|T*A*a(I-P  )A*'°v|  = |e*  I 

Lp[0,T;U]  Lp[0,T;U] 

S ck^  q a|v| 

L [0 ,T; D( A*_a) ] 

for  p = 1 , «,  0 < a < 1-q.  p 


Since  A ./?(U+H),  it  follows  from  Lemma  2.3  that 


(6.4.3)  B*A*  Q e £ (H-HJ ) 
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We  shall  be  using  this  fact  repeatedly  throughout  the  remainder  of  this 
proof  and  in  subsequent  proofs. 


t 

Let  e*(t)  = B*  j S*  ( t- ( j-1  )k  )P.v  ( Odt . By  (6.4.3.), 

(j-1)k  k 

t 

je*(t)  | = |B*A*'q  j A*q+aS* ( T“ ( j ~1 )k)A*"aP  v(x)dT  | 

(j-1  )k  k U 


< c J |A*q+aS*  ( t~  ( j - 1 )k)A*~aP  v(t)  I dr 

(j-1  )k  k H 


jk 


^ C j |A*q  + V(T-(j-1)k)A*'Y  v(T)|dT 

(j-1  )k  k H 


jk 


* c J 


1 . *- 


(j-1 ) k (x-(j-1)kq+“ 


|A  aP.v(T)|dT 


jk 

^ c (J 


(j-1  )k  ( T- ( j-1 )kq+a 


di)  |a  aPkv(jk)| 


< ck 1 -q_a  | A*_C 


V(jk)lH 


For  p = 1 we  have 


/ le2Ct) lu  = ^ i |e*(t)|dt 

0 j-1  (j-1  )k  ^ 
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For  p = 

ess  sup 
0<t<T 


< l k-ck1  q “|A  Vv(jk)  | 
j = 1 k H 


< ck 


N 

~q_a*k  I 
j=1 


|Pkv(jk}  li  * 
D ( A 


-a 


) 


£ ck 


N 

; l 

j = 1 


jk 

1/  V ( T ) d T | * 

(j-1  )k  D ( A ) 


.1  Q Ct  | n 
ck  |vj 


L1 [0 ,T; D( A 


l)] 


, we  have  similarly 


|e* ( t ) |y  = max 
j = 1 , — ,N 


{ess  sup  le*(t)|u} 
(j-1 )k<t<jk 


^ max  ck1  q a|A*  °‘p(<v(jk)  |H 

j = 1 , . . . ,N 


max  |Pkv(jk)|  * 

^ u \ 
M ) 


£ ck 


1 -q-a 


|v| 


L [0,T;D(A  “)] 
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We 


£ 

now  turn  to  the  estimation  of  e^(t) , where 


e*(t)  = B*(S*(t-(j-1 )k-I)  J S*(i-t)P  v(t)dT 


Again  by  (6.4.3) 


|e*(t)  |y  = |B  A “aA  q+a(S*(t-(j-1)k)-I)  J S*(T-t)A*'aPkv(T)dT| 


k 'U 


* c|A#q+a(S*(t-(j-1)k)-I)  j S*(T-t)A*"aPkv(x)dT|H 


t- (j-1 )k  T 

< c|A*q+0t  j A*S*U)d£  J S*(T-t)A*”aP  v(x)di| 

t k 

0 L 


H 


t— (j-1  )k  T 

^ c / |A*q+a+£S*(^)|d5  /|A*1_eS*(T-t)A*"aP  v(x)|dx 

t 


s 0 I -s^77  d«  / : ' i-e  I»'"“v(t)|<it. 


0 Cq*“*E  t (,-t) 


S ok'  q “ E / — frjlA*  “pk»(T)|dT 


t (l-t) 


For  p = 1 , we  have 


T T 


/ le3 ( t ) lydt  £ ok  q J J — |a  aPkv(i)|dTdt 


0 t (T~t) 
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By  reversing  the  order  of  integration, 


/ie3(t>lu 


J4.  , .1 -q-a-e 

ut  £ ck 


T T 

J (/ 


*- r 


0 0 (x-t) 

T 

< c^~q'a~e  / xe|A*"aPkv(x)|dx 
0 


dt)  ||A  Pkv(x)  ld't 


S 0TE/2  k1-q_a-e  J|P  Y<T>|  d, 

0 D(A  ) 


^ . 1 -q-a-e « i 

= ck  Ivl  #_ 

L1 [0 ,T; D( A “)] 


For  p = ®,  we  have  similarly, 


ess  sup  |e*(t)  1^  £ ck1  q a Gsup 


/ fr^lA  ~aPkv(x)|dx 


0£t<T 


0<t<T 


t (x-t) 


< ck1  q “ e(Te  j dx)  sup  |P.  v ( x ) | 


° 0<x<T  D(A  } 


£ ck 


1 -q-a-e 


|v| 


*-r 


L [0,T;D(A  “)] 


Now  we  let 


ej(t)  = B * 1 
n=j 


nk 

/ (S* ( x— ( j — 1 )k )-S* ( (n+1 -j )k) )P  v( x)dx 

(n-1 )k  k 
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As  before,  by  (6. A. 3) , 


N(t'lu- 


- - N 

|B‘A*'Vq  + a I 
n-j 


nk 

/ (S*  ( t — ( j — 1 )k  )-S*  ( (n+1  -j  )k  )')A*_aP  v(x)dT  |It 

( n— 1 )k  k U 


< c|A 


* 


q+a 


N 

l 


n-J 


nk 

/ (S*(t-(J-1 )k  )-S*  ( (n+1  - j )k ) )A*_aP  v ( t )d-t  | 

(n-1 )k  K h 


In  order  to  estimate  le* (t)  | we  shall  need  to  use  interpolation.  Let 


iM  ^ 

l I A 'e(S*(i-(j-1 )k)-S*((n+1-j)k))A*~Yv(T)dT, 
n=j  (n-1 )k  K 


and  furthermore,  let  a(t)  = a(jk)  = aj  for  (j-1 )k  < t < jk. 
Then 

|eJCt)|u  S o|A‘q*“*£a(t)|H. 

so  that 


- clal  * 

T;U]  Lp[0 ,T; D(A  q+a+e)] 


for  p = 1 , «>. 


It  is  clearly  sufficient  to  estimate 
do  so  by  estimating 


|a| 


Lp[0 , T; D (A 


*q+a+e 


)] 


and  we  shall 


lal  and  f a | 

Lp[0,T;H]  Lp[0 ,T ; D ( A* ) ] 

Then  we  shall  make  use  of  a standard  interpolation  theorem  [12,  Vol . I, 
p.  27]  to  conclude  the  argument. 


We  have 
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N nk 


|a.|  = | l j A ~e(S*(T(j-1  )k)-S*((n+1-j)k))A  ”aP  v(t)(1t| 

J H n=  i k 

n J (n-1  )k 


nk 


< l j |a  _£:  (S*  ( t- ( j- 1 )k  )-S*  ( (n+ 1 - j )k ) )A  _aP.  v(  t )dx  | dx 
n-j  (n-1)k  k H 


nk 


(n+1-j  )k 


I J |A*"e(J  A*S*U)d5)A"  aP.  v(t)  |„dx 

n=j  (n-1)k  T-(j-1)k 


*-r 


nk  (n+1-j)k 


I J (J  |A*1_eS*(0|d0  |A^aP  v(T)  | HdT 

n=j  (n-1 )k  t— ( j — 1 )k  k H 


nk  (n+1-j)k 


l / (J  -jzi  dO|A*'aP  v(t)|  di 

n=j  (n-1  )k  (n-j)k  E,  k h 


nk 


£ I “[(  (n+1-j  )k)e  - ( ( n- j ) k ) G ] f |a  aP.v(nk)|  dx 
n-j  e (n-1 )k  k H 


N 

£ ck  l [ ( (n+1-j )k)’ 

n-j 


( (n-j ) k ) e ] | A aPkv(nk)|H 


Therefore , 


|a| 


L1 [0 ,T;H] 


k 1 ‘ajlH 

J-1  J 


N N 


^ ck  * l l [ ( (n+1  -j  )k ) e 
j = 1 n=j 


( (n-j )k ) e]  |a  "aPkv(nk)|H 
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By  reversing  the  order  of  summation, 


N n 


|a  | S ck  1 ( l [((n+1-j)k)G-((n-j)k)e])  |a  aP  v(nk)| 


L^O.TjH]  n = 1 j = 1 


< ck2  l (nk)G  |A  ~aP  v(nk)  | 


n=1 


< ckTe(k  1 |A*~aP  v(nk)|  ) 


n=1 


ck  |v  | 


L1 [0 ,T; D ( A “)] 


Similarly,  for  p = °°, 


|a|  = max  |a  | 

L-[0'IiH]  j., N ' 


N 


< ck  max  I [( (n+1-j )k)G-( (n-j )k)G] |A*~aP  v(nk) | 

n=  ■ i K H 


j=1 N 


n=J 


< ck  max  |a  aP  v(nk)|  l [ ( (n+1 -j )k)G- ( (n-j )k) G] 


j— 


n=1 


< ck(Te) |P  v | 


L [0 ,T; D(A  ~a) ] 


ck 


'l  [0 ,T ; D ( A*_a) ] 


The  next  step  is  to  estimate  |a.| 


D (A* ) 


. We  have 
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l*ji 

3 D(A*) 


■ I* 


nk 


-III  * 

n=j  (n- 1 )k 


* i 


I -E 


(S*  ( t - ( j -1  )k )-S* ( (n+1- j )k ) )A  "aP  v(  t)c!t  I 

k 'H 


S 1 S jA*1  e(S*(T-(j-1)k)-S*((n+1-j)k))A*"aP,  v(T)Ldx 
n=j  (n-1 )k  k H 


We  estimate  this  expression  crudely  by  using  the  triangle  inequality 


N 

Ja  . | £ 1 

3 D(A*)  n=j 


| |A*1_CS*(T-(j-1)k| 

(n-1  )k  H*H 


+ I J |A  1 eS*((n+1-j)k)  ( d t | A * UP  v(nk ) | 
n=j  (n-1 )k  n H K H 


N 

< l j ( 2 

n-j  (n-1 )k  (t-(j-l)k) 


( (n+1-j )k) 


T=7)dTlA*’aPkv(nk)|H 


N 

^ c l [( (n+1-j )k)e 

n-j 


- ( (n-j )k)  + 


( (n+1  - j ) k) 


— ]|A*  aPkv(nk)jH 


For  p = 1 , we  have 


M = X l |a  . | 

L1[0,T;D(A»)]  j=1  J D(A«) 


N N 

-ck  I l C(  (n+1-j  )k)E-(  (n-j  )k)e+  * — ] |A*'aP  v(nk)  I 

J=1  n-j  ((n+1-j)k)1'e  k 'H 
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Sck  1(1  [ ( (n+1  - j )k)c-(  (n-j  )k  )E+ — ])|a*  aP  v(nk)| 

n=1  j=1  ( (n+1-j )k)1_e  H 

N n 

£ ck1+e  l (ne+  1 + j ■ dx)  | A aP  v(nk)| 
n=1  1 x1  G k H 


N 

£ c(Nk)e  l k |P  v(nk ) | 

n=1  D(A  a) 


= c lvl  *_ 

L1[0,T;D(A  “)] 


Similarly,  for  p = », 


la I = max  |a  | 

L.[OfT;D(A»)]  J D( A* ) 

J I | ) IN 


2 o max  j [ ( (rn-1-j )fi)  -( (n-j )fi)  * S — iJa01  P„v(nk)| 

J-1 N "‘J  C (n+1-j  )k)  ~e  k H 

#_  N 

£ ck  max  |a  "aPkv(jk)|  l ( (n+1-j  )£-(nk)E+ 1— — -) 

.1-1 N n=1  E 


< c(Nk)e|p  v| 

L^CO , T; D ( A “)] 


- clvl  « 

L^CO.TjDCA  a)] 


Let  us  summarize  the  results  we  have  just  obtained: 
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la  | * ck|v|  « . P - 1 

Lp[0,T;H]  Lp[0,T;D(A  “)] 


lai  s c |v | * , p = i,  ». 

Lp[0 ,T; D( A* ) ] Lp[0,T;D(A  “)] 

By  applying  the  interpolation  theorem  to  each  pair  of  estimates,  we  get 

1 _ g 

M * ok  |v|  , p = i,-t  o < e < i. 

Lp[0,T;D(A  D)]  Lp[0,T;D(A  )] 

We  are  concerned  with  estimates  in  Lp[0,T; D(A*q+a+e)] , hence  we 
take  S = q + a + e to  obtain 


|ej| 


Lp[0,T;U] 


la  I 


Lp[0 ,T ; D (A 


*q+a+e 


< ck1-q‘a_e 


)] 


|v| 


Lp[0,T ;D(A  a)j 


for  p = 1 , oo. 


This  completes  the  estimation  of  ejj . 

Finally,  we  let 
N 

e*(t)  = kB*  l (S*( (n+1-j )k)  - E n+1  J)P  v(nk) 

n-j  k k 

*_  N 

= kB*A  'Q  l A*q+a(S* ( (n+1 -j )k )-  E*n+1 ~j )A*~aP,  v (nk ) 

n=j  k k 

By  using  (6. A. 3)  as  before, 
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|e* (t) ^ S ck  l |A*q+a(S»((n+1-j)k)-E*n+1  j )A*~aPRv(nk) |H 
n=j 


In  order  to  derive  an  estimate  for  |e*f,  we  shall  use  interpolation  once 
again.  We  have  by  Lemma  4.1: 


„*n+1 -j  *-a 


*-  ■ 


j (S*  ( (n+1  - j )k  )-Er11’  ' J )A  uPkv(nk)|H  < |a"  aPRv(nK)|H 


and  also 


|A*(S*((n+1-j)k)-Ekn+1  j)A  aPkv(nk)|H  < -j  yR-  |A*  aPkv(nk)|H 


By  application  of  the  interpolation  theorem,  we  obtain 


|A*B(S*((n+1-j)k)-Ek^,-J)A""uPkv(nk)|H  < ^|A""°Pkv(nk>  |H 


,*n+1-j . *-a. 


ck  B . “-a. 


By  taking  6 = q+a,  we  get 


^ ok“q_c' 


le5Ct>lu  s ok  Z HTirrl**"\v("k)li 

n=j 


^ ck 


For  p = 1 , 


1-q-a-e  r , e 1 ■ „ *~a„  , , . , 

Jj  k • snrj|A  V(nk)l„ 


/|e«(t)|  dt  - k l |e«(jk)| 
n j = 1 b U 


£ ck*k 


. N N , 

1-q-a-e  ^ v ,.e  1 i.*~ar 


‘J,  Jj k h^i|# " pkv<nk)iH 


£ ck 


, N n 

•k,_<r“_E  i ( i 


n 


= 1 j-1 


H7PJ)jA  PRv(nk ) JH 


< c(kGln  N)k1_q  a"e  l k|p  v (nk)  | , 

n=1  D ( A a) 


ck 


1 -q-a-e 


v I 


L1 [0 ,T; D(A 


l)] 


For  p = °°,  we  have 

ess  sup  je*(t)j  = max  |e»(jk)ju 
0St<T  j-1,...,N 

- ck1  q “ e max  £ ke^-l:-T|A*"aPkv(nk)  | 

J-1 N n=j 

- ck1  q a e(keln  N)  max  |a  aPRv(jk)|H 

j-1 N 


ck1  -q-a_e'|v  j 


L [0 ,T; D(A  “)] 


The  proof  of  Theorem  A. 2 is  now  completed,  since 


(T*  Tk)v(t)  - e^t)  + e2(t)  + e^(t)  + e*(t)  + e^(t). 
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6.5  Proof  of  Theorem  5.1 


Let  t € [0,T]  and  let  j be  the  unique  integer  such  that 
( j — 1 ) k < t S jk.  We  recall  that  by  definition 

(6.5.1)  T*v(t)  = Tjphv(t) 
and 

(6.5.2)  T*kv(t)  = (ThkPhPkv)(t)  = (T*kPhPkv)(jk)  for  te[0,T]. 

We  shall  use  the  following  decomposition  for  (T*-T^k)v(t) : 

(T*~Thk)v(t)  = (T«-Tj)v(t)  + (VThk)v(t) 

= (T*-Tj)v(t)  + Tjv(t)  - Tjv((J-1)k) 

+ Tjjv((j-1  )k ) - T^kv(jk). 

Remark.  The  apparent  inconsistency  in  the  expression  just  above 

(jk  vs.  ( j - 1 ) k ) is  purely  a matter  of  notation.  As  it  turns  out, 

* 

Thkv(jk)  does  depend  on  the  values  of  v over  the  interval  ( (j-1 )k,T] 
(see  the  definition  of  PK  in  section  4.2). 

// 

In  the  course  of  the  proof  we  will  make  crucial  use  of  the  inverse 
approximation  property  (3-4),  which  we  recall  here: 
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(6.5.3)  lB*vh I u 5 ch_2qlvhiH 


Wh€  Vh. 


for  some  q < q. 

Let  e^t)  = (T*-T^)v(t).  By  Theorem  3.1,  the  following  estimates 

hold 


(6.5.4)  je.  J < ch2^1  q^ln  r|v|  for  p = 1 ,®. 

Lp[0,T;U]  Lp[0,T;H] 

Now  let  e2(t)  = T*v(t)  - T^v((j-1)k).  As  we  have  noted  in  (6.5.1), 
we  can  assume  without  loss  of  generality  that  v(t)  = Phv(t)  = v^(t)  for 
t [0  ,T] . Then 


eo(t)  = B*[  j S*(T-t)v.  (T)dT  - J S*(t-(J-1 )k)v.  (T)dT] 
2 t h h (j-1)k  h h 


= -B*(sJ(t-(J-1)k)-I)  J S*(T-t)vh(T)di 

t 


•B*  J 


(j-1  )k 


Sh ( t— ( j — 1 )k)vh(t)dT 


By  (6.5.3)  and  by  differentiability  of  S^(t),  we  have 


|e2(t)|u  * |B* (S*(t-(J-1 )k)-l)/s*(T-t)vh(T)di|l 


|B*  j S* ( t- ( j -1 )k ) v ( x )di  | 

(j-1  )k  n h U 


t— ( j — 1 )k  T 

s oh'2ql  I Ahsh<5)d«  | sh(’-t)''h<T)dIl„ 
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+ ch~^l  / S*(x-(j-1  )k)v  (x)dT  |, 


(j-1  )k 


Let  p = °°.  We  need  to  consider  two  cases. 


(i)  t + k < T 


(ii)  t + k > T. 


In  both  cases,  we  shall  make  use  of  the  commutativity  of  Ah  and  S 
We  also  fix  a,  with  0 < a < In  case  (i), 


|e*(t)  | < ch  2q|J 


t- ( J-1 )k  t+k 

*r.  * r *1  -a  * 


II  W«)d«  / #h  '“4 Sh<T-t)Vh<T)dTlH 

0 t 


t-  (J-1  )k  T 

• J S*(5)dE  / Vh(r-»>vh(t>H 


+ ch"2q|  / S*(T-(j-1)k)v  (t)<1t|h 

(j-1  )k  n h H 


_ t+k 

S oh'2q  J |A*1_aS*(T-t)|H<„d,|vh| 


L [0,T;H] 


k T 

+ oh’2q  |lsh(E„H-Hde  [o  3 


4 0t'"2<1  / |S*(T-(j-l)k)|  d,|v  | 

(J-1 )k  h H4H  h L [0,T;H] 


zr  * 
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By  analyticity  of  S^(t),  we  have 


t+k 


|e*(t)|,,  < ch“2q  j l-d^  \ 


1 


o C 


t (x-t) 


7 dljv  I 

1 -ct  h 1 


L [0  ,T; H ] 


_ k T 

+ ch_2q  I ^ i 7~jrdTlvhl 

0 t + k L [0  ,T;  H ] 

00 

_ Jk 

+ ch  2q  f dxjv  | 

(j-1  )k  n L^CO ,T;H] 


< ch 


-2a 


1 - a a » « 

k k jvh| 


ch“2qk 


In 


L [0 ,T; H] 


rlvhl 


L [0,T;H] 


+ ch_2q  k|vh| 


L [0,T;H] 


< ch  2qk  in  l|vh| 


L [0,T;H] 


As  for  case  (ii). 


t- (j-1 )k 

-2q,  r *a  * 


|e2  ( t ) |u  < ch_2q  | J AhaShU)d5  j Ah1_aSh(T-t)vh(x)dT|H 

0 t 


+ ch  2q|  J S.  ( t-  ( j-1  ) k )v,  ( x )dx  |H 
(j-1  )k  h h »H 


k T 

-2q  ^ fi.*1-a„*, 


s oh-  /|Ah\(5)|H<Ho«  /l*n  [0iT;HJ 


jk 


+ ch  2q  J |S*(x-(j-1)k)|  dx|v  | 

(j-1  ) k h H H ^ L [0,T;H] 
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_ k T 

< Ch"2q  j l_dt  j — L_ _dTjv  | 

0 f t ( T-t  ) 1 “ 


L [0 ,T;H] 


_ Jk 

+ ch"2q  j dt|v  | 

(j -1  )k  LjO.TjH] 


ch~2q  k1  a(T-t)ajv  j 


+ ch 


■2qt 


L [0,T;H] 


kKi 


L [0,T;H] 


< ch  2q  k|v  | 

LjO.TjH] 


Thus  we  have  proved 

ie*l  * ch_2qk  In  £|v  | 

L [0,T;U]  K L [0,T;H] 

00  CO 


£ ch~2qk  In  i|v| 

LJO.TjH] 

For  p = 1 , we  have 

T T T 

J|e*(t)|u  dt  S J|B*(S*(t-(j-1  )k)-I) / S*(T-t)vh(T)dx|lJdt 

T t 

+ JlB*  J Sh(x-(j-1 )k)v  (i)dx  | dt 
0 ( j - 1 ) k n n U 

_ T T 

< ch"2q  J / |(S»(t-(J-1 )k)-I)S*(T-t)vh(x) |Hdxdt 
0 t 

_ T t 

+ ch  2q  J J |s*(i-(j-1 )k)v  (x)  | didt 

0 (j-1)k  n h H 
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By  reversing  the  order  of  integration, 


T 


T t t-(j-1  )k 


j|e*(t)|dt  S ch”2q  j l\(j 


A*S*U)dOS*(i-t)v  (t)j  dtdr 

0 2 U 0 0 0 h h h h H 


T nk 


+ ch  2q  / / |S*(t-(j-1 )k)vh(T)  |HdtdT 


0 T 


where  n is  the  unique  integer  such  that  (n-1 )k  < t i nk.  Furthermore, 


T _ k t t-(j-1  )k 

J |e*(t)  | dt  < ch_2q  j / 1 (/  A*S*(OdOS*(T-t)v  (T)LdtdT 
0 0 0 0 


_ T T-k  t- ( j- 1 )k 

*o h-2MJ  kj  A*S*(C)d^)S*(t-t)vh(T) |HdtdT 

k 0 0 

T t t-(j-1)k 

♦ =»'2qn  K/  A»S*(^)d^)S*(T-t)vh(T)|HdtdT 

k t-k  0 


T nk 

ch  2q  / ( / dt)  |S*( t- (j-1 )k)vh( t)  |Rd t 

0 T 


k t k 

< oh_2q  J /(J|A.“SJ(C)|H<Hd5)|»*,-'IS.(T-t)|H<Hdt|vh(,)|Hd, 


T t-k  k 


* 0h-2q  J J (J|S.(5)|H<Hd5)|A.S.(T-t)|H<Hdt|vh(,)|HdT 


_ T T k 

* oh"2q  l 1 i„.„“iy  <>  In*' 

K T~K  U 
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* Oh'2,5k  /|S»(T-(J-1)k)|H<H|vh(,)|Hd, 


k t 

s oh'2q  j k ’■“(/  — 5-^  dt)|»h(T)|HdT 


0 0 (T-t) 


_ T x-k 

» dh'2q  / k(J  -i-  dt)|v  (,)|  d, 

k 0 T C n H 


_ T T 

+ ch_2q  I K1  ”“(/  dt)|v  (t)|  dT 

k T-k  (T-t)  a 


T 

+ ch  2qk  J|vh(t)  lndT 


S oh'2qk|v  | * oh'2qk  in  l|v  | 

L1 [0,T;H]  K n L1 [0 ,T;H] 


ch"2qk|vhj 


L1 [0fT;H] 


+ ch  2qk|vhJ 


L1 [0,T;H] 


^ ch  2qk  In  f|v  I 

k 1 h 1 


L.j  [0,T;H] 


Thus  we  have  proved 


b*| 


L<1  [0,T;U] 


ch  2qk  In  j- 1 v,  I 

k1  h1 


L1 [0  ,T; H] 
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We  now  undertake  the  estimation  of  e^(t)  « Thv((j-1)k)  - Thkv(jk),  which 
we  decompose  in  the  following  way: 


N 


e3 ( t ) = B*[J  S*(x-(j-l  )k)Phv(T)dx  - k j E*k+1  Jph^Pkv)  (nk)] 


(j-1  )k 
nk 


n=J 


N 


= B*  l [J  S*(t-(J-1  )k)P,  v(x)di  - kE*"+1  JP.(P,  v)(nk)]. 

n=j  (n-1 )k  h h hk  h k 


As  in  the  estimation  of  e2(t),  we  shall  assume  without  loss  of 
generality  that  v(t)  = Phv(t)  = vh(t).  We  also  make  the  following 
observation,  based  on  the  definition  of  Pk: 


nk 


nk 


j S*((n+1-j)k)vh(T)dT  = S*((n+1-j  )k)  / v^ijdx 


(n-1  )k 
(6.5.5) 


(n-1  )k 

kS*((n+1-J)k)(Pkvh)(nk),  n=1 ,2,. . . 


Hence  we  have  the  further  decomposition 


nk 


nk 


* / . 

e3(t) 


B*  l [/  S*(t-(j-1 )k)v  (x)dT-J  S*( (n+1-j )k)v  (i)dt] 

n=j  (n-1 )k  h h (n-1 )k  h h 


+ B*  I [kS*((n+1-j)k)(Pkvh)(nk)  - kE*"+1  J (P„vJ  (nk) ] 


n=J 


hk  k h' 


N 


nk 


= B*  I / (S*(T-(j-1)k)  - S*( (n+1-j )k) )v  (T)dT 

n=j  (n-1 )k  n h h 


* 

B*  l k(S*((n+1-j)k)  - Ehk+1  J)(Pkvh)(nk) 
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By  application  of  (6.5.3),  we  have 


- N 


nk 


|e*(t)|  S ch"2q  l \ j (S*(T-(J-1)k)-S»((n+1-j)k))v  (t)c!t|h 


n=j  (n- 1 ) k 


n * i _ • 

+ ch  l k|(S*((n+1-j)k)-Eh"+1'J)(Pkvh)(nk)|H 


n=J 


Jk 


■ oh‘2q(J_1)k  (|sh(T-(J-,)k,U  * 


(t) 


-2q  N 


nk 


+ °h  " n ' + l{  1)J(Sh(T"(j"1)k)"Sh((n+1"J)k)JH^Ivh(T)lHdT 


- N 


* oh’2qkn^  |ss(  (n41  )k)-<r 1 _J  u i<pkvh  > <"k)  i„ 


jk 


* ch"  qk(i  J |v  (T)|  dT) 
(j-1)k  n H 


(t) 


- N 


nk  (n+1-j)k 


ch"  Q I / 1/  AhShU)d^H+HlVh(T)lHdT 

n=j  + 1 (n-1 )k  t— ( j — 1 )k  n n H H n H 


N 


* 0l,'2q,<  Jj|S5((n*1-j)l!)-E*^'-J|H<H|(P1(»h)(nk)|H 


( t)  This  term  is  absent  if  j=N.  (Same  convention  on  next  page). 
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By  definition  of  Pk,  by  analyticity  of  Sh(t),  and  by  Lemma  4.1 


follows  that 


jk 


|e  (t)|  < ch'  q k(-  / lvh(x)lHdT) 

d U K( j-1  )k  n H 


(t) 


-2q  N 


nk  (n+1-j)k 


+ ch  q l j (j  -dZ)  |vh(x)  |H d t 


n=j+l(n-1)k  (n-j)k 


* ch’2<1  k L l(pkvh)(nk)  Ih 

n=J 


jk 


ch"  Q k(l  j |v  (t)|  dT) 
(j-1)k 


(t) 


nk 


+ ch"  Q l lvh(T>lHdT 

n=j  + 1 ^ (n-1)k 


+ ch 


■2q 


N 

I 

n=j 


1 


n+1-j 


l(PkVh)<nk)lH 


For  p = 


we  have,  since  e^Ct)  = 


e3(jk) , 


\<\ 


L [O.TiU] 


= max  |e*( jk)  < ch  2qk |v^ | 


1 -j-N 


L [0,T;H] 


+ ch  2qk  In  N |v  | 


L [0 ,T; H] 


+ ch  2qk  In  N |?kvh  | 


Thus  we  have  proved 


L [0,T;H] 


1-3 1 


L [0,T;U] 


ch  2qk  In 


SM 


L [0,T;H] 


/ 


, it 
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< ch  2qk  In  ~ | v | 


L [0  , T ; H ] 


Similarly,  for  p = 1 , we  have 


\A\ 


N 

= k l |e*(jk)  | < ch 

[0,T;U]  j=1  5 


2q 


N Jk 

X l I 

j-1  (j-i  )k 


lvh<T 


) Ih0t 


- N-1  N 

,"2q 


nk 


+ ch“2qk  l l In  |v  (T)|  dT 

j-1  n-J  + 1 J (n-1  )k  n ri 


i,HTei<Vh)!nk)l„ 

J-1  n«j 


By  reversing  the  order  of  summation, 


[0,T;U] 


ch  2qk|vh| 


[0,T;H] 


N n-1 


nk 


+ ch"2qk  1 ( l ln  J |v  (Oj  d-r 

n=2  j=1  nj  (n-1 )k  n H 


— „ N n 


+ HTT=jH(Pkvh)(nk)|H 


< ch~2qk|v  | 

L1  [0,T;H] 


+ ch 


-2q 


N 

k I 


n=2 


nk 

inn/  lvh( t)  lHdT 

(n-1 )k  n H 
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-2q,  2 
ch 


N 

l 

n=1 


In  n|(Pkvh)(nk) |H 


Ch‘2qk|»h| 


ch  2qk 


L1 [0 ,T;H] 


In  N |v  | 

L1 [0,T;H] 


+ ch  ^qk  In  N(k  l l(Pkvh)(nk)iH) 
n=1 


ch  2qk|v 


I 

L1  [0  , T; H] 


ch  2qk  In 


k , Vh  I 


L1 [0,T;H] 


Thus  we  have  proved 


le3 1 


L^O.TjU] 


£ ch  2qk  In  ^|v,  I 

k*  h1 


L1 [0,T;H] 


£ ch  2qk  In  -r-|v | 

L1 [0 ,T; H] 

Since  (T*-T^)v(t)  = e^(t)  + e2(t)  + e^(t),  this  concludes  the  proof  of 
Theorem  5.1 . 


6.6  Proof  of  Theorem  5.2 


4.2. 


Before  starting  with  the  proof,  we  give  a corollary  of  Theorem 
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Corollary  6.3.  Let  T*  be  defined  as  in  (2.37),  and  let  T*  be  defined  by 


T*v ( t ) = T*v(jk)  = B* 
Then 


N 

l I 

n=j (n-1  )k 


S*  ( (n+1  - j )k ) (Pkv ) ( x)dt . 


(6.6.1)  | (T*-T*)v  | < ck1  q e|v| 

L [0 ,T; U ) LjO,T;H] 


p = 1 , «. 

// 


Proof  of  Corollary  6.3.  The  result  follows  immediately  from  the 
observation  that  (T*-T*)v(t)  = e*(t)  + e*(t)  + e^(t)  + e*(t),  with  e*(t) 
as  defined  in  the  proof  of  Theorem  4.2.  In  particular,  by 
taking  a = 0 in  Theorem  4.2  one  obtains  (6.6.1). 


Coming  back  to  the  proof  of  Theorem  5.2,  we  set 
(T  _Thl<)v(t)  = e*(t)  + e^(t)  + e^(t),  where 


// 


e*(t)  = (T*  - T*)v(t) , 


e*(t)  = (T*  - T* )v ( t ) , 


e3(t)  = ^ ~ TAk)v(t)’ 


nk 


with  T*v(t)  = B * 1 j S*((n+1-j)k)P.  (P.  v)(T)di. 

n=j (n-1 )k  n n k 


The  estimation  of  e^ (t)  is  precisely  the  conclusion  of  Corollary  6.3, 
i.e.,  (6.6.1).  The  estimation  of  the  other  two  terms  e^(t)  and  e^(t) 
will  require  the  following  result. 
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Leiruna  6.4.  For  0 < h < 1,  let  {V^}  be  a family  of  subspaces  of  D(A*) 
for  which  properties  (3-1)  and  (3.^')  hold,  and  let  { A^}  be  a family  of 
approximations  of  A such  that  properties  (3.5),  (3.6)  are  satisfied  for 
both  A and  A*.  Then 


A*A^  ]Ph  e £ (H+H) 


// 


Proof  of  Lemma  6.H.  For  simplicity  let  vh  = 


Phv  for  veH.  We  have 


^A*Ah  Vh " lAh  Vh ^D( A* ) 


" H(PhA  Ah  Ph)vh^D(A*)  + iPhA  ViJd(A») 


Using  the  inverse  assumption  (3.^'),  we  derive 


l#*C1vhlH  S oh"2|(Ph#*’''Ah"lph)vhlH  * °lVhlH 


S o»'2[|(A*'1-A*',Ph)vh|H  . |d-Ph)»*'’vh|H.h2ivh|H] 


We  now  use  (3.1)  and  (3.5)  to  obtain 


.*-1 


lA*Ah  \Ih  s =h'2C=l>2|vh|H.  oh2|vh|H  * h2|vh|H] 

s °KIh 


and  thus  Lemma  6.M  is  proved. 


// 
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Coming  back  to  the  proof  of  Theorem  5.2,  we  are  now  in  a position 
to  estimate  e£(t)  and  e^(t). 


From  the  observation  (6.5.5)  we  made  in  section  6.5,  we  have 


(T*  - T*)vCt)  = (T*  - T*)v(jk) 


nk 


= B*  l J [S*((n+l-j)k)-S*((n+1-j)k)P.  ](P  v)(T)dt 
n=j  (n-1  )k 


= B*  l k(S*( (n+1-j )k)-S*( (n+1-j )k)Ph) (Pkv) (nk) 


n=J 


Thus 

* 


I e2 ( t ) fu= Jb*A*  q l kA*q(S*( (n+1-j )k)-S*( (n+1-j )k)P  ) (P.  v) (nk)  L 
n=j 


S c I k|A*q(SM(n+1-j)k)-S*((n+1-j)k)Ph)|H+H|Pkv(nk)[H. 


We  recall  from  Lasiecka  [11]  that 


21 

(6.6.2)  |s*(t)-S*(t)Ph|H+H  < c , 0 < £ < 1,  t > 0. 

t 


By  Lemma  6.4  and  by  analyticity  of  S*(t)  and  S£(t),  we  also  have  that 


|A*(S«(t)-S.(t)Ph)|H<H  < |A»S«(t)|H<H  ♦ |A*A*_1|H,„l»J;S*<t)Ph| 


H+H 


(6.6.3) 


t 


t > 0. 
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We  now  apply  the  same  interpolation  theorem  [ 1 3 ] we  have  used  several 
times  before  to  (6.6.2)  and  (6.6.3)-  Thus  it  follows  that 


. 2£( 1-q) 

|A*q(S*(t)  - S*(t)P  ) | u < c , . — 

h h 'ri+H  t£(1-q) 


. — , 0 £ q < 1 


Hence  by  taking  £ = 1- 


1-q* 


, N 2(1 -q-e) 

2(t)|u  < c I k ^— | (pkv)  (nk)  iH  • 

n=j  ( (n+1-j )k) 


For  p = we  have 


ie2l  = max  |e*(jk)| 

L [0 ,T; U]  1 < j < N 


N 


Sch2<’-"->(  ik*.  -f^py-i 

n=1  n 


L [0,T;H] 


£ ch2(1_q"e)|v| 


L [0,T;H] 


For  p = 1 , we  have  similarly 


Ie2  I = k l le2(Jk)ln 

L*,  [0  ,T;U]  j = 1 


£ ckh 


2( 1 -q-e) 


N N 


1 lk  ’ -fr^l(Pkv)(nk)li 


j=1  n=j  (n+1 - j ) 
N n 


< ckh2(1_q_e)  l ( l ke.  — :— ) |(P  v)(nk)  | 

n=1  j=1  (n+1-j )'  e k H 


* ckh2(1_q"£)  l (nk)e|(P  v)(nk)l 
n=1  K H 


< ch2(1_q_e)(k  l |(Pkv)(nk)  I ) 
n=1 


< cn2(1“q_e) jv | 


L1 [0 ,T; H] 


The  estimation  of  &2  is  completed,  and  now  we  turn  to 


e3(t) 


(T*  - T*k)v(t) 


B*  I k(S*((nM-j)k)-E*k*1'J)Ph(Pk»)(nk) 


In  fact,  by  (6.5.3) , 


|e»(t)|u=|B*A*"q  I kA*q(S*((n+1-j)k)-E*U+1“J)Ph(Pkv)(nk)|u 


5 °n^l**,(s;((n*,-J)k)-E*r,'J>PhlH.Hl‘V»<^)lH' 


From  Lemma  4.1  applied  to  the  semigroup  S^(t),  we  have 


(6.6.1)  |(S*<<n*l-j)k)-E*"+,-J)PhL  „ S 


hk  h 'H+H  “ n+1-j 


and 


(6.6.5)  |A*(S*((n+1-j)k)-E*”+1  j)P  j 


hk  h 'H+H  ” (n+1 -j  )k ' 


Hence  by  Lemma  6.4, 


(6.6.6)  |A*(S*( (n+1 -j  )k  )-E*f*+1  J)P  I 

n nk  h 1 


H+H  (n+1 -j )k ' 
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Interpolation  between 


(6.6.4)  and  (6.6.6)  yields 


*n+1-jvD  ■ 


|A*q(S*((n+1-j)k)-Ehk  ;phlH+H  “ ’ kQ 


. 1-  , 0 S q £ 1 


Therefore, 


. , Y kc  1 rS(P,v)(nk)L  , e>0- 

|e|(t)  ly  s c L k (n+1  -j  ) ' K H 


n=J 


For  p = ”, 


|e«|  * .k’-O-t  ! k‘.  i)IVl 

3 L [0,T;U]  n=1  L_LU , 1 , n J 


s ek^M 


L „[0,T;H] 


Similarly,  for  p = 1 , 


* 3 L^O.TjU] 


|e;|  < ok1'q'Ek  I I.k‘nn=jl<V'(nK,lH 

i^i  ...  j=i  n=j 


N n 


< ck1^  I ( I ke^)|U  (nk)|H 

n=1  j=1 


< Ck’^^k  l |(Pkv>Cnk)lH 


n=1 


< ck1  q_E|v| 


L^O.T-.H] 
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Since  (T*-T^k)v(t)  = e*(t),  e^U)  + e*(t),  the  proof  of  Theorem  5.2  is 
complete . 


CHAPTER  VII 
NUMERICAL  EXAMPLES 


7.1  Opening  Remarks 


In  this  chapter  we  shall  implement  the  fully  discrete  approximation 
scheme  described  in  Chapter  V in  the  canonical  case  of  the  heat 
equation.  We  shall  also  present  numerical  results  from  a computer 
program  which  implements  this  discrete  scheme  in  the  one-dimensional 
case  with  P = (0,1).  These  results  will  be  presented  in  the  form  of 
tables  and  figures,  and  we  shall  moreover  explain  the  nature  and  scope 
of  these  results  so  as  to  facilitate  their  interpretation 


Consider  the  canonical  heat  equation  with  nonhomogeneous  Dirichlet 
boundary  conditions 


7.2  Implementation  of  the  Fully  Discrete 
Scheme  for  the  Heat  Equation 


(7.1  ) 


ut(t,x)  = Au(t,x) 


u(0,x)  = uQ(x) 


on  Q = (0  ,T)xP 


on  I = (0,T)xf 


on  0 
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Let  A - -A  + zero  Dirichlet  boundary  conditions.  Then  by  Green's 
formula,  the  bilinear  form  a(u,v)  associated  to  -A  is 

a(u,v)  = (Au.v)^  = - (Vu , Vv  )^  + <f^,v>r 

= -(Vu.Vv)  Vu,veH°(p). 

We  now  recall  from  Chapter  III  the  continuous-time  approximation 
algorithm: 


find  uh(t)  e Vh  such  that 


9vh 

(?uh(t),Svh)0  - <g(t).  5^->r. 


M 

If  ^i^i  = 1 f'0rm  a basis  for  the  subspace  Vh,  and  if  we  let 
M 

uh(t)  = I then  (7.2)  can  be  reformulated  as 


(7.2) 


(“h(t)'Vh)D 


uh(0)  given 


M 


9 ^ . 


t)(*i,«J)Q-  -.Iiai(t)(V0>.,V<0j)  - <g(t),g^>r,  j = 


= 1 M 


(7.3) 


a. (0)  given,  i = 1 , . . . ,M. 


Let  a<t>  ' (a1<U a„(t))T.  L-  «V*jVi.J-1 
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9 4>,  t 

K * ((V<0  V*  ) ) and  F = (-1 , . . . ,-H)1 , 

1 J w l »J  = 1 M av  9v 


Then  (7.3)  has  the 


following  compact  form: 


f La(t)  = -Ka(t)  - <g(t),F>r 

(7.4)  J 

J a(0)  given 

We  are  dealing  here  with  a linear  system  of  ordinary  differential 
equations.  Thus  the  discrete- time  approximating  scheme  described  in 
section  5.2  can  be  applied  to  (7.4)  with  minimal  modifications.  If  we 
set  gn  = g(nk),  where  k is  the  parameter  of  discretization  in  time,  and 
let  an  be  the  discrete-time  approximation  of  a(nk),  then  the  fully 
discrete  approximation  of  (7.1)  takes  the  form 


k L(an+l"an)  Kan+1~<gn+1 ,F>r  n 0 » 1 » * - * 

a = a(0)  given 

0 

or  by  factoring  an+1  , 

[ (L+k*)an+1  = Lan  - k<gn+1,F>r,  n-0,1.... 

(7.5)  / 

1 aQ  = a(0)  given 

( 

As  we  pointed  out  in  Chapter  V,  we  are  now  faced  with  a sequence  of 
linear  systems  of  algebraic  equations.  This  is  theoretically  a simple 
task,  but  in  practice  it  can  be  quite  challenging  for  very  large  values 
of  M.  The  difficulty  also  increases  with  the  condition  number  of 
L + kK,  which  itself  is  usually  difficult  to  evaluate. 
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7.3  Numerical  Computations  and  Results 


In  this  section  we  shall  calculate  the  matrices  L and  K as  well  as 
the  vector  F in  the  case  of  the  one-dimensional  heat  equation 
with  q = (0,1),  and  where  the  approximating  functions  are  linear 
splines . 

More  precisely,  if  the  unit  interval  is  partitioned  into  M = 1 /h 
equal  subintervals,  then  the  "hat"  functions  defined  by 


(x) 


Jx  - (i-1)  if  (i-1 )h  < x < ih 

x + (i-1)  if  ih  £ x < (i  + 1)h,  i = 1,...,M-1 

0 otherwise 


span  a space  of  linear  splines  on  [0,1]  of  dimension  M-1  . We  note  that 
these  splines  satisfy  zero  Dirichlet  boundary  conditions. 

It  can  be  readily  seen  that 

(<f>i  » 4>i_i  )n  = ^ h for  i = 2 M-1 


and  that 


**i  ^ 


for  i = 1 , . . . ,M-1 


Also  = ° 

Furthermore 


(V*.  , V <j> 


i-1 


h* 


for  | i-j  | >2 


i = 2 M-1 


(vVvVq 


2 

h 


and 


i = 1 


9 •••  9 


M-1 
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anu  also 


(V4*i  )Q  = 0 


if  i-J  * 2. 


Finally,  we  have 


. I h 


3v 


if  i-1  ,M-1 


if  i*1  ,M-1 


Remark . We  can  see  that  the  matrices  L and  K are  both 
tridiagonal.  This  property  is  very  advantageous  from  a computational 
point  of  view. 

// 

Thus,  after  simplification,  the  fully  discrete  approximation  (7.5)  takes 
the  form 


(7.6) 


S' a , = $>a  + cG , n = 0 , 1 , . . . 

n+1  n 


aQ  given 


where  o - 6-g.  - * * 2c,  - 1^,.,  - 1-c. 


*1.1  - "•  *1.1-1  • *1-1,1  ■ ’•  G(,)  ■ «„.1<0)’  <3(M-’)  ■ Vi(1)’ 


and  all  other  entries  are  equal  to  zero. 

We  shall  now  present  the  results  of  a computer  program  which  solves 
the  sequence  of  linear  systems  (7.6)  for  given  input  data  aQ 
and  { gn } such  that  the  formula  for  the  corresponding  solution  of  problem 
(7.1)  is  known.  The  program  then  calculates  the  relative  error  of 
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discretization  ehk  in  the  norms  of  [0 ,T; (£3)1  (resp.  L^CO.TjL^fl)] ) 
for  g in  L [0,T; L_( f ) ] (resp.  L [0 ,T; L„( D ] ) t i.e., 

i 2 00  2 

\u~\k\ 

L [0,T;L2(Q)] 

ehk  = — > P = 1 * 00 - 

Lp[0,T;L2(D] 

The  values  of  ehk  as  a function  of  h and  k are  calculated  for  two 
particular  solutions,  say  u-|  and  u2,  where 

2_ 

u^t.x)  = x(x-1  )[X- — + t]  + t2  , t > 0,  0 £ x < 1, 

and 

2 

u2(t,x)  = eX+t^  T ^cos  Tr(x+2t)  , t > 0,  0 i x i 1, 

In  tables  1 and  2,  the  values  of  ehk  corresponding  to  u1  are  shown 
in  both  the  cases  p = 1 and  p = <*>  respectively;  the  same  is  true  for 
tables  3 and  where  this  time  ehk  corresponds  to  u2.  More  precisely, 
the  values  of  log-|g(10  ehk)  are  displayed  in  terms  of  log2  - and 
l°g2  which  increase  by  1 whenever  h and  k are  divided  by  2. 

If  we  let  k = ch^  for  some  c > 0,  then  h ^2  + k ^ = h ^2 (i  + e ^ ) . 
So 

log10(h1/2+  k /j4  ) = log1Q  h1/2+  log1Q  (1  + c /i4  ) 

= ~ V2  log-]  o £ + c 
= - V2  .(log102)(log2  + c 
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Thus  the  rate  of  decrease  of  log.]Q  (h  ^2  + k ^ ) with  respect  to 

log?  ^ is  constant.  Moreover,  we  can  see  that  the  observed  rate  of 

-0.15  (asymptotic  rate)  is  very  close  to  the  theoretical  value  of 

-l/^log-iQ  2 - - 0.1505.  Our  theoretical  rates  of  convergence  are 

therefore  confirmed  by  these  numerical  examples. 

The  information  contained  in  tables  1 — is  displayed  in  graphic 

fashion  in  figures  1-4.  These  figures  provide  a more  intuitive  and 

qualitative  look  at  the  numerical  results.  For  example,  the  figures 

2 

clearly  indicate  that  the  theoretical  restriction  k < ch  is  also 
necessary  in  practice  in  order  to  obtain  the  best  rate  of  convergence. 

It  is  also  apparent  from  figures  3 and  4 that  the  strategy  of  fixing  one 
parameter  of  discretization  while  decreasing  the  other  is  not  only 

suboptimal,  but  can  also  cause  serious  instability  problems. 

The  theoretical  rates  of  convergence  derived  in  Theorems  5.3  and 

5.4  are  valid  for  "rough"  boundary  input,  i.e.,  for  input  in 

Lp[0,T;L2(r)] . In  fact,  our  last  example  shows  calculated  values  of  the 

solution  in  the  case  of  discontinuous  boundary  input.  These  values  are 

calculated  for  two  choices  of  the  boundary  input  as  the  mesh  size  h = ^ 

M 

is  decreased  while  respecting  the  condition  k = ch  . As  can  be  seen  in 
Table  5,  the  calculated  values  of  both  solutions  exhibit  a clearly 
convergent  behavior.  This  confirms  the  stability  of  our  scheme  (5.7). 

Remark . The  solutions  u1  and  u2  restricted  to  the  boundary  are 
clearly  infinitely  differentiable.  It  is  also  clear,  however,  that  the 
rates  of  convergence  are  not  improved  although  one  might  expect  this  to 
happen.  An  alternative  to  this  situation  may  be  provided  by  an 
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extension  of  Nitscne's  method.  In  fact,  in  the  semidiscrete  case  it  is 
known  that  Nitsche's  method  [17]  has  rates  of  convergence  which  improve 
with  the  smoothness  of  the  boundary  data. 

// 
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